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PREFACE 


In days gone by, hypervelocity flight in planetary atmospheres primarily concerned entry into the 
atmosphere and descent to the ground. But in recent times, the discipline has been broadened in two 
directions. The first concerns reusable flight configurations that are primarily space vehicles which use 
the upper layers of an atmosphere to perform hypervelocity aerodynamic maneuvers to change orbital 
altitude, orbital plane inclination angle, aerocapture on return from another planet, or rendezvous with 
another space vehicle. The second direction involves vehicle concepts that fly from the ground into 
orbit by use of air-breathing powered lift. The former is often referred to as an aeroassisted orbital (or 
space) transfer vehicle (AOTV) (or ASTV) which skips out of the atmosphere and returns to space after 
the aeromaneuver that was performed during the shallow dip into the upper atmosphere. These vehicles 
can greatly enhance space payloads because they replace a very costly (in terms of fuel weight) 
propulsive maneuver with an aerodynamic maneuver. These aeromaneuvering concepts enhance the 
utilization of space for scientific, commercial, national, and international purposes. 

Alas, however, it has come to be that hypervelocity flight is almost a dead language — to a con- 
siderable extent, the national institutional memory is failing. In order to guard some of what we have 
learned, a series of graduate level courses has been reintroduced. This monograph is an outgrowth of 
that effort. It contains material that I presented in a series of lectures in 1967 at Stanford University, and 
began to repeat in 1986. It was apparent to the students in 1986 that the material was reassembled on 
short notice, consisted largely of work that I was involved with over the years (as a matter of conve- 
nience — with no intent to slight the work of many others), and was attended by considerable difficulty 
for me and surely for the students. It was written in barely legible longhand on the chalkboard. More 
serious, over the years notebooks and technical papers had been misplaced because of carelessness; the 
burden of administrative dutyjhat d ealt in personpower, budgets, and partially realized plans; and 
because of many moves by me within NASA, including the worst case — a 2-year assignment to NASA 
Headquarters in Washington, DC. Thus, for example, my notebook on nonneutral gas effects was found 
by happenstance in the garage of a friend — a week before the lectures began. 

The lectures were televised in 1987 and 1988, which limited my use of the chalkboard and 
allowed me to be seated most of the time. Moreover, two previous students, Bertrand Delmas-Marsalet 
and Frederic Laurent Chalot, had prepared a succinct version of the class material by use of a word 
processor, which I tried to flesh out by use of a portable typewriter. Thus in every way, the lectures in 
1987 were improved — or so it seemed to me. 

The current text includes the previous material, but has been expanded to include some brief 
considerations of one-dimensional compressible flow without transport phenomena early in Chapter 2, 
flight in atmospheres of other planets, very severe flight environments (Chapter 9), explanatory material 
inserted where appropriate throughout the text, and additional material about hypervelocity flight at high 
altitude in Chapter 11. 


John T. Howe 
1989 
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. . .keep that which is committed to thy trust, avoiding the profane and vain 
babblings and oppositions of science falsely so called which some by 
professing have erred concerning the faith. . . 
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CHAPTER 1 


INTRODUCTION, CONTENT, AND APPROACH 


The purpose of this monograph is simple to state. We intend to develop the equations that 
describe the flow over an object in hypervelocity flight; assess the relative importance of all of the terms 
in the equations (and neglect those that are of lesser importance); describe the terms that treat the rele- 
vant, highly energetic, real-gas phenomena; solve the set of equations to obtain results, and examine 
those results to form some mental notions of trends and effects important to hypervelocity flight. We 
will examine the very severe flight environment of the Jupiter (or other planetary) atmosphere. Finally, 
some consideration will be given to flight at very high altitude where the flow field is relaxing chemi- 
cally and even thermodynamically at finite rates. It is important to state this because there is some dan- 
ger of forgetting what we are about in the process of developing the subject. 

Thus we will also find ourselves trying to map flight domains in which particular phenomena are 
important. To the extent possible, we will be forming notions about those processes that will be most 
sensitive and those that will be relatively insensitive — before the problem is solved. That is a good thing 
to do. 


Some of our flight experience in the Earth's atmosphere is depicted in figure 1 - 1 , which shows 
flight regimes in terms of altitude and velocity. On the lower left is a line that represents the flight of the 
supersonic Concorde. By contrast, as the monograph title suggests, we are primarily concerned with 
flight toward the right of the figure— hypervelocity flight. We say hypervelocity rather than hypersonic. 
The former includes both compressible and energetic real-gas effects, while the latter primarily denotes 



Figure 1-1.- Comparison of vehicle flight regimes in Earth's atmosphere. 
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compressibility effects — which can occur in a very cold gas that has little energy content per unit mass 
of gas (low enthalpy, if you will). For reference, low earth orbit circular satellite speed is approximately 
7 km/sec. The figure shows domains, trajectories, and some lines that approximately bound some 
physical and chemical flow-field phenomena associated with hypervelocity. Note four such lines. To 
the right of the lines labeled O 2 and N 2 , 10% or more of the oxygen and nitrogen molecules in the air 
behind a normal shock wave are dissociated to form atoms. To the right of the double vertical line at 
about 9 km/sec, electrons are stripped from the atoms and molecules to a significant extent, and ioniza- 
tion effects become important. Finally, at an altitude above about 50 km, the frequency of collisions 
among air particles is diminished to the extent that the rate at which physical and chemical processes 
proceed must be considered. That is, the air may be out of equilibrium chemically, thermodynamically, 
or both. This is the real world of hypervelocity flight. We will have occasion to refer to and modify this 
figure in more detail subsequently. 

Phenomenologically, in describing real-gas flows over objects in hypervelocity flight in an 
atmosphere, we must deal with these and other effects of the dissipation of the kinetic energy of the 
probe in the atmosphere — energy that is often sufficient to vaporize the probe. Thus the flow field may 
be complicated by the presence of gases that have vaporized from the probe surface, and in extreme 
cases by solid particles that have spalled from the surface because of thermal stress in the probe material. 
Even at that, the surface temperature is generally much lower than the gas temperature behind the shock 
wave. Thus a flow-field energy equation which treats real-gas thermal phenomena is essential. Often 
the energy equation is strongly coupled to the momentum equation because velocity terms are common 
to both, and the energy level is related to velocity squared. On the other hand, the solution of the 
momentum equation is not always greatly affected by the corresponding solution of the energy equation. 
Sometimes momentum equations can be solved without an energy equation to yield shock shapes and 
pressure distributions that bear some resemblance to real flight — depending on the domain and how par- 
ticular we are. 

Since the kinetic energy imparted to the atmosphere by a hypervelocity vehicle is large enough to 
break molecular bonds or displace electrons, we must consider chemically reacting flows (and also the 
possibility of electric fields and their effects) even at some distance removed from the displaced elec- 
tron. Moreover, in some regions of the flow, the gas temperature exceeds the apparent temperature of 
the Sun, and we suspect that the emission and transfer of thermal radiation should be considered in the 
gas — which may be complicated by chemical reactions that give rise to local sources and sinks of 
energy. 


Nevertheless, it is very important that we not consider the analysis of a hypervelocity flow field 
to be intractable; it actually involves a number of scientific disciplines that are tractable — fluid dynam- 
ics; chemical kinetics; electrodynamics; transport of mass, momentum, energy, chemical species, and 
surplus charge; and radiative transfer. Further, we must consider local thermodynamic nonequilib- 
rium — do we need more than one local temperature to characterize the gas mixture? It appears that it is 
necessary under some circumstances. 

The task of developing each discipline or each line of thought in detail in one small monograph 
is out of the question. Thus we are compelled to introduce each discipline at the subjective point where 
it suits our purpose, combine these interrelated physical and chemical phenomena in a tractable way, and 
terminate it at a reasonable point. We will leave further development of this continually emerging sub- 
ject to those with greater perception — to the students during their careers. 



Specifically, we adopt the following approach. Rather than begin with the Boltzmann equation, 
we begin with the hydrodynamic equations, or "the equations of change" which can be derived from the 
Boltzmann equation (refs. 1, 2). Nor do we take moments of that equation to obtain transport expres- 
sions, but use the Chapman-Enskog result (refs. 1-3)— or something akin to it. Radiative transfer theory 
is not derived from quantum notions, but rather from developments from stellar literature (which sur- 
prisingly can be shown to have been developed independently by the paint and paper industry). Almost 
exclusively, we adopt the macroscopic approach obtained from the results of kinetic theory— and say 
very little about that theory. When at last we begin to invade the inner privacy of atoms and molecules 
explicitly, we will say a few words and draw a few lines and withdraw with our sensibilities partially 
intact. That invasion is a subject to be dealt with separately. All that we need to know is where in the 
flight domain it becomes important, and understand some features of that domain in an approximate 
way. Thus if the reader or student has a nodding acquaintance with compressible fluid mechanics, and a 
reasonable tolerance for pain, he or she is qualified to begin. 


REFERENCES 


Hir^hfelder, J. O.; Curtiss, C. F.; and Bird, R. B.: Molecular Theory of Gases and Liquids. Wiley, 

2. Bird, R. B.; Stewart, W. E.; and Lightfoot, E. N.: Transport Phenomena. Wiley, 1960. 

3. Chapman, S.; and Cowling, T. G.: Mathematical Theory of Non-Uniform Gases. Seconded., 

Cambridge Univ. Press, 1951. 
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where Ci = pi/p is the local mass fraction of species i in the mixture. An equation of state relates the 
thermodynamic properties 


I 


i=l 


i=l 4 


where R is the universal gas constant and Mj is the molecular weight of species i. 


(7) 


A shock wave normal to the one -dimensional flow in figure 2-1 can be described for present 
purposes by use of the density ratio across the shock, which we will call e. 


Poo 

Ps 


e 


( 8 ) 


where the subscript °° refers to conditions ahead of the shock, and s refers to conditions behind the 
shock. Note that e is affected by a change of chemical composition across the shock, which will be 
treated subsequently. Also, it is notable that E is usually much less than unity. Thus the compressible 
flow relations (1), (3), and (4) can be used to form the relationships 


u s = eu^ 

(9) 

Ps Poo Poo^oo( 1 £) 

(10) 

2 2 
h h U - U s 

hs-hoo- — -y 

(11) 


These relations will be used in Chapter 10 where flows that are not in chemical equilibrium will be 
examined on a one-dimensional basis. 

Flow of a compressible gas over an object in flight is usually not one-dimensional, especially if 
the object tends to be blunt. Moreover, the shock wave over an object in supersonic flight will generally 
tend to be a curved shock which is not everywhere normal to the flow nor parallel to the surface of the 
flight object or body. In Chapters 3 and 4, it will be helpful to have quantitative relationships for all the 
thermodynamic and flow-field variables across such a shock wave for compressible flow, as well as 
some "strong shock" approximate expressions for hypervelocity flow. The reason for the latter is that 
some of our order-of-magnitude arguments for sizing terms in the flow-field equations can be facilitated 
by these expressions. 

For either a two-dimensional or axisymmetric object flying at hypervelocity, U, the bow shock 
wave in the stagnation region of the body is shown schematically in figure 2-2. The stagnation region 
could be the leading edge of a wing or a fin, or it may be the nose of a fuselage. It could be the inlet to 
an engine as well. Significantly, it is a region of a nearly normal shock wave, which implies strong 
compression of the atmospheric gas and potentially serious thermal effects in the flow field. Since the 
shock wave is generally not parallel to the body, the geometry of each, and the components of velocity 
normal and tangential to both the shock and the body are shown in figures 2-2 and 2-3. Note that 


6 


CHAPTER 2 

SOME USEFUL COMPRESSIBLE FLOW RELATIONS FOR AN IDEAL GAS 


^ wl i! L be helpful subsequently to have some simple relationships concerning compressible flow 
to refer to. The steady one-dimensional flow of an ideal compressible gas in a constant-area duct 
wherein there is no transport of mass, momentum, energy, or chemical species should be familiar to the 
reader (ref. 1). That flow is depicted in figure 2-1, and is described by equations of the continuity of 
mass (where p is gas density and u is the gas velocity). 


PoO‘ “<*, Poo' h OO. "Too 

UNDISTURBED ATMOSPHERE 
COMPOSITION N 2 


P v u s . P s . T, 

SHOCKED GAS 
n 2 . n. n 2 , N + , e‘ 


SHOCK WAVE 


Figure 2-1.- Thermodynamic and chemical change across a normal shock. 


pu = constant 


a momentum balance (where p is the pressure and x is the one-dimensional coordinate) 


(its integrated form is) 


pu 


dn = _dp 

dx dx 


( 1 ) 

( 2 ) 


pu 2 + p = constant ( 3 ) 

and an energy equation (where h is the static enthalpy, and H is the total enthalpy) 

2 

h + ^- = H = constant ( 4 ) 

The static enthalpy for a component species of a perfect gas is related to the temperature, T, by 

c Pi dT + hj 0 ( 5 ) 

where Cpi is the specific heat at constant pressure, and hj is the heat of formation of that component 
of the gas at a reference temperature (zero, for example). The enthalpy of the gas mixture of v 
components or species is simply 
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Figure 2-2- Bow shock wave — not everywhere parallel to the body. 



Figure 2-3 - Velocity components behind shock; shock not parallel to body locally. 

shock-related velocity components are here denoted by the symbol q for clarity. Moreover, the local 
radius of curvature of the body is denoted by R here. 

Figure 2-3 resolves the velocity components in shock coordinates behind the shock into those in 
body coordinates. The small angle between the directions normal to the body 0b and the shock 0 S is 
denoted by co. It can be shown from the ordinary Rankine-Hugoniot compressibility relations for 
curved shock waves (similar to those shown previously for normal shock waves). From considerations 
of the conservation of mass 


Pooqn*, - PsQn, 


(12) 



Poo — 

qns= p7 qn ~ = eUcos 0S 

* ^ 


qt = U sin 0 S 

( 14 ; 

(q t unchanged across the shock) 

u s = qt cos 0) + q ns sin co = U(sin 0 S cos co + e cos 0 S sin co) 

( 15 ) 

-v s = qn s cos 0 ) 

- q t sin co = U(e cos 0 S cos co - sin 0 S sin co) 

(16) 

Also, from a momentum balance, 


P« + p®q noo - Ps + Ps^ng 

(17) 


which is similar to the development for equation (10) in the direction normal to the shock. The strong 
shock approximation 


P® « Ps 

used in equations (12) and (13) leads to 

Ps = P®U 2 (1 - e)cos 2 0 S 

The energy relations across the shock can be expressed as 

ho ° + 2 ( q n~ + q ?) = hs + + q?) 

or 


( 18 ) 


( 19 ) 


( 20 ) 



Use of equation (13) in (21) yields 


h s - = 1 q2_a _ e 2) 

Again, for a strong shock, hoo can be neglected compared with h s , and 

hs = j U 2 cos 0 S (1 - e 2 ) « U 2 cos 2 0 s 


( 21 ) 


( 22 ) 


( 23 ) 


An equation of state, needed to calculate e, will be considered subsequently. 

Some simple geometric relationships will be useful later. Figure 2-4 shows a portion of the body 
surface (y = 0), and a line, y = constant. 
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Figure 2-4 - Segment of body surface. 
On the body (y = 0), from the right side of the figure, 

(dr) y= o = sin ct dx 


(r) y =o = 



sin CT dx 


For a surface where y = constant (left side of figure), 

r(x,y) = r y _o + y cos CT 

r 

r(x,y) =1 sin G dx + y cos CT 

dx = -R dCT 

£& = _!_ 
dx R 

These relationships will suffice for subsequent reference purposes. 
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CHAPTER 3 


CONTINUITY AND MASS FLUX CONSIDERATIONS 


This chapter begins the formulation of the differential equations that describe real-gas flow fields 
over objects in hypervelocity flight. Simplifications will be made by use of relationships from Chap- 
ter 2. Although the simplification is not great, it is indicative of how it is done for all of the flow-field 
equations. The development will not be repeated for the other equations. 

Continuity relations will be written for any of the individual species present which will include a 
term that represents the formation or depletion rate of that specie by chemical processes. It will be 
summed to yield the familiar global continuity equation for compressible fluid flow. It will be special- 
ized to describe only chemical elements, which will add a term that represents the transport of each ele- 
ment relative to the mass averaged motion of the mixture. Further specialization can be made for each 
chemical specie. Charged species are considered; we will obtain the surplus charge diffusion equation. 
From that we will obtain the electric current and electric field. The latter gives rise to forces that affect 
overall motion and individual transport. We will do this in some detail because it is simplest to do so at 
this point rather than later, when the equations and concepts are more complicated. There is some dan- 
ger that this development may appear to be tedious. However, perhaps the redeeming feature of this 
chapter is the wealth of useful relationships that can be derived from continuity considerations applied to 
real-gas flows. 


GENERAL CONTINUITY EQUATION 


An element of volume in the coordinate system x,y,z is shown in figure 3 - 1 . The coordinates 
are orthogonal at their intersections, but are not necessarily Cartesian. The coordinate units themselves 
do not have to be of physical length, but the product of the coordinate unit and its respective metric 
(hl.h2.h3) is an arclength in that coordinate direction, and the square of the elemental arclength from the 
point x,y,z to x + dx, y + dy, and z + dz is expressed in terms of hi, h2, and h3. 

hi=hi(x,y,z) etc. (1) 

arclength (ds) 2 = (hi dx) 2 + (h2 dy) 2 + (h3 dz) 2 
For convenience, consider the chemical specie i. 

uj = absolute velocity of species i 


uj = u + Uj (2) 

average relative 

The components of the absolute mass flux of species i in the coordinate directions are 

PiUj = pju + Pi Uj (3) 
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Figure 3-1.- Volume element. 


or 


similarly 


pjUj = piU + Jix 

(4) 

PjVj = PjV + J iy 

(5) 

Pi W i = PiW + Jiz 

(6) 


Summing equation (4) over all species i = 1 to v yields the mass average mass flux in the x direction. 

v 

pu = I PW (7) 

i=l 


Thus, 


X P* u i = 4 j Pi + X J ix = PU +X J ix 


i=l 


i=l 


i=l 


i=l 


( 8 ) 


0 


where the Jj are components of the mass flux vector relative to the mass average in each direction. 
Note that each component of the species mass flux vector relative to the global mass flux, summed over 
all species, is zero. 

The mass flux of i normal to a surface is 


mi = cjm + Jj 


iy 


(9) 


(as before Ci = mass fraction = pi/p.) 


u i r 
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A mass balance of species i in figure 3-1 is 


0 0 -N 

— (pih i dx h2 dy h3 dz) = - — (piUih 2 dy I13 dz)dx (pjVihi dx h3 dz)dy 


dx 


- ^ (PiWihi dx I12 dy)dz + Kjhi dx h2 dy 113 dz 


( 10 ) 


where the rate of change of the mass of species i per unit volume is shown as the difference of inflow 
and outflow, plus the rate of mass production per unit volume, that is 

Kj = mass rate of production of species i per unit volume. 

Thus 


dPi- + --l 

dt hih2h3 


0 0 0 
— (piUjh2h3) + — (pjVil?ih3) + — (pjWjhih2) 


-K i= 0 


(11) 


which is the conservation equation for species i. This form of the species conservation is basic, from 
which other relations will be derived. First let us derive the global conservation of mass flux. 


Global Continuity Equation for the Steady State 

Assume d/dt = 0 (steady state) and d/dz = 0 ( 2 D or of revolution); for an axisymmetric body, z 
represents the cross-flow direction, and the cross flow is taken to be zero. Then substitute equations ( 4 ) 
and ( 5 ) into equation (1 1): 


d d 

— [h2h3(piu + J ix )] + — [hih3(pjv + Jj y )] = hih 2 h3Ki 

0 0 r) r) 

— (h 2 h 3 piu) + — (hih 3 piv) = - (h 2 h 3 Jix) (h^fy) + hih 2 h 3 Ki 


( 12 ) 

( 13 ) 


The summation over all the species i gives 


d_ 

dx 


h 2 h 3 u X Pi + 
i=i / d y 




V \ 


hlh3V^ Pij = -^|h2h3S J ixj-^|hlh3E J iyj + h lh2h 3 X K i (14) 


Equation (8) shows that 


v v 

X J ix = X J iy =° 
i=l i=l 


( 15 ) 
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as was mentioned previously. Since 


^ pi = p => — (h2h3pu) + — (hih3pv) = hih2h3 j Kj 


i=i 


dy 


(16) 


i=l 


Consider the following chemical reaction: 

A 2 — » 2A 

The rate of mass loss of A 2 equals rate of mass formation of A. We have 


i=l 


(17) 


(18) 


in general. Then equation (16) becomes 


9 9 

— (h 2 h 3 pu) + — (hih 3 pv) = 0 
dx dy 


(19) 


hi=K=l+|- 


The metrics which convert coordinate changes to lengths are for the x, y, and z direction respectively 

( 20 ) 

h 2 =l (21) 

h 3 = ri (22) 

where j = 0 (two-dimensional) or j = 1 (axisymmetric). (The last assumes no cross flow.) Thus 

(23) 


J- (purJ) + (Hpvri) = 0 
dx dy 


which is the steadv-state global continuity equation . 


Further Simplifications 

Equation (23) can be simplified to the "thin shock layer" approximation by an order-of- 
magnitude argument as follows: 


dr* : d(pu) ar : d(pv) 
pu — + rJ - + HrJ -^-2- + ptf 
dx dx dy 




\ dy dy 

From equations (27) and (29) of Chapter 2 and equation (20), above: 


(24) 


: ! 1 I I : 


14 



1 + = rt sin a 

R 


( 25 ) 


dr* 

dx 


da 

= sin a - y sin a = sin a 
dx 


For order-of-magnitude estimates, we proceed as follows. Assume that shock and body are 
essentially parallel. From 0 S = 0b (to = 0) and equation (15) of Chapter 2 where we assume e varies 
slowly with respect to x, 


^ (PsUs) = U cos 0 S 



(from fig. 2-2 of Chapter 2) and 


R d0 s = R d0b * dx 



PooU cos 0 S 
eR 


From equation (16) of Chapter 2: 



O 


Ps v s ~ Pw v w 

5 


o 


Po„U cos 0 S 

5 


(26) 

(27) 


(28) 


(29) 


where the + corresponds to p w v w — PsV s ~ +p«>U cos 0 S (that is, surface mass addition as high as the 
free- stream mass flux is allowed), while the - sign corresponds to p w v w = 0 (no mass addition is 
allowed as well), and 8 is the shock standoff distance. 

From equation (26) of Chapter 2: 


— = cos a (30) 

dy 


Say, 8 ~ eR, with e = Poo/Ps- 

Substitute into equation (24) using equation (25), behind the shock and employ the right side of 
equation (14) of Chapter 2 for u s : 


0 = 


PooU 


sin 0 S K sin a + ri cos 0 S ± HrJ cos 0 S - p^U cos dirt cos a + ^-1 
£K £K \ R/ 


PooU 


(31) 


Let 


H = 0(1) 

(32) 

O 

II 

(33) 
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At the most 


sin 0 S ~ 0(1) 
cos 0 S ~ 0(1) 
sin a and cos a ~ 0(1) 

Multiply equation (31) by e/pooU; at the most, the terms are of the order 


(34) 

(35) 

(36) 


(l) + (l)±(l)-e(l + 1) = 0 

Since e « 1, the last term is negligible and so is the last term in equation (24), which it represents. 
Thus r can be regarded as rb(x), and equation (23) can be written as 


|^(pui J b ) + «r J b |-(pv) = 0 ( 38) 

M ^ thg Steady-State global continuity equation for two- d imensional or axisvmmetric flow 
reduced by order-of-magnitude arguments for the thin layer approximation, becomes, from 
equation (23), 




(39) 


where 

n> = ibOO (40) 

Now reduce the species continuity equation for the steady state in a similar way. Rewrite equa- 
tion (13) using equations (21) and (22), and the definition 


Cj = pj/p ( 41 ) 

— (ricjpu) + — (rJHcipv) = -~- (rjj*) - -t- (HrJj iy ) + HdKi (42) 

Ci ~ (pmJ) + purJ . + Hpi vri ^ + q (Hpv rj) = - (rij ix ) - 1- (HrJ J iy ) + HrJKi (43) 

The underlined terms sum to 0 because of equation (39). 

We then divide by rJ, consider rJ ~ rJ(x), and (H = 1). We neglect the mass-flux derivatives in 
x-direction, and obtain the steady-state species continuity equation for species i for two-dimensional or 
axisymmetric flow 


pu^+pv|^ = -A(j iy) + Ki 


dx 


(44) 


n '■ r 
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This is the counterpart of equation (11), for steady-state, two-dimensional, or axisymmetric flow for the 
"thin layer," where as before 

Jjy = mass flux vector in v direction 

K; = mass rate of production of species i per unit volume 


Elemental Continuity 

We next consider the element p. 

Define : 

M p = atomic mass of element p 

Mj = molecular mass of species i 

Ppj = number of atoms of element p in molecule i 

Also define a set of numbers apj as 


(45) 


where a pi is the mass of element p in species i. Multiply equation (44) by a pi and sum all species: 

(46) 


V dcj vp 0Cj ^ dJiy Y' V 

P" 2- “Pi 57 + P v Z «Pi3- = - E “Pi 3^ + 2, “Pi K ‘ 


i-i 3* r ft F ‘3y ft * 3y ft 

Further define the elemental mass fraction 


and the elemental relative mass flux 


Note that no elements are created 


c'p = X a Pi c i 
i=i 


J' p - X a Pi J ‘y 
i=i 


(47) 


(48) 


° = X a Pi K « 

i=l 


(49) 
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( 50 ) 


The last expression is true for ionizing flow, as will be demonstrated. Thus we have 


pu 





which is the simplified elemental diffusion equation for the thin shock layer. 

Thus we have defined two sets of constants, a pi and (3^, which are the mass of element p in 
species i, and the number of atoms of element p in species 1 , respectively. We have also defined the 
elemental mass fraction and elemental mass flux as c' p and J' p , respectively. Moreover, we have 
asserted the relationship in equation (49) which states that the mass of element p is unchanged by the 
production of species i, which will be demonstrated subsequently. The steady-state elemental diffusion 
equation for two-dimensional or axisymmetric flow reduced by order-of-magnitude considerations is the 
result, as is shown by equation (50). 


It is significant that the production term which appeared in the species continuity equation (44) 
does not appear in the elemental continuity (or diffusion) equation (50). There are two advantages that 
the elemental diffusion has: (1) the specie production rate need not be specified (which is of practical 
utility when the flow is in chemical equilibrium) and (2) there are usually fewer elemental diffusion 
equations than there are species diffusion equations. 

As an exercise, to demonstrate equation (49), consider the following reaction: 


0 — » 0 + + e - 

(51) 

where element p is oxygen, and the species are 0, 0 + , and e~. Then 


0t Oq = 1 

(52) 

Mn 

0o+ = M 
iV1 o + 

(53) 

«'Oe- = 0 

(54) 

V 

X«PiK, K 0 + Mot V = M ot (K ° M O* tK O* Mo) 

(55) 

From equation (18), 


V 

£ Ki = 0 = K o + K +K_ = 0 

i=l e 

(56) 


Electrons and ions are produced in the mass rate proportions 


By use of equation (57) in the last equality of equation (56), 

M K °* + ^ K 0 *H 0 { 1+ iv) 

Use equation (58) in the last equality of equation (55) 


i-pft-M 1 


M \ 

" M o + V 

1 + m‘ +M °V 

i=l 0 + 


0*/ J 


Factor Kq+ 

rr KM 0* +M-> +Mo] = 0 = £ «P,K, 

0 + i=l 

which demonstrates equation (49) even for ionizing processes. 


Charged-Particle Continuity 

Finally, let us consider electrically charged particles. 

Again, define a set of numbers: 

C 0 i = charge number of a particle = 0 neutral 

-1 electrons 

+1 ions 

+2 doubly ionized 

Multiply equation (44) by (-ojj) • (M e _/Mi) and sum over all species i: 

a 

pU 3T 

Define : 


"v M _ 


"v M _ 


"v M _ 

Z (“^i) TT" c i 

i=l Mi 

o 

+ pv — 

dy 

S Mi 

_ <7 

3y 

2 (_<0i) w Ji y 

i=i 111 


M 


+ S(-«i)Tf K 

" Mi 


i=l 


V 

i=l 


r.-X 

i=l 


I M 

i* 

COi 


Ci 


1 Mi w 
/ M 


©i 


1 mT ,y 


(58) 

(59) 

(60) 


i (61) 

(62) 

(63) 
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Rewrite equation (61) 


pu 


9c' e 
Bx 


+ pv 


9c' e 
By 



which is the surplus charge diffusion equation. 


We should take note of five things: 


1. We have set 


v 


Z 


M 

(-coi)-f Kj = 0 
Mj 


which can be demonstrated. 


2. g'p, is proportional to the surplus charge density, o^ . as follows : 

v v v 

Pc = e X ®i n i = en ^ coj — = en ©ixj 
i=l i=l n i=l 


(64) 


(65) 


(66) 


where e is the electron charge, tt>j is the charge number for species i, nj is the number of species i per 
unit volume of fluid, and xj = nj/n. From the relations of equations (62) and (66) and Avagadro's rule, 




m; = ^ , n^i = pj , DM = p 


, v* , s M a _ M A 
C e = £ (-0)i) jfr <H = -jf £ (-tOiXi) = - 1 


i=l 


M 


i=i 




e PL 


! Pc 


where L is Avagadro's number. For a neutral gas, even though it is ionized, 

c' e = 0 


(67) 

( 68 ) 


(69) 


or there is no surplus charge locally. 

3- is proportional to the v co mponent of electrical current J y . From equation (5) we have 

Ji y = p i Vi = n i m i Vi = ^Ml- Vi (70) 

From equation (63) we have 


v M niVi 

J e — X (~®i) 
i=l 


eL 


M v 
e~ V 


eL 


i=l 


M _ 

ecoinj Vj = J v 

eL J 


(71) 
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which also defines the y component of the current as 


J y = X ecOjnjVi (72) 

i=l 

4. By summing over all species, it can be shown that 

c’e + S c 'p =1 (73) 

P=i 

5. By summing the elemental diffusion equation (50), adding the surplus charge diffusion equa- 
tion (64), making use of equation (73), and integrating leads to the useful relation 

n 

J’ey + Z J 'py = 0 ( 74 > 

p=i 

Thus we have derived several relations which will prove to be useful from continuity considera- 
tions for the conservation of mass, species, elements, and electric charge. Subsequently, we will com- 
bine these with the equations of momentum and energy. 
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CHAPTER 4 


THE EQUATIONS OF CHANGE; THE HYDRODYNAMIC EQUATIONS 


In the previous chapter we derived some equations which properly belong in combination with 
those of this chapter. That is, in chapter 3 we have derived the chemical species continuity equation, 
equation (1 1), in generalized orthogonal curvilinear coordinates; specialized that to two orthogonal 
coordinates for the steady state, equation (13); summed that over Jill species to yield the corresponding 
global continuity equation, equation (19); specialized it to two-dimensional or axisymmetric flow, 
equation (23); and finally simplified that by an order of magnitude analysis which allows mass addition 
at a rate comparable to the free-stream mass flux. From this reduced form of the global continuity equa- 
tion, equation (39), we rewrote the corresponding reduced species continuity equation, equation (44); 
summed to yield the elemental continuity equation, equation (50), and derived the surplus charge diffu- 
sion equation, equation (64). Moreover, we established some useful statements expressed by equa- 
tions (18), (49), (65), (73), and (74); related charge density to surplus charge, equation (68); and 
obtained expressions for the electric current, equations (71) and (72). All of this was obtained from 
continuity considerations, allowing for mass addition at a boundary, transport phenomena, and species 
production at a finite rate. 

In this chapter, we will combine statements concerned with momentum and energy considerations 
with the above. But we will not derive these in detail as we did in chapter 3. Rather, let us simply state 
the procedure for obtaining these expressions, and write the result (refs. 1-3). 


Governing Differential Equations 

When the equations of change (ref. 4, p. 698) are written in orthogonal curvilinear coordinates and 
applied to the hot gas between a body in hypervelocity flight and its bow shock wave (fig. 4-1), the 



Figure 4-1.- Typical bluff configuration. 
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form^m leads to two momentum equations and an energy equation, each of which contains about 
00 terms. To simplify, we have eliminated all of the higher order terms by an order of magnitude 
analysis (ref. 1). The body configuration used for this analysis typifies a bluff flight object, or the nose 
or leading edge of a vehicle in which aerodynamic heating may be severe. The individual terms in the 
equations were sized in six regions of the flow: both behind the shock wave and near the surface; in the 
stagnation region, near the curved shoulder; and on the conical flank. Basic assumptions are that the 
shock layer is thin and that mass addition at the body surface can be as large as the ffee-stream mass 
x T UX 'u emi j order e, (eRe)* 1 , and (eRe 1 )" 1 are neglected compared with unity (Re is Reynolds 
Number and Re’ is based on bulk viscosity). Details of the simplification are too voluminous to dis- 
play. However, the resulting hydrodynamic equations (which can be compared to simple hydrodynamic 
equations (ref. 5, p. 319)) are 

x-momentum 


_ du du d 

pu dT pv aT^ 


*d y dx 


( 1 ) 


where p is viscosity. 
v-momentum 


dp _ pu 2 
dy R 


XniYj 

i=l 


(2) 


where Yi is the y component of a body force (such as an electrical force which arises from an induced 
electric field caused by charge separation). 


energ y 


pu ^ + Hpv ^L = -hA- (q ) + JL 
dx ^ dy dy Hy dy 


pu 


du 

dy 


divq r + y' Y,^ 

tl m > 


(3) 


conduction radiative 

transfer 

where the underlined terms concern transport of energy by conduction, and energy transfer by gaseous 
and surface radiation. An alternate form of the energy equation will be presented subsequently. To 
these we add the global and species continuity equations, equations (39) and (44) from chapter 3: 

continuity , the elemental version 


^(punJ + ^(pvnO-o 


(4) 


or the species version 
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( 5 ) 


5cj 3ci 3 

pu ar +pv a7 = -a7 <J ^ ) + Ki 


These are essentially the boundary-layer equations for a reacting gas, but importantly, they apply 
throughout the flow about the body and are valid for strong ablation. Because ionized species are pres- 
ent at high-speed flight, we have retained force terms arising from electric fields in the y-momentum 
and energy equations. We write the elemental continuity equation for element p and for electrons 


pu ^ +Kpv 5 = -|r (MJ p > 


(6) 


which is the elemental form, and 


9c' 


9c' 


pu — +pv — = 
9x 9y 



(7) 


which is in the surplus charge continuity form. Equation (6) is useful for the case of chemical equilib- 
rium, for which it would replace equation (5); the advantage being that there are fewer elements than 
there are species. Thus species concentrations would be calculated by other means. Equation (7) would 
be used if there is significant electrical charge separation. 


To equations (1)— (5), we add the equation of state (for the sixth unknown): 


„ piRT Ci ^ nitiiiRT RT , ^ 

Pi= -M- = s pRT= ir = " i r = ' ilT 


( 8 ) 


where R is the universal gas constant, L is Avagadros number, nj is the number of species i per unit 
volume, and k is the Boltzmann constant. Dalton's law relates partial and static pressure 


p-Zb-p-tZ^ 1 ^ 

i=l i=l m * M 


(9) 


Other useful relations are 


H 



h = X c ibi 

i=l 




Cp dT + Hj 


( 10 ) 

( 11 ) 


( 12 ) 


and for thermochemical equilibrium 
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(13) 


®i — Cj(p,T, Cp) 

In these expressions, H and h are total and static enthalpy, R is the universal gas constant and M is 
the mixture molecular weight (defined by eq. (9)). Equation (13) is for chemical equilibrium, and can be 
expressed by the law of mass action or by free energy minimization relations. For chemical nonequilib- 
num, reactions proceed at rates that can be expressed in the Arhennius form: first the reaction rate 

k = AT«e-^ (14) 

secondly 


K e = BTP e" 4 ^ 


(15) 


The latter is the equilibrium coefficient or law of mass action. 


Energy Equation Expressed in Temperature Form 


It is often convenient to write the energy equation, equation (3), in terms of temperature as follows 
From equation (1 1): 


0 h = y 

m 


dhi 8cj| v< 

Ci 3“ +hi 3^ = X 
3y 3y 


3§T 

'3T3y "‘dy 


3T y 3hi -r. 3cj 9T ^ Y* 3ci 


or 


dh 

3y 


- 3T v , 3cj 
c p^- + 2, h i3- L 

dy S 3y 


where c p is called the frozen specific heat, and is 


(16) 


(17) 


v 

X CjCp 
i=l 

Differentiating equation (10), where the partial derivatives of v are negligible compared to derivatives 
of u and H leads to 


9H ah , a u 

aTaT^ ^ 


vi : it ; 
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( 19 ) 


3H _ 3h 3u 
3y dy U dy 


Substituting equations (17)— (19) in equation (3) with H = 1 yields 


pd 


_ 3T I v . ^ c i 3u 

'^l + ^ hi a 7 +u sr 


3x 


i=l 


+ p'i 


_ 3T v> , 3c; 3u 

Cp a7 + £ hi 3h + u 37 


= — — q v + — |}lu— — divq r + 5' Yj — 

3y 4y 8yr 5y| 4 ~ mi 


We rewrite equation (20): 

7 \ 

pul 


_ 3T 

Cp 


dx 

LV / 


X h i 

i=l 


dCi 

dx 


H — 3T v -1 , 3cj 


2 3u 3u 
PU 3^ puV 3y 


d 3 

'dy’* 57 


1 3u' 

“ u 57 

\ / 


-divqr + X Y i^f 

i-1 mi 


Multiplying (1) by u, we obtain 


2 8u , 3u d 3u . 

PU ^ + PUV ^7 =U ^7 P ^7 | ■ U, 


u \ 

ap 


3x 
v / 


Subtracting equation (22) from equation (21) gives 


pd 


/ \ V 


_ 3T 
Cp 3 — 
v 3x / 


X h i 

i=l 


dx 


j - V 1 l ^ c i 

pv|cp a7 + ^ hi a7 


- 57 ^ 
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i=l 


m; 


n 


which is the temperature form of the energy equation . 


(20) 


(21) 


( 22 ) 


(23) 


This completes the formulation of the differential equations reduced by order-of-magnitude argu- 
ments. They apply throughout the flow about the forebody shown in figure 4-1, and are representative 
of the flow wherein severe heating may occur, such as forebodies in axisymmetric or two-dimensional 
flow. This would include forward stagnation regions, leading edges of lifting surfaces (and perhaps 
engine inlets), and bluff inclined surfaces. Among other things, the equations allow for chemical reac- 
tions (equilibrium or nonequilibrium), electrical charge separation (and the electrical fields that are 
induced), and ablation at rates up to the free-stream mass flux rate. Subsequently, various terms in the 
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equations pertaining to heat conduction, radiative transfer, electric Helds, and the transport of mass by 
four kinds of diffusion (concentration gradients of species, thermal, pressure, and (electrically) forced 
diffusion) will be described and assessed as to the relative importance of each. This will be done in the 
next chapter. 

Thus we have the governing differential equations concerning mass conservation (eq. (4)), state- 
ments of the x and y components of momentum (eqs. (1) and (2)), the energy equation expressed in 
terms of total enthalpy (eq. (3)), the species continuity equation (eq. (5)), and the equation of state (eq. 
(8) or (9)); effectively six equations in six unknowns, p, p, u, v, H, and Cj (at least for an electrically 
neutral gas). Alternatively, equation (6) can replace equation (5); the energy equation, equation (3), can 
be expressed in terms of temperature and can be replaced by equation (23); and for a nonneutral gas, the 
surplus charge equation, equation (7), would be used along with an electrodynamic or electrostatic 
equation which will be presented later. 


Illustrative Boundary Conditions 

For the balance of this chapter, we consider some illustrative boundary conditions for these differ- 
ential equations. These are shown in equations (24)-(34) for y = 8 (behind the shock wave), and y = 0 
(at the wall). Recall that e is required to calculate all the unknowns listed previously behind the shock 
wave. That requires an equation of state (eq. (8) or (9)), which in turn requires some knowledge or 
assumptions about the chemical state of the gas. For example if we assume the gas to be frozen across 
the shock, the species are unchanged, equation (29); while for equilibrium we would have equation (30). 
Similarly, temperature and electron concentration are treated according to whatever assumption we make 
(eqs. (31), (32), and (34)). It can be seen that these boundary conditions can be written in various forms 
depending on the assumptions made and the choice of alternate dependent variables, as long as we are 
consistent. 


Boundary conditions at the wall are also shown. Equation (36) illustrates steady-state mass addi- 
tion at the wall in response to flow-field heating; that is, the excess heating is compensated by the 
appropriate mass injection or ablation at the wall. Figure 4-2 is a schematic diagram corresponding to 
equation (36). 

Boundary conditions - At y = 8 


u = u s 

(24) 

V = v s 

(25) 

P= Ps 

(26) 


H 



(27) 


Cj — c is 


(28) 
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for chemical nonequilibrium 


Cj_ — Cj s 

(29) 

while for chemical equilibrium 


Cis =/ (p.T.Cp) 

(30) 

For chemical equilibrium with species diffusion 


T = Ts(h s ,p Sl Ci s ) 

(31) 

or with elemental diffusion 


T = Ts(h,p,c^c^ 

(32) 

c = c - c 
P ps 

(33) 

c = c' = c 

e es e«, 

(34) 

At y = 0, with no velocity slip at the wall 


u = 0 (no-slip condition) 

(35) 


The steady state ablation may be written 

. . -qyw-qradw-YwOTl 

m w = (pv) w = — —r r (steady state) (36) 

— ^interior 

where y w is the wall emissivity, and h w includes the enthalpy of all gaseous species at the wall. Thus 
the total enthalpy at the wall is 


V 

H w — h w = (Cjhj) w (37) 

i=i 

where 


r 

hiw — I Cp dT + 


h 


(38) 


From equation (5) of chapter 3 written for atmospheric species which we say do not penetrate the wall 

miw = cjm w + Jiyw = 0 (39) 

The equivalent statement for atmospheric elements would be 
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rn p = CpTTi w + J pw = 0 (40) 

th^waS” (4 °^ W ° Uld n0t CqUal Zer ° if element p ori g in ated in the wall. For no flux of electric charge to 


me c ( xn w + Jg w — 0 

Steady State : (proof of eq. (36)) 

~ q yw T' w oT * (pv) w h w -q radw 



Figure 4-2 - Ablating material in coordinates relative to the surface (wall). 
For steady state ablation, the energy out = energy in 

(pv) w h w + Y w oTi = -(q y + q r Xv + (pvh) interi or 
The mass balance is expressed by 


(41) 


(42) 


(p v )w — (pv)interior (43) 

Thus, we obtain the relation of equation (36) given above: 


(pv) w = 


(q y + q r Xv ~ YwCtT^, 


h w -h 


interior 


(44) 


It may be noted that the boundary condition (eq. (39)) states that chemical species which originate 
in the atmosphere do not penetrate the wall. Equation (40) is not zero for elements which originate in 
the wall, and is zero for elements that do not originate in the wall. Of course for transpiration suction in 
which the boundary layer is being sucked into the wall, the zero does not pertain. Note also that equa- 
tion (40) does not require the wall to be neutral electrically. The boundary condition (eq. (41)) states 
that there is no absolute flux of charge into the wall. These statements are illustrative only and can be 
rewritten to describe any particular requirement. 


In summary, the governing differential equations are typically equations (1>— (5) (with alternatives 
eqS " /o\ 3nd ^ re P lac ing eq. (5), and eq. (23) replacing eq. (3)) and the equation of state; either equa- 
tion (8) or (9). Illustrative boundary conditions behind the shock wave are equations (24)-(28), with 
alternate choices equations (3 1)— (34). At the wall they are equations (35)-(39), with alternate choices 
equations (40) and (41). 
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CHAPTER 5 


TRANSPORT PROCESSES AND EXPRESSIONS 


At this point, the differential equations concerning conservation of matter, momentum balances in 
two coordinate directions, and the exchange of energy in the flow field have been developed in their 
most simple forms. The equations contain terms that represent viscous effects (|i, in eqs. (1) and (3) of 
chapter 4), the transport of mass by diffusive processes (eqs. (5), (6), and (7) of chapter 4), forces that 
are induced electrically (eqs. (2), (3), (5), (6), and (7) of chapter 4), and the transport of energy (eq. (3) 
of chapter 4). In this chapter, we will describe those terms and assess their relative importance. Other 
terms also appear that describe the transfer of energy by the emission and absorption of radiation, and 
the formation and disappearance of species by chemical reactions. These will be treated in separate 
chapters (6 and 10). 

With respect to transport, we will not derive the transport expressions in this monograph. They are 
a subject unto themselves. An outline of the derivation appears in references 1, 2, 3, and 4. The expres- 
sions for the transport terms were derived independently by Chapman and Enskog. They were able to 
derive the Navier-Stokes equations from the Boltzmann equation, and in so doing arrived at formulas for 
the transport expressions as well as for the coefficients that appear in the expressions. Generally, the 
coefficients are very difficult to calculate. We will not do that either, but will simply write the transport 
expressions, and make some evaluations concerning their application. 


General Transport Expressions 


Our first object is to examine the diffusive transport terms which appear in the hydrodynamic 
equations to assess the importance of diffusive processes in the flow fields of flight objects in severe 
thermal environments. The general Chapman-Enskog expressions used include three vectors: the 
macroscopic gradient vector, the mass flux vector, and the energy flux vector (ref. 2, pp. 483, 485, 489, 
516, and 522). 


The macroscopic gradient vector of species i is defined as 


dr dr p 


v 

-^Xi-^kXk 

m; . . 


( 1 ) 


k=l 


The mass flux vector in a v component gas mixture is the mass flux of species i relative to 
coordinates moving with the mass averaged velocity. Its expression (general case) is 

Ji = — X “ D I = ni m iYi = piYi 

P “J dr 


( 2 ) 


(note that Du = 0) 

The energy flux vector is (the subscript c for "conductive" will be omitted for q in this chapter) 


PRECEDING page 
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The dj expression contains three kinds of terms: concentration gradient, pressure gradient, and 
external force. It appears in the mass-flux vector multiplied by a multicomponent diffusion coefficient 
Djj along with a thermal-diffusion term containing the thermal-diffusion coefficient D| Thus, four 
kinds of diffusive transport (concentration, pressure, forced, and thermal) are represented in the mass- 
flux vector, the diffusive velocity Vj, and the energy flux vector q. 

The nomenclature is 

Xi = mole fraction of species i 

ci = mass fraction of species i 

Xi = electrostatic force on particle i 

Dy = multicomponent diffusion coefficient of species i in species j 
Dj = thermal diffusion coefficient 

Yi = diffusive velocity of species i relative to mass average velocity 
X == "modified" translational thermal conductivity 

Importantly, (1) Dy is a function of p, T, and local composition. Thus it is very difficult to evalu- 
ate for more than a ternary mixture. (2) Moreover, for argon, at 1 atm at 17,000 K, X differs by only 5% 
from the ordinary X\ the "ordinary” translational thermal conductivity. 

Let us now use the basic relations, equations (1)— (3), to derive some transport expressions that are 
in common use, and note their limitations. 


Effective Thermal Conductivity 

The assumptions are pressure and forced diffusion are neglected; dp/dr = 0, the gas mixture is neu 
tral, in chemical equilibrium and the elemental composition is fixed. Equation (1) becomes 


di 


= axi 

8r 


(4) 


Equation (2) becomes 


Ji = 



mimjDjj 


axj._D[aT 

dr T fa 


Under the assumption of fixed elemental composition and chemical equilibrium (ref. 5) 


( 5 ) 


Xi = Xi(P.T) 

Differentiation of the mole fraction by the chain rule yields 

3r \ 5T l p dr [dp jj dr 
Substituting equation (7) into equation (5) gives 




p,-. 


3T dr dp dr 


\ t 

2.3T 

T 0 r 


Since dp/0r = O 


Ji = — V mjmPil 

P j=l 4 

Substituting equation (9) into equation (3) leads to 


I Di_dT 
dT dr T dr 


q = - 


h - y hj 

i=l 


n 2 V 3 x j _ x' Di dXi 

7 - 2, m 1 m J D ij -J---j- + nkT ^ ~-L- — 

P j=i dT T £? n ‘ m ‘ dT 


3T 

9r 


The brackets can be called an effective thermal conductivity . Thus 

- dT 

q = — A-eff 

dl 


which is a Fourier conduction format with 


(6) 

(7) 


( 8 ) 


(9) 


( 10 ) 


(10a) 


^-eff = ^-eff(P>T) (11) 

Although the application has limitations noted above, this approach has been very useful, as will be 
apparent when the equations are solved subsequently. It is noted that 


Cj = Ci(p,T) 


( 12 ) 


h(P»T) = X Cjhj 
i=l 


(13) 
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(14) 


c p = 


9h 


V / 




dh i\ . V u l d °i ] 


\0T/ p i= i \dTjp “ "\dTj 
or if c p represents the "frozen" specific heat 


C P = 


I — Cp + ^ hij 

\9T/p w \3T. 


Partially differentiate equation (13): 


Since dp/9i is zero, we obtain 


Then equation (11) becomes 


3h/3h] 3T + /ah| 3p _/ah\ dT 
dr \3T/ p dr [dp) T dr [dTL dr 


dh _ 0T 

3r Cp 3r 


Xeff^h^ p ah 


We define an effective Prandtl number as 


°p a r ^effar 


Pr eff = ^ 
^eff 


which will appear subsequendy as well. 


(15) 


(16) 


(17) 


(18) 


(19) 


Binary Mixture Approach; Fick's Law 

The assumptions are forced, thermal, and pressure diffusion are neglected (dp/dr s 0, and the gas 
mixture is neutral). This leads to the following expressions for equations (1) and (2): 



D[ = 0 


Ji = mim 2 Di2 — 
p dr 


( 20 ) 
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where D 12 is the binary diffusion coefficient. Moreover 


Xi 



1 = 1,2 


(note that T Xi = 1, so that ^1= and di = -d 2 ) 

£1 3r dr 

2 

M = X XiMi = X 1 M 1 + X 2 M 2 

i=t 

Differentiate equation (21) with respect to the space coordinate 

3Xi _ M dci ci 3M 
dr Mi 3r Mi 3r 

Using equation (22) 


dr 


M dci ci 
Mi 3r Mi 


Mi^- + M 2 
dr 


dll 

dr 


From the note above equation (21): 


3c 1 
~dl 


Mi 

M 



i 


M 2 

Mi 


dXi 

dr 


3c 1 
~dr 


1 jXi 

M 3r 


(M 2 ci + Mi -C1M1) 


3 Ci _ MiM 2 3xi 
3r m 2 dr 


or 


3c 1 

d7 


n 2 mim 2 3xi n 2 dx 2 

— n — = — o mim 2 — 
n 2 m 2 dr p 2 3f 


Substituting equation (28) into equation (20) yields 


t dci 

Ji = -pDi2 — 
df 


This is the binary diffusion approximation, or Ficks Law. It is useful for either equilibrium or 
nonequilibrium flow. Then equation (3) can be written (for no thermal diffusion) 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 
(29) 
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From equation (17) in chapter 4: 


Equation (30) becomes 


5h - 9T 

lT Cp * + £ h, ¥ 


3T = J_ 

& Cp 


3h 


_y h ^£i 

* "s '* 


x 

q = __ 


9h _ 
dr 

i 9h 


Define the frozen Prandtl number 


P? = ^E 
X 


and the frozen Le wis number 


(31) 


(32) 



(33) 


(34) 


(35) 


Le = 


_ pPi 2 gp 


(36) 


Thus for a variable elemental composition (even for nonequilibrium), we have 

P p i2 = l _ M- 
Le Cp pj- 




V 


\ dT -/ 


l ~t 


(37) 


(38) 


This is the energy flux vector for a binary gas mixture without pressure gradients, electrical fields, 
or thermal diffusion. But it allows chemical nonequilibrium with variable elemental composition. Of 
course, the y component of q is of most interest, and thus are the partial derivatives of h and q with 


respect to y. This form of the energy-flux vector has been used extensively from the earliest days of 
spaceflight (refs. 6 and 7). The relationship is often approximated by allowing Le = 1, such that 


q = 


p 3h 

Pr a r 


(38a) 


for a reacting gas in or out of chemical equilibrium. Although it applies to a binary mixture, it has also 
been applied extensively to multicomponent mixtures (ref. 8) that are "sortable" into two classes wherein 
each class has a "common" molecular weight and collision cross section among the species. For exam- 
ple, species "1" may include O 2 , N 2 , C 2 , CO, CN, and NO; while species "2" may include O, N, and 
C — but must exclude such species as H 2 , H, e*, etc. 


The Stefan-Maxwell Approach, and the Bifurcation Approximation 

Generally, the hot gas cap about a hypervelocity flight object is truly "multicomponent" chemi- 
cally. Pressure gradients, temperature gradients, electric fields, and variable elemental composition may 
all be important. If we apply equations (l)-(3) to determine transport phenomena, we find that Dy is 
very difficult to obtain. For any two species i and j, Dy is not a function of T and p alone, even for 
chemical equilibrium. It is also a function of the concentration of all other species (ref. 3). It should be 
calculated at every point — even for a temery mixture, it is an enormous task. It would be desirable to 
calculate the transport fluxes and avoid the complication. The Stefan-Maxwell approximation does that 
(ref. 2, p. 718, and ref. 3, p. 570). Moreover, the more recent "bifurcation" approximation extends that 
concept conveniently. 

Recall equation (2), and consider a binary mixture for the moment. Replace pi with nimj, and the 
right-hand equality of equation (2) such that 


Yi 


= Jl 


mmi 


n 2 v , dT d In T 

X-£Z"Wk-sjjr 


dr 


which in the case of a binary mixture is 


\t n 2 . DT 9lnT 

Vi = Jl_ m 2 Di 2 d2 - nTmT — — 


nip 


Yi = m 2 D 12 d 2 - 

nip 


V 2 = ^-miD 21 |di- 


dr 


Djp 

9 In T 

n 2 mim 2 Di 2 

dr . 


CL 

3 In T 

n 2 mim 2 D2i 

dr 


(39) 

(40) 

(41) 

(42) 

(43) 
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Since £ Jj =0, Di 2 = D 21 , thus 

i=l 


Therefore 


which gives 


mini V, + m 2 n^Y 2 = 0 


0 = 


min 2 m 2 Di 2 

' Djp d In t] 

P 

U 2 

n 2 mim 2 D 12 9 f 


m 2 n 2 mi l ) 12 


dl- 


Djp 9 In T 


n 2 mjm2Di2 J 


di + d2 = (d[ + d|) £ ii£l 

n 2 mim 2 D 12 9r 


from the definition of di (ref. 2, p. 470) 


Idi = di+d 2 = o => d|=-dJ 

i=l 


Subtracting equation (42) from equation (43) gives 

-.2 


y 2 -y 1 = n m * n l D 12 

pn 2 ni 


di 


Dip 

9 In T 

n 2 mim 2 Di 2 

3r . 


n 2 m 2 n 2 , 
P n in 2 


""ft + 


D iP 9lnT 


n 2 mim 2 Di 2 dr 


= — (P,+P2) 
pnin 2 


di- 


^ 2 P 9 In T 


n^im^Du di 


9 In T 


^r J -(Y 2 -v 1 ) = di- 

n 2 D i2 9 r 


n 2 n ip 

( ° 2 ) 

,n 2 miniDi 2 

In2m2/ 


n 2ni 

n 2 D 


CV 2 -Vi)=di- 


i 2 


9 In T 

n 2 ni 

( Dj )' 

dr 

_c,n 2 Di 2 

\ri2m2/ 


(44) 


(45) 


(46) 


(47) 


(48) 

(49) 

(50) 

(51) 
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where di includes chemical concentration, and thermal, pressure, and forced diffusion. Curtiss and 
Hirschfelder (ref. 2) have generalized this Stefan-Maxwell relation to multicomponent gas mixtures. 
The following discusses that development. The generalization leads to 


njnj 


j=i n2D ij 


(Yj-Vi) = di- 


9 In T 


ninj 

f D l 

n 2 Dij 

njmj 


n i m i 


( 52 ) 


Note that 


V- — li — 

njmi p . 


h_ 

CiP 


(53) 


and 




di=-X 


XiXj 


j=i pDij L 


t 3 In T t 3 In T 

+ ij + Dj 


9r 


dr 


(54) 


(55) 


from equation (52), where Dij(l) (which is explained in what follows) is written Djj (the binary diffu- 
sion coefficient). 


Now concerning the Stefan-Maxwell equation (52) applied to multicomponent gas mixtures, the 
derivation is lengthy and goes deeply into transport coefficient theory — beyond the scope of this mono- 
graph. The derivation is outlined in chapter 7 of reference 2. Briefly, (1) the function that is used to 
perturb the Maxwell-Boltzmann distribution expression is expressed as a sum of integrals wherein the 
integrands contain infinite series of Sonine polynomials. (2) Only one term in the series is used which 
yields the "first approximation" to the multicomponent diffusion coefficient. (3) Then by relations 
shown, and steps outlined in reference 2, equation (52) is derived. We should note that equation (52) is 
valid for a multicomponent mixture, but uses only binary diffusion coefficients which are comparatively 
easy to obtain, and are insensitive to the mixture composition. Equation (52) is generally called the 
Stefan-Maxwell relation. It is an equation for Vj or Jj in terms of binary diffusion coefficients and 
replaces equation (2). 

An interesting development by Kendall, Rindal, and Bartlett (ref. 8), called the "bifurcation 
approximation," extends this development and has been in use in recent years. It makes use of equa- 
tions (53)-(55). They have found that over a broad range of particles fry can be represented semiem- 
pirically to remarkable accuracy by the expressions in equations (57) and (58). The development is par- 
ticularly noteworthy because it can be used for mixtures which contain the chemical components of 
reacting air as well as species that were produced by an ablating heat shield in the extreme thermal envi- 
ronment of the gas cap over a hypervelocity flight vehicle. The development continues as follows. 
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From equation (55) with thermal diffusion retained, define Ji and Jj so that 


Substituting in equation (56) 


j=ip!V Cj Ci 

D = D(p,T) 


F; = constant such that Djj = — = Dj; 

FiFj 


_ y Xi F iXj F j ft _ Ijl 
ft pD \ c j c »/ 


Therefore 


Equation (1) becomes 


*UM 

c i Mi 


V 

j 1 V Mp m t m_m_ t 

di JBSrV % * _ V % * 


or 51 p m S 


v - v 


■ M?. CjFj v F jJj FJi y c . 

«ri 1U *-> M. M- “ 1 


pD |_ Mi p' Mj Mi “ Mj 


Multiplying by Mj/Fj and summing over all i 




P D i=l ( ft Mj ft Mj 


and noting that 
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(64) 


XdT = o 


2 li = o and 2 c » = X Xi = 1 


i=i 


i=l 


i=l 


i=l 


Equation (63) becomes 


V fill- PR Y 

]fr M J 


r 1 ! 3 r ri 


3 In p 

dr 


M £L PX i . + 
Fi P m i 


Ma 

F iP 


X n k2£k 


k=l 


v 

We now solve equation (62) for 3xi/3r and substitute equation (65) for 2 Fjlj/Mj . 

j=i 


(65) 


Sr Mi £ Fi sj pD Mi jfj Mj Mj £ Fj| P £ k - k Sr J P m i| 


_Cj_ 

Mj 


dr pmj " ^ 


p k=l 


( 66 ) 


Further, define 


j=l h J 


and 


(67) 


Zi=^, ±34-1 

F **2 S 

Then, the chain rule of differentiation gives (from eq. (68)) 

.. azj a^ 2 Mi axi 

M-2 + Zi = 

3f 3f P* dr 


Solving for d^Jdr. 


dXi _ F; 
Sr Mi 


„ ^i . 7 ^2 

M-2 + Zj 


ar ar 


(68) 


(69) 


(70) 


Combining equations (70) and (66), multiplying by Mj/Fi, using X Zi = 1 and X dZi/dr = 0 from 
equation (68) produces 
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3r 3r «i pD^Mj 


* " k=l 


(71) 


where L is the Avogadro's constant. 
Define 


j»I J 


(72) 


We solve equation (71) for ^L=i + Dj T ain T/3r, using MjXj = Ci M: 

(M-Mj) 


i+D Ti£LZ,_J2./| 

9r Mm j 1 


u - dZi . 7 <Hl2 

1*2 + Zj 

«r 5r 


or j=i r q_ 


dr v p k=i 


Cj 

Mi 




3r 


k=l 


Further, define 


(73) 


1*3 = X c" (M - Mj) 
j=l *j 


E 4 = X p - -LpXj + Mj nfcXk 


(74) 


(75) 


k=l 


Equation (73) becomes the bifurcation approximation, generalized here to include all four types of 
diffusion. 


KRB code 


Ji = -D‘ 


T 5lnp pD 


Mpi 


az; duo 

m — + (Zi-Ci)-^- 
dr dr 


L\ 


M-Mj 

Fi” 


-P3 


d_r Fi 


-LpXi + Mj £ nkXk - E 4 1 

k=l 


(76) 


Equation (76) is an expression for the mass-flux vector of species i which removes all j depen- 
dence for thermal and chemical diffusion. The portion concerning transport by chemical concentration 
gradients is labeled KRB and is from reference 8. To this we have added the transport terms associated 
with thermal, pressure, and forced (electrical) diffusion. It is not a pleasant appearing expression, but it 
is profoundly simpler than the approach which calculates multicomponent diffusion coefficients, or even 
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that which uses only one term of the Sonine polynomial series in the integrand. The mass-flux vector 
expressed in equation (76) can be used to replace equation (2), and can be used with equation (1) in 
equation (3) to obtain the energy-flux vector. 


Summary of Transport Expressions 

This is as far as we will pursue the various transport expressions and their approximations. These 
are summarized in table 5-1, with some comment on applicability. It should be said that the comments 
partly reflect convenience. That is, the Xeff approach, though strictly limited to a fixed elemental 
composition, is employed for an elemental composition that is assumed to vary negligibly. Thus, for 
convenience, the properties have been calculated, tabulated, and published for the fixed elemental com- 
position of air without reference to a flow-field solution by Hansen (ref. 5), Yos (ref. 9), Moeckel and 
Weston (ref. 10), and others. Some of these results have been expressed as analytical correlations by 
Viegas and Howe (ref. 11), and subsequently have been used as input to flow-field solutions by a num- 
ber of analysts including Howe and Viegas (ref. 12), implying that the elemental composition was fixed 
or that its variation was of small consequence. The analytical expression for equilibrium air showing the 
normalized ratio of density-viscosity product to the Prandtl number (the latter expressed in terms of 
Xeff) is presented as a function of normalized static enthalpy for three pressure levels in figure 5-1 
(adapted from ref. 11), wherein the correlation expression is a ghastly-looking formula. The circles 
noted as references 1 and 2 on the figure are references 5 and 10 of this monograph. The function will 
be seen subsequently included in the solution of the flow-field equations. As assessment of the sensitiv- 
ity of errors in transport properties as they affect surface heating is presented in reference 13. 


Diffusion Assessment, and Compatibility of a Charged Gas With Linear Flux Theory 

Four kinds of diffusive transport have been considered; that driven by chemical concentration 
gradients, pressure gradients, temperature gradients, and forces arising from induced electric fields 
(corresponding to electric charge separation). Each of these phenomena adds to the complexity of the 
treatment, so it is appropriate to ask, what is their relative importance — can some of them be neglected? 
If so, under what circumstances? Ahtye (ref. 1) showed that thermal diffusion is important in ionizing 
gases, and chemically driven diffusion is generally considered to be important (and is so treated). But 
what of pressure and forced diffusion? Are they negligible? They have not been treated to any large 
degree in the literature. Moreover, Meador and Staton (ref. 14) have argued that the Navier-Stokes 
equations are not compatible with a charged gas, which would preclude charge separation and thus 
forced diffusion as well. We should assess these things. Much of what follows is drawn from Howe and 
Sheaffer (ref. 15). We begin by assessing diffusion caused by pressure gradients relative to diffusion 
caused by chemical concentration gradients for various stages of chemical dissociation and ionization in 
some regions of the flow field. Then electrically forced diffusion is assessed in an approximate manner. 


Pressure Diffusion 

Since the pressure and concentration gradient terms appear together in the expression of the 
macroscopic gradient vector, we compare them to evaluate the role of the pressure diffusion. For that 
purpose, the electrical-force terms are neglected in this paragraph. Using equation (53), equation (1) can 
be written as 
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TABLE 5-1.- TRANSPORT EXPRESSIONS 







TABLE 5-1.- CONCLUDED 
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*Our X is "modified translational" conductivity (called X' in ref. 2, p. 491), but it is very nearly equal the translational conductivity (for Argon at T up to 17,000 K). 

















h/h r 

Figure 5-1.- Normalized ratio of density- viscosity product to effective Prandtl number — equilibrium air. 

d /m \ . in-, n;m; \ 3 In p 


^ ln ' n v I 

dr v 

2, n i m i 
\ i=i 


dr 


(77) 


a d /iij i n; 

d i* "" "" 


1 


n i m i ) 

amp 

v I 

X n i m » | 

i=l 1 

dr 


(78) 


We look at the y-component of the pressure gradient with the approximation that the density ratio 
across the shock wave is small, £ <<.„! . From equations (14) and (18) in chapter 2, and from equation (2) 
in chapter 4: 


a p_ 

pu 2 

(79) 

9y 

R 

u s = U 

sin 0 S 

(80) 

p s = p JEJ 2 cos 2 0 S 

(81) 


Behind the shock, the order-of-magnitude analysis of the Navier-Stokes equations applied to the 
shock layer of a blunted cone indicates that the transverse pressure gradient is such that 


3 In p _ y 3p 
9y ^ 9y 


O 


PsU 2 sin 2 e 5 
pJJ 2 cos 2 9 s Rj 


= 0 


/ tan 2 6h | 
» eR I 


( 82 ) 
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With the assumption 


ir £ 

(83) 

Equation (81) becomes 


/ o \ 

0 In p _ tan 2 0b 1 

dy 8 1 

(84) 

In the stagnation region 0b = 0, so that from equations (84) and (78) pressure diffusion is negligi- 
ble (electrical forces are neglected). But in the shoulder region, tan 2 0b = 0(1), and equation (84) yields 

0 In p _ 1 3p _ 0 |l| 
dy p 0y ^8 j 

(85) 

Consider a dissociating gas and the reaction 


A 2 — >2 A 

(86) 

that we will write 


M — >2 A 

(87) 


Then 


From equation (40), with Dj T = 0 


From equation (77) for atoms: 


n = nA + nM 

mM = 2 mA 


V 



p 


0 / n M 
dy\ n 


m M n ^ 

n A m A + 8 


From equation (88) 



Substituting equation (92) into equation (91), 


( 88 ) 

(89) 

(90) 


(91) 


( 92 ) 
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( 93 ) 


Similarly, 



Vm,- . n2mA , Dam 
(n - n A )p 


1 M , n A 


+ — 1 - 


n 


6 \ n Js n 2n - nA 




1. With very little dissociation near the shock 

- 0(e) « 1 

From equations (93) and (95) 


V*,.^Dam 




P 5 


which can be large 


V Hy -! Dam r+. 1 (e) + ( f ) 1] = 2 SS!a ^ 
P l 8 ' 2, gJ 2p5 


which is small . 

2. Half-dissociation gives 

Equations (93) and (95) become 


n A = n M » if = i 


V A , = ^ DAM [-l(l) + l(l)] 
v My =^ DAM [i(i) + i(i)] 


(94) 

(95) 

(96) 

(97) 

(98) 


(99) 

( 100 ) 


Thus the diffusive velocities of atoms and molecules are comparable, and pressure diffusion is important 
to both. 
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3. For nearly complete dissociation, pressure diffusion can be shown to be negligible for both 
atoms and molecules by similar arguments. 


4. For an ionizing gas wherein we consider A — > A + + e' for simplicity, or replacing the symbol 
A + by I for the ion, assess the importance of pressure diffusion as follows: 

From equation (40), with D{ = 0 (recall that Dy = 0): 


Ya= n (itiiDasIi + meD A< de) 
n A p 

(101) 

- 1 = nip ( meE>Ie ^ e + m A D IAdA) 

(102) 

Ye = ^-(m A D eA dA + mjDejd,) 
n e p 

(103) 


From equation (78) we have, for example, 



men 

n e me + nimi + nAm A 



(104) 


Noting that m A » mj » me and for n A » nj ~ n e (slightly ionized gas) we see that pressure gradi- 
ent is not important for d A and di, but is important for de (exercise). 

It can be shown, from what we have seen so far, that pressure diffusion is 

1. Not important in the stagnation region (unless charge separation induces 3p/9y) 

2. Important in the shoulder region for the diffusive velocity of 

a. Atoms and molecules in a slightly dissociated gas 

b. Atoms and molecules in a half-dissociated gas 

c. Atoms and ions in a slightly ionized gas 

d. Atoms, ions, and electrons in a highly ionized gas 


Forced Diffusion 

The electric force on a particle is 

Xi = co,eE 

where 

X = dynes 

e = 1.602x10-12 erg/V (the electron charge) 

E = V/cm (the induced electric field) 


(105) 
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(recall that C0i = number of electron charge on species i) 

We allow charge separation, so that, from electrodynamics 


div £ = — * 0 
V 


( 106 ) 


where 

Pc “ erg/(Vcm 3 ) 

V = 8.855x 10* 7 erg/(V 2 cm), which is the dielectric constant. 

For convenience, useful conversions of physical units are 

Joule = Coulomb-Volt = 10 7 erg 
Farad = Coulomb/Volt 
Amp = Coulomb/sec 

We consider the electric field to be primarily a function of the coordinate normal to the surface, 
and obtain 


OFy P(. 

dy = V 


The boundary condition at y = 0 is 


(Fy)w = 


V 


where F w = surface charge density . The integration of equation (107) gives 


(107) 


(108) 


= — f p c dy' + — 

V Jo V 


(109) 


Assuming a grossly neutral, locally nonneutral gas we write for a strip in the y direction of cross- 
section dA, 


(dA)\j/E y = 0 = dAT w + dA 



Pc dy = 0 


(HO) 


which leads to 


(Ill) 


Pc dy 

From equations (109) and (1 1 1) we obtain 




( 112 ) 


For convenience, we assume that a thin stratum (of thickness A) of charged gas (with average 
charge density p c ) exists behind the shock as shown in figure 5-2. 



Figure 5-2 - Locally nonneutral flow field. 
The field behind this layer is, from equation (112) 


Ey(5 - A) = 



P.dy = -|[6-(8-A)] = -£^ 


(113) 


Since 


Pc = e £ (Dinj 


i=l 


Ey(8-A) = --&-A. 

y V6 


e X ®i n i 

i=l 


(114) 


We consider only A, A + , e', and the ionization reaction 

A — » A + +e~ 

We define the degree of charge separation S for a ternary mixture of atoms, ions, and electrons as 


S = 1 " = n^ (nA+ “ n<r) = S ®i n i 


(115) 


i=i 


Equation (1 14) becomes 
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(116) 


Ey Hv)(s 

From equation (105), the y component of Xi is 


enA + S 


n A*S 


(In order to compare the forced and concentration diffusion terms, we will neglect pressi 
The y component of equation (1) can be expressed as 


H - 3Xi 
■ iy 

dy 


k=l 


Approximate 


^Xi __ Q Xis ~ Xiw 

5v l 8 


dy 

Then, substituting equations (117) and (1 19) into equation (118) yields 


a. .. Xis - Xiw Ci 






k=l 


/ 


.2.5 A 


' vf" A * s ' 


The equation of state is 


p i = P!5I, nimi J8r =ni EI = „ ik T 

Mi Lm i L 


where R is the universal gas constant, and L is Avagadro's number. 
Dalton's law for partial pressures is 


Thus, 


and 


P = X P‘ = nkT 

i=l 


n j _ Ilj _ Xj 
P nkT " kT 


aZ t 

p m, pm, kT 
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(117) 

diffusion.) 

(118) 

(119) 

( 120 ) 
( 121 ) 

( 122 ) 

(123) 

(124) 


i'i i r 


Using the notation 


A — > A + + e~ 


I 


c^jS _ CjnjS _ CjXiS _ I MA XiS 
“p nPf" kT “PM/kT 

And substituting equations (125) and (124) into equation (120) gives 


X is _Xi 


hy 



Mi 


-®i + T 7 5CiSJ e 2 n lS 


M 


5 A 


Y 8 


(125) 


(126) 


Forced diffusion can be neglected if the second term on the right is much smaller than the First. Multiply 
the right side of equation (126) by 5/Xis* Thus the criterion for species i becomes 


e 2$l-QL A £ [Mi -0)i 
\\f kT | g M 1 1 


«!-- 


(127) 


If Xiw/Xis « 1, and nj = Xjn, the criterion becomes 


\j/kT 




^XiS 

5 


« 1 


(128) 


where for purposes of estimating the effect, all quantities are evaluated behind the shock wave (at high 
altitude, large gradients exist behind the shock wave as well as at the body surface). For increasing Xi, 
the criterion is first violated by the electrons (Mj/M « coi) in the first parentheses followed closely by the 
ions. Generally, the most sensitive condition is for the electrons. Note that 


Accordingly, the criterion becomes 


The Debye length is defined by 


JL = _P_ 
kT (kT) 2 


e 2 a 2 P A 
¥ (kT) 2 8 


XiS « 1 


a 2 _ ykT 
e 2 ni 


In these terms, the criterion is 


(129) 


(130) 


(131) 
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(132) 


e 2 ni 5 2 A ^ 
\j/kT 5 


s^a 

A 2 5 


S« 1 


The pressure find temperature for air in chemical equilibrium behind a normal shock is given as a 
function of flight speed and altitude in reference 12. For a given p and T, ft for air can be estimated by 

Xl~X n *+Xo + +XN$ + --- (133) 

For assigned values of S, the locus of the criterion equation (132) is shown in Figure 5-3, which was 
abstracted from reference 15. Thus for 6 = 0. 1 cm and A = 0.5 5, forced diffusion can be neglected 
above the lines for constant S, but is not negligible below the lines. For the trajectories shown, a trace 
of charge separation can affect dg-, dj, and thus the mass flux vectors for atoms, ions, and electrons for 
most of the entry regime. The above argument could be made more rigorous, but it is sufficient to pro- 
vide the information that we wish to know. 



Figure 5-3 - Locus of electronically forced diffusion criterion (eq. (132)). 


Compatibility of a Charged Gas With Linear Flux Theory 

Finally, inasmuch as we are considering nonneutral or charged gases, we should say a few words 
about the compatibility of the Navier-Stokes equations with a charged gas. Several writers have raised 
doubt that they are compatible, and Meador and Staton have offered a proof that the concept is not com- 
patible (ref. 14). The proof argument is presented briefly here, along with my rebuttal (which also 
appears in ref. 15). 

The argument centers around the macroscopic gradient vector (eq. 1), which can be combined with 
equation (105) to yield 
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in i r 


The electrodynamic relation 


-■ E = — = en V 

3r V ft V 

was employed with an assumed relation for the electric current (I here) 


I = aiE+o 2 


9T 


9r 
\ '/ 


The current was set to zero, which resulted in the expression for the field 


E = — 

<*1 


9T| 

v dr, 


This expression was used in equations (134) and (135). The charge density, effectively 


(135) 


(136) 


(137) 


en ^ 0);%; 
i=l 


from equation (135) was used in equation (134), which resulted in an expression of the form 


A - ^ ♦ «i - o ^ + H [*» + v N V 2 t 1 
9r 9r lP a i /L m * l^i / J 


3T 

dr 


(138) 


The last product V 2 T(9T/9r) is "third order" in temperature, which is not consistent with linear flux the- 
ory. Thus, it was argued, (8/8i) E and therefore charge density p c must be zero in equation (135) to 
avoid the third-order expression. 


However, that conclusion seems to be a consequence of the use of two relations where one would 
suffice. There are hydrodynamic equations to calculate every quantity in equation (134) except E. One 
additional relation for E (eq. (135)) needs to be solved simultaneously with the hydrodynamic differen- 
tial equations. Equation (135) is not (in the hydrodynamic sense) an equation for calculating 

v 

CDi%i 

i=I 

because the mole fractions are calculated from the solutions of the diffusion equations, and C0j is a 
set of constants. However, reference 14 used equation (135) to calculate 


X 

i=l 

and used a second relation (eq. (136)) with ]_= 0 to calculate £. Of course, both equation (135) and 
equation (136) are valid and physical, but if both are used in the hydrodynamic framework, equation 
(135) should be regarded as a relation for E, and equation (136) a relation for one more unknown, the 
current, J. Clearly, within the hydrodynamic framework, equation (134) is first-order without requiring 
that pc be zero. The electric field is simply calculable from the electrodynamic differential equation 
(135), just as the vectors d lnp/3i and dXi/dr are calculable from the hydrodynamic differential equa- 
tions for insertion into equation (134). Thus there are no third-order terms and, on these grounds at least, 
a charged gas is compatible with the hydrodynamic equations. 

Alternatively, if it is assumed that there is no current, and the Chapman-Enskog expression for the 
current (eq. (71) of chapter 3) 

V 

T V Wi T 

i = e X-J i (139) 

i=i m > 

(where Ji is the mass flux vector of species i) is set to zero, there obtains for the electric field 

E = A ^- +B S 0J i m Pu^ i +C^ (140) 

where 


and 


X co.mjDijOCj - cj) 
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e cOiirijDjjCj 
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® cOjitijDijCj 

_ ij 
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(143) 
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Substitution of equations (140) and (105) into equation (1) yields 


dr 


di= — + Xi-Ci| 1 +eA(j^ 7 lco i -n X Xk®K 


k=l 


( 8 * 


dr 




(144) 


Again, there are no third-order terms, and the expression is compatible with linear flux theory (Navier- 
Stokes equations). It is self-consistent with the Chapman-Enskog solution of the Boltzmann equation. 
Coefficients A, B, and C can be computed with that theory. As before, the mole fractions of all species 
are obtained from the solutions of the diffusion equation (eq. (5) of chapter 4), rather than from equa- 
tion (134). For the zero current approximation, equation (135) is superfluous, and can be omitted. 
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CHAPTER 6 


RADIATIVE TRANSFER; RADIATING SHOCK LAYER 


In the development of terms that appear in the governing flow-field equations of chapter 4, we 
come to the final formulation, that of radiative transfer. In a gas mixture at high temperature, various 
atomic and molecular processes cause the gas to emit and absorb radiant energy. As we said at the out- 
set, we will not delve into the internal events of these species — that is a subject in its own right. Aside 
from a superficial description of these internal processes, we will confine our attention to the transfer of 
radiation in and through this gas mixture and assume that the coefficients that we need to quantify emis- 
sion and absorption of radiation are available to us. The physical processes have to do with the transi- 
tion from one energy level to another in a gas species accompanied by the emission or absorption of 
radiation (a photon of radiant energy at a particular wavelength or frequency); emission if the species 
drops to a lower energy state, or absorption if it is excited to a higher energy level. Thus there are rota- 
tional, vibrational, and electronic states of molecules and atoms that are affected. Moreover, for electron 
processes, there are free-free, free-bound, and bound-bound processes which correspond to a trajectory 
change of a passing electron, capture of a free electron, or the transition of an electron from one orbit to 
another within a species, respectively. These various processes give rise to photons in specific parts of 
the electromagnetic spectrum, including spectral lines, bands, and continuum radiation. The details of 
the physical processes are beyond the scope of this monograph. 

There are two phenomena that will be neglected because they are not important to the application 
to hypervelocity flight. First is the internal energy of radiation on the macroscopic level. The phe- 
nomenon of the instantaneous presence of photons of radiant energy in a volume of gas does not affect 
the internal energy of the gas at temperatures that concern us (up to 20000 to 30000 K). Secondly, radi- 
ation pressure will be neglected for our application. It is important for very large regions of very hot 
gas, but can be neglected in the equations of motion for present purposes. These two phenomena are 
discussed in references 1 and 2. 

The phenomenon of radiative transfer in a hot gas cap is illustrated in figure 6-1, which shows a 
disk that is 10 in. in diameter mounted normal to a hypervelocity test stream of air. The air in which the 
model is immersed is heated by an electric arc discharged into the air upstream of the model, after which 
the air is passed through an aerodynamic nozzle having an exit diameter of 3.5 ft. The electric arc used 
can operate at 60 MW to produce enthalpy levels up to 20,000 Btu/lb m on a continuous basis, and can 
accommodate models up to 20 in. in diameter. The test conditions for this figure correspond to a less 
severe environment intended to simulate flight at an altitude of 63 km and at a speed of 4.3 km/sec 
(recall that orbital speed is about 7 km/sec near Earth). The gas between the bow shock wave and the 
disk is luminous in the figure because it is emitting and absorbing radiation, some of which escapes the 
gas cap. (This figure is by courtesy of Warren Winovich and John Balboni of NASA Ames Research 
Center.) 

Thus the figure illustrates the phenomenon which we intend to describe. Our objective in this 
chapter is to provide expressions for the radiative heating rate, 7f r , and the term, div qf r , used in the dif- 
ferential equations that describe the flow field about a body in hypervelocity flight (eq. (3) of chapter 4, 
the energy equation). Eventually, we will solve the set of coupled integro-differential equations. 

For present purposes, a more detailed development of that which follows can be found in refer- 
ences by Chandrasekhar (ref. 3), or a very readable account by Kourganoff (ref. 2), wherein a very 
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Figure 6-1.- Radiating gas cap. 


Some Definitions 


detailed treatment of the development can be found, including some subtleties that are not elaborated 
upon in this monograph. 

As shown in figure 6-2, let P be a fixed point in an absorbing and emitting gas mixture and let S 
be a fixed line through P. Let da be an element of area containing P and let y be the angle between 
S and the normal Y' to da. 


The amount of radiant energy flowing through da, in a specified frequency interval (v, v + dv), in 
a direction making an angle y with the normal to da, within an elementary solid angle dco' about the 
direction S, in time dt, is 


dE v = I v cosy dadco'dvdt ( 1 ) 

It follows that the spectral intensity I v can be defined as the energy flowing at the point P in the 
direction S, per unit of time, of frequency interval, of solid angle, and of surface area normal to S (what 
is fundamental is not da, but its projection da-cos y). In the general case, 
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Figure 6-2- Radiation coordinates. 
Iv = I v (P,S,t) 


( 2 ) 


It may be noted that for all S directions that make an angle y with Y' the solid angle dco is 
(2k sin ydy). The differential spectral radiation flux due to contributions from one S direction making 
an angle y with the normal Y' is 


dF v = = I v cos y dco' 

da dt dv 


(3) 


Summing over all S directions making an angle y with the normal Y' (replace dco' by dco) and 
over all angles y yields the directional spectral radiation flux 


Fv = ( Iv 
Jy= 0 


cos y(2tc sin y dy) = 2k I I v cos y sin y dy 

Vo 


r 

Jy = ( 


(4) 


The interaction between radiation and the gas is expressed in terms of an absorption coefficient and 
an emission coefficient. With reference to figure 6-3, the intensity I v in the direction S becomes 
I v + 8I V along S' by interaction between the radiation and the matter dm contained in the volume 
da ds. Experiment and theory show that 


= -K v p ds = -|i v ds 
Iv 

where K v is the mass spectral absorption coefficient, 1 and radiative scattering phenomena are 
neglected. 


(5) 


Similarly, an element of mass dm sends in directions confined to an elementary solid angle dco', 
in the frequency interval (v + dv) and in time dt, an amount of radiant energy equal to 

dey = J v dm dv dco' dt (6) 


1 In the strict sense, the absorption coefficient, K v , deals with the Lagrangian loss of photons (ref. 2, pp. 4-5). 
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and is composed of that energy truly emitted by (created by) and that scattered into (from all directions) 
by dm (J v is the mass spectral emission coefficient). 

Denoting the intensity of black radiation of frequency v and temperature by B V (T), the local 
thermodynamic equilibrium (LTE) assumption is one of Kirchoff s laws 

J v = rl B V (T)K V (7) 

where Ty is the index of refraction of the medium. For simplicity, we restrict ourselves to media for 
which r v = 1. Thus 

J v = K v - B V (T) (8) 

for LTE. The explicit form for the function B V (T) is given by quantum statistics, and is called Planck's 
law: 


BvCD 


2hy 3 

1 ] 

C 2 

e hv/kT _ l 


(9) 


where c is light speed, h is the Planck constant, and B V (T) is also called the "source function." 


The Equation of Transfer 


So far, we have adopted a partially Lagrangian point of view and, by means of the absorption 
coefficient and the emission coefficient, have examined what happens to certain photons. From the 
Eulerian point of view, the variations in the intensity of the photons near a given point are now 
examined. 

Accordingly, consider the small cylinder element (shown in fig. 6-3) of cross-section da and 

length ds in the medium, and the radiant energies which cross the two faces normally in a time dt and 

in directions confined to a solid angle dco'. 



Figure 6-3.- Volume element of gas oriented in the S direction. 
The energies in and out are 

dE Vin = I v da da>' dv dt 
dE Vout = (I v + dl v )da dco' dv dt 
The energy balance can be written as 

-(dE Vta - dEv^,) = J v dm dv dco' dt - (K v I v p ds) (da dco' dv dt) 


( 10 ) 

( 11 ) 

( 12 ) 
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Noting that 


dm = p ds do 

we substitute equations (10) and (11) into equation (12) 

dl v dv da dco' dt = J v p ds da dv da)’ dt - K v I v p ds da dv dco' dt 

This yields a first expression of the equation of transfer in the S-direction 

_ 1 dly(s) _ y , v _ Jy 
pKy ds v( } Ky 

which for LTE becomes 


1 dly(s) 
pK v ds 


— I v (s) — By(T) 


(13) 


(14) 


(15) 


(16) 


Plane-Parallel Radiating Layer of Gas 

For our purposes, a blunt stagnation region or a stellar atmosphere is considered to be stratified in 
homogeneous plane-parallel layers in such a way that a single variable (geometrical depth y) is suffi- 
cient for the specification of a layer of the atmosphere. The atmosphere is supposed to be in strict 
radiative equilibrium, i.e., it is assumed that the heat interchanged by convection and conduction is neg- 
ligible compared with that interchanged by radiation. 

We shall therefore introduce the geometrical depth y and the direction parameter w = cos 0 (see 
fig. 6-4). 




EDGE OF STAR 
(OUR WALL) 


Figure 6-4 - "Plane-parallel" radiating gas. 

dy = ds cos y = ds cos(7t - 0) = -ds cos 0 (17) 

ds = — = -dy sec 0 s - J- dy (18) 

cos 0 

w = — 1 — = cos 0= -cosy (19) 

sec 0 

As before, the assumption of local thermodynamic equilibrium, J v /K v = B V (T), is made (eq. 8). Thus we 
obtain the modified equation of transfer in any S-direction in terms of y: 
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( 20 ) 


where 


JBL 


dl 


pK v dy 


*=I V -B V (T) 


I v = I v (w,y) ; p = p(y) ; K v = K v (y); T = T(y) 


( 21 ) 


The Intensity I v Expressed as a function of B v and the Optical Depth 
We introduce the function T v called the spectral optical depth 

■y 


■f 


T v = I pK v dy 

Jo 


( 22 ) 


thus, the element of optical depth 


dx v = pK v dy 


(23) 


represents the relative weakening of the intensity transmitted, normally, by the corresponding layer of 
the geometrical thickness dy. Then equation (20) becomes 


w^I v -B v (T) 

QTv 


where 


(24) 


I v = I V (T V ,W) 

We first want the solution of the homogeneous equation 

— = dx v 

Iv w Ulv 


(25) 


(26) 


This is 


Iv 



(27) 


where C is a constant. The particular integral of the nonhomogeneous equation is obtained by substi- 
tuting equation (27) into equation (24) (where C is regarded as a function of x v ). From equation (27), 


Thus 



(28) 


Cex pfe)+ w 



(29) 
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finally, we have the relation 


dC _ B v(T) 
dTv w exp (^) 


(30) 


Integrating equation (30) gives us C as 


C 



Bv(T)«p(-JL)A 


(31) 


where the variable t replaces t v . 

From equations (27) and (31), we deduce the relation we were looking for: 

i , = -' x p(w )| B v ex p(-w)w=-| ^''pj-^rjw < 32 > 


Intensity Direction Convention and the Radiative Flux Expression 

The determination of the constant C brings in the distinction between I\0( v , + ) (toward the wall), 
and Iy (away from the wall). 

We define Iy for 0 < w < 1 by switching limits, changing sign, and allowing C — » thus for 

rays directed toward the wall from layers between t v and «>, from equation (32). 

IJ=| B v (T)exp(-^)dt (33) 

Ay 

Moreover, we know from equation (22) that at the surface (t v = 0), rays directed away from the wall are 
from equation (32) and our convention 



But at t v , the rays directed away from the wall are from equation (32). 






(34) 


(35) 


therefore, at T v , 1^ is the sum of the contributions from the source functions corresponding to the layers 
below t v , each weakened according to its optical distance (x v - 1) from location t v , and of the surface 
term Iyo- 
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Writing the last integral of equation (35) first, the expression for rays directed away from the wall 
at the location t v is (by use of eq. (32)) 




ly(Xv) = (Ivo) ex P| w I - I B v (T)eXf 
Recall from equation (4) that 

I v cos y sin y dY : 


(36) 


F v = 2k (" 

Jy=0 


0 


-2„r 

Jn/2 


(I v )cos y sin ydy=2n I (I+)cos y sin y dY (37) 


But 


and 


cos Y = -w , sinYdY = dw 


cos y sin y dY = — w dw 

Substituting into equation (37) using equations (36) and (33) and noting that when y = 0 
w = -1, yields 1 


(38) 

(39) 


Define m by 


Thus when 


II 

- f BvCDex^- i^) * 

w dw 

r i 




-2tc 

Bv(T,exJ-‘ w T ”k 

w dw 


Jo 

A v 1 1 




m = -L , dw = - t w <j w - _ dm 
m 2 m 3 


(40) 


(41) 


w = ±1 => m = ±1 
w = 0* => m = ±oo 

Correspondingly, substitute m for w and change the limits in equation (40), 


( 42 ) 

( 43 ) 


68 


! I TP 


F v = -2 tc 


f 


Og - f MD^PHt-TvKm)] dt J Mm) 

Jo 


/T 


-271 


B v (T)exp[-(t - Xv)(m)] 


dt 


nr 


(-dm) 


( 44 ) 


where the first integral originates at the wall, the second integral is flux outbound from the gas between 
the wall and x v , and the last integral is inbound flux from the gas between x v and «>; and the entire 
expression gives the net spectral flux at the location x v . Interchange the order of integration, integrating 
over m first, noting that B V (T) is not a function of m, 


F v = 2n 




exp[-(t - x v )(m)] 



dt 


+ 2tc I B V (T) J I ■ 

J X v I J -F°° 


expKt - x v )(m)] 


dm 


m z 


dt 


(45) 


We introduce integro-exponential functions of order n (an integer), which are defined by 


En(Q = I 


m" 


(46) 


In the integrals of equation (45), we replace m by -m and dm by -dm, which yields 


[ - exp(x v m) 

j, 


J i m 


(47) 


and 




exp(-(t - x v )m) 


dm = - 


m z 


f 


exp[-(x v -t)(m)] 


dm = -E2 (x v - 1) 


mr 


(48) 


We rewrite equation (45) as 
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p _ 0 f - exp(-x v m) f Jy [*” 

Fy ” 2n I (V 5 dm + 2 tc I B v (T)E 2 (x v - t)dt - 2n I B v (T)E 2 (t - x*) dt (49) 

Jl W ; 0 A v 

We define the radiosity of the wall (which is assumed to be diffuse, i.e., the intensity does not 
depend on the direction) 


J vo = Kl v (0) 


Then, the first term of equation (49) becomes 


• +oo 


2k | (1^) CXP( Tvm) dm = 2J V0 E 3 (t v ) 

m 


(50) 


(51) 


We now assume that the wall is opaque with emissivity e v and reflectivity R v 

, 4 - 

Jv 0 = e v KB v (T w ) + R v 2k I B v (T)E 2 (t)dt 

Jo 


(52) 


last integral of equation (49) 
with t v = 0 


Substitute equation (52) into equation (51), then into equation (49) integrated over v, with the notation 
Qr 


■f 

Jo 


q r = I F v dv 

/o 


(53) 


we find at x v , the radiative flux is 

wall emission wall reflection 


■1 


Or I 27ce v Bv(T w )E 3 (x v )dv + I 2 R v E3(Tv)) 


) 2R v E 3 (t v )) [ 

Jo Jo 


27tB v (T)E 2 (t)dt 


dv 


IT 


+ 2 k | | B v (T)E 2 (t v - t)dt dv 

/o 


Poo Poo 

- 2 "I j 

Jo Jt v 


B v (T)E 2 (t - x v )dt dv 


(54) 


gas layer below x v gas layer above t v 

That is, for a plane-parallel atmosphere, equation (54) is the radiative flux at x v (y) in the +y-direction. 
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Grey Gas, Grey Surface 


The grey gas assumption implies that the radiation is not a function of wavelength, i.e., properties 
such as K v and I v are independent of the frequency. The grey surface assumption yields e v = £ (which 
here represents the grey surface emissively), and R v = R (which here represents the grey surface reflec- 
tivity). Thus to simplify our development to a grey radiation problem. 


I. 


BCD = B v (T)dv 


(55) 


Then, the grey flux in the +y-direction is (where K and x are independent of frequency, v, and t 
is a dummy x(y)) 


f 

Jo 


q r (y) = 27teB(T w )E 3 (x) + 47tRE 3 (x) | B(t)E 2 (t)dt + 2ic | B(t)E 2 (x - t)dt - 2 ju | B(t)E 2 (t -x)dt (56) 
where 


■f 

Jo 


■f 


The divergence of the radiative flux is 


B(T) = B(t) 


div q r = ^- q r (y) 
- 9 y 


Recall that 


1 


and 


x = I pKdy 


dx = pK dy , •$-= P K — 

K J dy K dx 


We differentiate E 2 and E3, by use of equation (46) to obtain 

dE 2 (x - 1) 


dx 

dE 2 (t - x) 
dx 


= -Ei(x-t) 
= E,(t-x) 


(57) 


(58) 


(59) 


(60) 


(61) 

(62) 
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dE 3 (X ) 


( 63 ) 


We now differentiate equation (56) with respect to y, using equation (60), and assuming that the 
wall is black, e = 1, and R = 0. First write 


■f 


q r (y) = 27tB(T w )E3(x) + 2n I B(t)E 2 (x - t)dt - 2rc I B(t)E 2 (t - x)dt 


f 


By Bx By 


1 

f T 

-B(T w )E 2 (x) + 

B(x)E 2 (0) + j B(t) dE2 ^ _t) dt 

Jo dT 


-B(x)E 2 (0) + j B(t)E 2 (t - x)dt 
Using equations (60)-(63) and E n (0) = l/(n - 1) (for n > 1), 

-B(T w )E 2 (x) + 2B(x) - f B(t)Ei(x - t)dt - I B(t)Ei(t - x)dt 


— = 2?cpK 

ay 


f 


Combining the integrals. 


aq r 

By 


= 27tpKl 


-B(T w )E 2 (x) + 2B(x) - I B(t)Ej(lt - xl)dt 


f 


(64) 


(65) 


( 66 ) 


(67) 


Finally, we use Stefan's Law, 7tB(T) - oT 4 (where o is the Stefan-Boltzmann constant) and the 
fact that there is no radiation originating beyond x = x s (at the shock): 


3q r 

— = div q r = 2opKj 


J /»Tsor« 
0 


-TtE 2 (x) + 2T 4 - I TtOEKIt -x l)dt 


( 68 ) 


Equation (68) is the term to be used in equation (3) of chapter 4. 

Transparent Gas Approximation for the Emission of Radiation 

It is convenient for emitting gases which absorb little, or for purposes of making radiation esti- 
mates, to consider the extreme case of a transparent gas (fig. 6-5). 



Figure 6-5 - Transparent emitting gas cap. 

The radiative emission rate per unit mass, per unit time, in the frequency interval (v,v + dv), within 
the solid angle do)', is 


or 


de v = J v dm dv dco' dt = J v p dV dv do)' dt 


(69) 


dVd[ = J vP dv dco' (70) 

In all 4k steradians, the gas emits radiation at the rate per unit volume 

4tcJ v P dv (71) 

Over all the frequencies, the emission rate per unit volume is 


E t = 4;tp I J v dv = 4jtp | K v B v (T)dv 

/ o 

We define the Planck mean mass "absorption" coefficient as 


r r 

I J v dv = 4jtp I 

Jo Jo 


(72) 


K p = 


J K V B V dv 71 j 

_o Jo 

[\ dv " 

Jo 


K V B V dv 


(73) 


Then, the energy emission per unit volume is (K p is obtained from tables) 

E t = 4oT 4 K p p = div q r (transparent gas) (74) 

In the one-dimensional approximation, half the radiation is toward the wall and half is directed 
away from the wall, thus 
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(75) 




E t (y')dy ,+ 


where the last term is the outward emission from a radiatively black wall. At the wall, the flux is 


qr(O) = aT t-^ j E t (y)dy (76) 

and if the wall emission is neglected, 

qr(0) = ~ 2 J E t(y) d y (77) 

A combination of shock tube-experiments and theory by Kivel and Bailey (ref. 4) estimated the radiant 
emission rates per unit volume of air as a function of temperature and density. The air was assumed to 
be in chemical equilibrium with the corresponding chemical species concentrations. These results were 
correlated by Howe (ref. 5) in the form 


div q r = E t = c 



(78) 


The correlation is shown by the solid lines in figure 6-6, where the dashed curves are the results from 
reference 4. The ordinate is in terms of Et/2 because it was used to estimate the heating rate to the wall. 
Thus Et/2 multiplied by —5 is an estimate of the radiative heating to the wall from a transparent gas at 
uniform temperature in accord with equation (77). The correlation expression (eq. (78)) was received 
with more graciousness than it deserved (ref. 6, p. 398). 


The Nearly Opaque Gas Approximation for the Transmission of Radiation 

The other extreme, with respect to the transparent gas approximation, is the strongly absorbing gas 
that is nearly opaque. The latter approximation leads to a flux-divergence term in the energy equation 
which has the appearance of a radiation "diffusion" term, rather than the complete integral radiative- 
fransfer expression. The corresponding differential equation is much more appealing than the full 
integro-differential equation. The approximation is called the Rosseland diffusion approximation. 
Although it appears in the literature, and a terse description of its derivation can be found in reference 7, 
I am unaware of a meaningful derivation in the literature. Thus a brief derivation is given in what fol- 
lows (see fig. 6-7). 


mi ■ n 
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Figure 6-6 - Correlation of radiative emission rate for equilibrium air. 



Figure 6-7 - Three-dimensional orthogonal coordinates. 

For any point in the gas that does not "see" the wall, the first two terms of equation (54) can be 
omitted. For present purposes, consider the gas in a three-dimensional, orthogonal coordinate system 
shown in figure 6-7. Any point x,y,z, receives a portion of the radiation emitted from other locations 
x',y',z' a distance L away. We have defined |i v = pK v (eq. (5)), and J v (eq. 7), the mass spectral 
emission coefficient, or its equivalent K V B V (T) (eq. (8)). It is now convenient to define 
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By — 4|C|XvBy 


which is the spectral energy emission rate per unit volume. Thus in the absence of radiative scattering, 
the x component of the spectral flux at x,y,z which originated at x',y',z' is 


F Xv (x,y,z) = - 


ByCx'^z') 

4;tL 2 


cos(x,L)exra-| p v dL' dx' dy' dz' 


where the brackets are the intensity at radius L without absorption, the cosine term gives the x com- 
ponent, and the exponential factor accounts for absorption. Also 


cos(x,L) = 


and as mentioned above 




ex \ J HvdL ; (82 > 

is the absorption along the path L. For this special case of strong absorption, we assume that the tem- 
perature varies slowly spacewise, and that radiative equilibrium exists locally. The local energy density 
of radiation at the frequency v is characterized by a black body at the temperature T and is (see ref. 2, 
eq. (2.8)) 

uv -7 ^ r (83) 

(c 3 [exp(hv/kT) - l]j 

The energy absorbed at that frequency per unit volume per unit time is 

cuvfiv (84) 

which is equal to the rate of energy emission locally. Thus 

B v = cuvpy (85) 

Assume that along a path length l/|iy, Uy|iy changes little for all frequencies, so that equation (85) 
characterizes T(x ,y',z') locally. The variation of |i v is small with distance, and py is kept constant in 
the integrals of equations (80) and (82). Substitute equations (85) and (81) into equation (80) to obtain 


F Xv (x,y,z) = -cp v 


0 (x'-x) 


e~^L dx' dy' dz' 
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Expand u v in a taylor series 


u v (x',y',z') = u v (x,y,z) + 


9u v 


\ 


dx 


( ' u ^ Uv 
(x - x) + — 

x,y,z \ a y I x,y,z 


(y* - y) + 


du v 

dz 


(z'-z) 


(87) 


fx.y.z 


Substitute equation (87) into equation (86) to obtain 

rrr 


Fx v ( x >y>z) = -c^v 


JJJ- 


u v (x,y,z) + 


du v 


y ^X Jx, y,z 


( x '~ x ) + 


du v 

. dy 


(y'-y) 


fx.y.z 



(x'-x) 


exp(-p v L 


j dx'dy'dz' 
1 471 L 2 


( 88 ) 


The terms that are underscored integrate to zero because they are antisymmetric with respect to (x' - x), 
etc. Thus 


Fx v (x.y»z) = -c|i v 



(x' - x) 2 exp(-pvL) 
4 jiL 3 


dx' dy' dz' 


(89) 


Since 


cos(x,L) = cos 0 = — — — 

L 

this becomes 


(90) 


Fx v — c|i v 


3u v 


9x 


fx,y,z 


i 


^ e~^ vL dx' dy' i 
4 tcL 


(91) 


Now replace the volume element dx dy dz by an element in spherical coordinates centered at x,y,z as 
shown in figure 6-8. Thus the volume element is 


(L sin 0 d<|))(L d0)dL = L 2 sin 0 d0 d<t> dL 


(92) 


and 
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cos 2 0 
4 jiL 


L 2 sin 0 d0 d<> dL 



VOLUME: (L sing d0)(L dfl)dL 
OR L^sinS d 6 dtfdL 



Figure 6-8.— Spherical coordinates and volume element. 
We integrate with respect to 0 first, and then with respect to <|>. 


and finally 



F Xv (x,y,z) 


-city 

3 



\ /x.y. 




w 

— 

3ji v 

i 9x / 


Since u v is a function of T alone 


Fx v (x,y,z) = - 


c du v dT 
3^v dT 3x 


Thus the total flux component in the x direction is 


( 97 ) 



Rosseland defined the "mean free path," Ar, such that 


F x 



dErad 

dx 


or since 


(98) 



(99) 


where c is the Stefan-Boltzmann constant. Then 


F x =- I ^^rT 3 — (100) 

3 dx 

The y component is obtained by replacing x by y, and the divergence of Fy or q f y is simply 
(3/3y)(q ry ). Of course, Ar can be derived from standard tables of data by the expression obtained by 
combining equations (97), (99), and (100) 


Ar = 



1 dUy 
|iv dT 


dv 



( 101 ) 


and has the units of length. The approximation is somewhat limited in flow-field computations because 
it tends to break down near the surface of the vehicle where temperature is diminished to the extent that 
Ar becomes large; that is, the gas becomes less opaque. It could be terminated arbitrarily at a distance 
from the surface that is some factor of Ar. 

Before proceeding to use radiative-transfer expressions in the flow-field differential equations, it 
should be noted that in the treatment which follows, the gas will be assumed to be grey. That is, the 
dependence on wavelength will not be included. That is for convenience and not because of necessity. 
As a practical matter, a single integral over physical space rather than a double integral which includes 
wavelength appears in the integro-differential equation for energy. The latter integral could be included 
with a modest increase in complexity. Our treatment uses the Planck mean mass absorption coefficient, 
equation (73), where K v weighted by B v was integrated over the entire spectrum, which smoothes out 
any strong or weak absorption effects. It is simply denoted by K in the next chapter. Another approxi- 
mation which considers portions of the spectrum to be transparent, and other parts to be strongly 
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absorbed (ref. 8, for example) is sometimes called the "picket fence" approximation. It is a simple spec- 
tral form and is easy to use; it may be particularly useful for ablation products which, although cold 
absorb strongly in some portion of the spectrum. 


Having mentioned spectral features, it is useful to mention some terms that have been used in this 
chapter, that is, a single photon of radiation at a particular frequency or wavelength 

AE = hv = ^ (102) 


here Planck’s constant is 


h = 6.6252x10 27 erg sec (103) 

The speed of light is c, and he = 1.9862x10-12 e rg p. An angstrom is 10-8 cm or 1(H p. Thus the 
units of a photon can be shown to be ergs, or can be expressed in terms of electron volts, a unit of energy 
corresponding to a particular wavelength. We should also mention the particular wavelength corre- 
sponding to the peak in the black-body function By The peak is a function of temperature and the 
expression is called Wein's Law 

X. m = 2890/T (micrometers) 

Thus for a temperature of 14000 K, 

K = 2890/1. 4X10 4 = 0.206 urn 

and the photon energy can be calculated to be 6 EV at that wavelength by use of equation (102). 


(104) 

(105) 
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CHAPTER 7 


SOLUTIONS OF THE RADIATING SHOCK LAYER 


The differential equations listed in chapter 4 can now be solved. Air in the hot flow field is 
assumed to be in chemical equilibrium. The structure of the flow field will be shown and quantities 
derived from the solutions will show radiative and convective heating as well as the effects of transpira- 
tion or mass addition. What follows is largely drawn from Howe and Viegas (ref. 1). Thus the equa- 
tions for the nonadiabatic shock layer in the stagnation region of bluff bodies are solved for flight speeds 
up to 50,000 ft/sec in air. The effects of energy transport by conduction, diffusion of dissociated and 
ionized species, and gaseous radiation (with reabsorption) are included in the analysis. The effects of a 
foreign species on the radiant energy transfer in the shock layer are also investigated. Convective and 
radiative heat-transfer rates in the stagnation region of the body are obtained from the solutions and are 
compared with the results of others. Coupling between the two modes of heat transfer is examined. A 
simplified method for predicting the effect of radiative transport on convective heat transfer is discussed. 
Shock standoff distance determined from the solutions is presented and compared with the results of 
other investigations. Results indicate that coupling among the energy transfer processes may reduce the 
heat transfer by as much as 50%. 

At flight velocities above about 30,000 ft/sec two phenomena may importantly affect the aero- 
dynamic heating rate in the stagnation region of bluff bodies. The first of these is the ionization of the 
dissociated air, which may affect heat transfer because it increases the total thermal conductivity of the 
air. The second phenomenon is the emission and absorption of radiant energy by the air in the gas cap, 
which for large blunt bodies may impose a much larger heating load than convective processes do. 

Each phenomenon has been treated separately by others. The effects of ionization on convective 
heat transfer have been examined in boundary-layer analyses by Hoshizaki (ref. 2), Cohen (ref. 3), and 
Pallone and Van Tassel (ref. 4) for equilibrium air; Adams (ref. 5) for frozen flow; and Scala (ref. 6) for 
equilibrium nitrogen. All of these investigations except reference 6 predict that ionization will have a 
modest influence on convective heat transfer. Stagnation-region heating caused by radiative emission 
from the isoenergetic equilibrium shock layer has been studied by Kivel (ref. 7). 

At higher flight speeds, it can be expected that a coupling exists between convective and radia- 
tive heat transfer, and that the coupling is more pronounced for larger bodies. For example, the convec- 
tive heat transfer to the stagnation surface can be expected to be diminished if the shock-layer enthalpy 
is diminished by radiant emission. Also, if radiation produces enthalpy gradients in the shock layer, 
thermal conduction and diffusion of species may become important in that region. 

As noted previously, for purposes of analysis, it has been customary to separate the flow field in 
the stagnation region of blunt bodies into a viscous nonadiabatic boundary layer and an inviscid isoener- 
getic shock layer. The present analysis, however, treats the entire shock layer and includes all three 
energy transport processes (radiation, conduction, and diffusion), as well as momentum transport by vis- 
cosity, and the solutions extend from the vehicle surface to the bow shock wave. The analysis is formu- 
lated in this way in order to study the coupling between the radiation and the convection. 

The effects of mass transfer from the body surface to the flow field are also considered insofar as 
they affect energy transport by radiation. In this regard, some of the simplifying assumptions of a previ- 
ous study (ref. 8) are relaxed. In that analysis, the hot air in the shock layer was only permitted to emit 
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radiation, while the cold foreign species near the body was only permitted to absorb radiation. Another 
specialization of that analysis was that both the shock layer and the foreign species in the boundary layer 
were required to be at almost a black condition insofar as radiation phenomena were concerned. That 
situation may never be realized for flight in the Earth's atmosphere, and the results may be quite different 
if it is not. In the present analysis, all of these restrictions are lifted. The air and the foreign species are 
each permitted to both emit and absorb radiation and neither is required to be black. 


ANALYSIS 


Regime of Analysis 


The flight regime in which the analysis is valid is presented in order to support some of the 
assumptions which will follow, and to relate the analysis to specific space missions. The regime is 
shown as the shaded area in figure 7-1. The upper boundary of the regime is that of equilibrium flow 
and appears as the heavy curve, taken from reference 7, and its extension to higher velocity. The 
boundary shown is for a 5-ft nose radius; for a 1-ft nose radius it would move down roughly 40,000 ft in 
altitude. The boundary on the lower right of the figure corresponds to a temperature of 15,000 K behind 
a normal shock. It is only a limit in that transport properties are presently not correlated at higher tem- 
peratures. Representative entry trajectories from references 7 and 9 (where results of ref. 10 are used in 
conjunction with ref. 9) enter the figure from the top. The ballistic coefficient is (m/CDA) and the 
lift/drag ratio is L/D. Curves for several estimated (ref. 1 1) values of the ratio of radiative to convective 
heat transfer for a 5-ft nose radius are shown to indicate the relative importance of each. Estimated 
values of the same ratio for a 1-ft nose radius would be about 10% of those shown. The heavy vertical 
curve is the approximate ionization boundary. Results will be presented from just to the left of the 
vertical ionization curve to about 50,000 ft/sec flight speed at shock-layer pressure levels of 0. 1 to 
10 atm. 


Differential Equations 

The flow model for the analysis is as before (fig. 7-2). It is also patterned after that shown in 
reference 12. The velocity components u,v at the point x,y are parallel and normal to the body sur- 
face, respectively. The coordinate x is a measure of distance along the body surface only, the distance 
of the point x,y from the surface is y, and from the body axis (for x/R « 1) 

r(x,y) = Hx ( 1 ) 

where 


H = 1 + (y/R) (2) 
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Figure 7-1.— Flight regime. 


U 



Figure 7-2.- Flow model. 

The Navier-Stokes equations for a thin shock layer (8/R « 1) for x ~ R were reduced by an order 
of magnitude analysis in reference 13 (which was completely reformulated as described in chapter 6 to 
assess specific limitations and effects). In the present analysis, these equations are reduced still further 
by an order of magnitude analysis where it is assumed that x ~ 8, restricting the results to the stagnation 
region. In addition to the transport of energy by ordinary conduction, considered in reference 13, trans- 
port by diffusion of reacting species and by gaseous radiation is included. For the radiative transfer the 
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assumptions of a grey gas, black surface, and plane-parallel shock layer are made. Detailed discussion 
of these assumptions is given in chapter 6 and in references 14-16, chapter 1. The resulting set of equa- 
uons for conservation of mass, x and y momentum expressions, and conservation of energy are 
respectively, * 


where 


dH 

9x 


8H 

dy 


_ ( p u r j ) + _(Hp V rj) = 0 

(3) 


(4) 

|=R 2+ i”< Y * 

J i=l 

(5) 

ii 

£ 

dT v 3ci I 

k w +p ^ Dijhi w) 


+ HpK 

1 2aT 4 (t)Ei(lt - tl)dt - 40T 4 + 2oT^E 2 ('t) 

Jo 

(6) 


U 2 _|_ y2 

H-h + u : v 
2 

(7) 


V 

h = Xcihi 

i=l 

(8) 

hi(T) = ( c pi dT + h? 

Jo 

(9) 


x = J pK dy 

Jo 

(10) 


and 


E n (0= dCO 

co n 


( 11 ) 
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The boundary conditions for equations (3)-(6) are at y = 0; 


at y = 5; 


where 


u — 0, v — , H — — h w 


Ux 

u = u s = — 

v = v s = -eU 
p = p s = pooU 2 (l -e) 


U 2 

H = He = 


( 12 ) 


(13) 


e 


Poo 

Ps 


(14) 


For flow near a blunt stagnation region in the absence of electric fields, it is convenient (although not 
imperative) to eliminate equation (5) entirely and replace it with 



(15) 


The advantages of using equation (15) are not only that it results in a saving of computation time 
but, primarily, it permits a simplification in the energy equation (6). That is, conduction and diffusion 
processes in the multicomponent mixture of reacting species (including effects of ionizing reactions) can 
be combined in the manner described by reference 17. Thus the transport of energy by diffusion of 
reacting species alone can be represented by a "reaction conductivity," as formulated by reference 18, 
which summed with the ordinary thermal conductivity yields an equivalent total conductivity. As a 
result, the energy equation (6) becomes 1 


3H i(J 3H 

pu '5r Mpv 


or B (p dh 1 
By (Pr By 


+ JtpK I 2aT 4 (t)Ei(lt-tl)dt-4aT 4 + 2aT^E2(t) 

Jo 


(16) 


where Pr contains the effective thermal conductivity as presented in reference 19 and chapter 5. 


'it is noted that for a binary mixture of atoms and molecules, in which diffusion and conduction effects are accounted for 
separately, the same equation (16) would result if the Lewis number were assumed to be unity (not assumed here) and 
equivalent Prandtl number were replaced by the frozen Prandtl number. 
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The integro-differential equations to be solved are now equations (3), (4), (15), and (16) subject 
to boundary conditions (12) and (13). Before proceeding with the transformation and solution to these 
equations, it is pertinent to mention the additional equations and assumptions required for studying 
effects of foreign species on radiative transport. ~ 

Foreign species effects.— Ideally, it is desirable to study the effects of a foreign species on both 
thermodynamic and transport properties in order to obtain the effect on heat transfer. However, there are 
two major difficulties associated with such a study. These are, first, that the chemical behavior of the 
mixture of air and ablated gas is usually not known, and, second, that the radiative emission properties of 
the reacting foreign species and its products are not known. For these reasons, some simplifications are 
introduced that are applicable to this portion of the analysis only. In particular, attention will be focused 
on the gross effects of the foreign species on radiative properties alone to the exclusion of its effects on 
other transport properties; that is, transport properties of the foreign gas are taken to be identical to those 
of air, except for the radiation absorption coefficient. In spite of this, a considerable amount of general- 
ity is retained as being intrinsic in the formulation of the problem. 2 The diffusion equation applicable to 
the foreign species is 


where cf is either the local mass fraction of an inert foreign species or is the pseudo mass fraction 
(refs. 20 and 21) of a reacting foreign species (which accounts for both the pure form and its local com- 
bination in reaction products). The symbol Df represents the "diffusion coefficient" of the "foreign" 
species relative to air. 


The boundary condition at the surface is obtained from a mass balance, making the assumption 
that the air does not diffuse into the surface, and is at y = 0 



v w / 1 , 

D fw (1 Cfw) 


and at the shock, the boundary condition is at y = 8 


(18) 


c f = 0 (19) 

The influence of the foreign gas on radiative transfer is obtained by its influence on the Planck 
mean mass absorption coefficient K (which will be defined subsequently), which can be expressed for 
the mixture of gases as 


K = (1 - Cf)K ai r + c f Kf 


( 20 ) 


2 For gas mixtures without ionization or light-gas elements, Lees (ref. 20) shows that it is often possible to lump all chemical 
species, including foreign reacting species, into two components insofar as diffusion transport is concerned. He shows that as 
long as the resulting Lewis number is approximately unity, heat-transfer conclusions are independent of the reaction flow 
model employed, of magnitude of reaction rates, and of any restrictions on the momentum equation. In the present work, the 
bulk of the foreign gas will be near the relatively cold wall where at least the air is not ionized and we can appropriate the 
generality of the Lees result for this portion of the analysis. 
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This is a good representation if the foreign species retains its identity (is inert) and its spectral 
characteristics are similar to those of air. It is noted that both Kaj r and Kf are functions of temperature 
and their partial pressures so that equation (20) can be written 


K - K^p^ T) 


1 + Cf 


Kf(p f ,T) 
Kair(p a ir T) 


-1 


( 21 ) 


Equation (21) indicates the need for either knowing or specifying considerable detail about the radiant 
properties of the foreign species as well as those of air. Although the latter is known, the former is not, 
so we adopt a simpler approach. We simply inquire as to the overall effect of a foreign gas that is more 
absorbing or less absorbing than air and say 

K = Kair(Ps,T)[l + Cf(a-l)] (22) 


Expressing K in this way takes advantage of the fact that K a ir(ps» T) has been correlated at several 
pressure levels in reference 22, by use of experimental data obtained by reference 23. The significance 
of a is obtained by a combination of equations (21) and (22) 


a = 


1 + 


Cf 


Kair(Pair> T) 

Kair(Ps» T) 


(1 


-Cf)+Cf 


Kf(pf, T) 
K a ir(ps> T) 


1 


For strong injection of foreign gas, 


Cf w ->1, Pf w -> Ps 


and at the wall 


a = 


Kf( P s, T) 


Kair(Ps> T) 


w 


(23) 


Thus we consider a to represent in an approximate way the ratio of the absorption coefficient of the 
foreign gas to that of air at the shock-layer pressure and wall temperature. Expressing K in the 
approximate form of equation (22) retains all significant absorption phenomena; that is, when a = 1, 

K = Kair(p s , T). When Cf is 1, K = aKaj r (p s , T). When cf is zero, K = Kai r (p s , T). Finally, when 
a = 0 (injection of a nonradiating gas), K is given by the approximation K = (1 - Cf)K a ir(Ps. T). 

Having coupled the energy equation (16) to the foreign gas diffusion equation (17) by use of 
equation (22) to obtain the influence of a foreign species on radiative transport in the gas, we proceed to 
transform all of the differential equations and their boundary conditions. 


Transformations 

Examination of equations (3), (4), (15), and (16) shows that except for the factor H appearing in 
all the equations, and the radiation term in equation (16), the equations are like the usual ones for a 
compressible laminar boundary layer. It is noted that although 3-C is essentially unity for a thin shock 
layer, its derivative with respect to y is not negligible. Thus, the appearance of H in the continuity 
equation is the most significant difference between the present equations and the boundary-layer 
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equations. For these reasons the present approach* is to modify the usual boundary-layer transforma- 
vSles"^ 6 ^ effeCt ° f Thetransformsusedtochan g e x and y to % and T] as independent 


= 1 Ps4 fe 


u s dx 


A stream function is defined such that 

|^ =purJ > |^ = -«pvri 

which satisfies the continuity equation (3). Other definitions are 

V2jf 


Then from equations (25)-(27) 




— T 
T = ^- 
T s 


PsHsH^ 


u _ df_ 
u s dr\ 


Applying the transforms and definitions to equations (4) and (16), noting that (u 2 + v 2 )/2 « H in the 
stagnation region, yields 


3 An alternate approach in which the continuity equation is reduced further is presented in reference 12, pp. 


pp. 390-393. 
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L*| +f £L. 2t (i?a 

rail 2 / an 2 i Us 


|L) - 2 2=(l-e) 
°nl p 


t df_ 9 2 f df\ 
^1 dr\ 2 d$l 


( 33 ) 


and 


_8_U9g\ + f 9£ _ I 
dri \Pr dri / dr\ 


4RaT^K s 


(j + DU 3 


K 


f 

Jo 


T^OEiflt - xl)dt - 2T 4 + T^x) 


2 p(df dg df 9g| 

s ia^ aria^j 

Finally, in order to reduce these equations further, it is assumed that 

Psji s = constant 


(34) 


(35) 


and that similarity exists; that is, all dependent variables are functions of r\ alone. Equations (33) and 
(34) become 


(cpf')’ + ff" = - 


1 


(j + DL 


(f) 2 _ 2^(1 -e) 
P 


(36) 


and 


£g'l+f«' = - 


= Q 


4RcT^K s 

K 

f -rifOEiflt - xi)dt - 2T 4 + T^E 2 (x) 

a + i)u 3 


Jo 


(37) 4 


In order to relate x to q conveniently, it is assumed that p = p(y) and thus from equations (10) and (25) 

f- 

x = p I K dp (38) 

Jo 

where 

P = K s / JPs^sR (39) 

V U(i + 1) 


The boundary conditions, equations (12) and (13), become at TJ = 0 


4 For the case in which the flow field is emitting radiation, but is transparent to radiation, the right side of equation (37) is 
expressed differently. That case is treated briefly subsequently. Except where stated to the contrary, the integral expression 
for radiative transport as shown in equation (37) was used throughout the analysis. 
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at r| = T| s , 


f — f w Pw v w 



R 

PsM- s U(j + 1) 


f = 0 


g = gw 


2h w 

U 2 


( 40 ) 


/ RU 

f — f S — Poo y\f — — 

V psPsO + 1) 

f = f's= 1 t 

g = gs = 1 


(41) 


It is noted that q s is actually determined where the above equations for f s and f s are simultaneously 
satisfied. The diffusion equation (17) is transformed in like manner and becomes 



+ fc'f = 0 


Its boundary conditions are from equations (18) and (19) at 11=0 


(42) 


c’f = c'f w =(^-) fw(l-Cf w ) (43) 

19 / w 

where Sc is the Schmidt number |i/pD. At q = q s 

cf = c fs =0 (44) 

The equations to be solved are (36) and (37) (and (42) if foreign species effects are considered), 
subject to boundary conditions (40) and (41) (as well as (43) and (44)). Heat transfer to the surface and 
shock standoff distance are derived from the solutions of the equations. They are discussed in the next 
section. 


Standoff Distance and Heat Transfer 

The shock standoff distance, obtained from equations (25), (26), and (41), is simply 
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Subsequently, it will be compared with the approximation given by Hayes (ref. 24) which is 

8 Re 
~ 1 + f2z 

The total heat-transfer rate at the stagnation surface is 


qj - qc + qr 

where the convective heat-transfer rate is 


(45) 


(46) 


(47) 



and the net radiative heat transfer at the surface is (for a black surface) 


q r = aTjTt-2| ^(t^todt 


r 

Jo 


(48) 


(49) 


Equations (47)-(49) can be combined and transformed by equation (25). The result is the total heat- 
transfer expression 


q T — 


|_ 9 j U 5/2 / PsM-sO + 1) 

\Pr/w 2 V R 


g'(0 )+ oT^tJ,- 2 f T 4 (t)E 2 (t)dt 

Jo 


(50) 


The radiation flux from the hot gas cap to the wall, qg, is represented by the integral in equations 
(49) and (50). For the case of a transparent isothermal gas cap, this integral, with the aid of equation 
(10), and the fact that in this instance E 2 (t) * 1 and T « 1 becomes 

q g = -2aTj f T 4 E 2 (t)p s K s dy = -2aT^p s K s 5 (51) 

Jo 

The total energy emission rate per unit volume of gas is shown by reference 16 (p. 31) to be 

E t = 4jtp I B v K v dv (52) 5 

Jo 


5 It is useful to note that E t for equilibrium air can be obtained from reference 23 by the relationship E t = 2oT 4 (e/L) where 
e/L (the emissivity per unit depth of gas) is the terminology of the reference and not of this chapter. 
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The definition of the Planck mean mass absorption coefficient is 


Ks 



( 53 ) 


From equations (51) and (52) the transparent isoenergetic gas-cap radiation is 


q g =- 



(54) 


This expression can be compared with that used in reference 1 1 for the plane-parallel isoenergetic 
transparent gas cap represented in figure 7-3. The elemental volume dV is emitting energy at the rate 
of E t per unit volume. The rate at which energy is received by the wall area dA from the gas volume 
dV is 


dqgdA = - 


E t dV 


4k 


I 


dA cos 0 


where the first set of parentheses is the energy emission rate per unit solid angle from dV, and the sec- 
ond set is the solid angle subtended by the wall area dA. Expressing the volume element dV in spheri- 
cal coordinates with the origin at dA, and integrating overall gas space yields the net radiant heating rate 
per unit wall area 6 


6 If, instead of the plane-parallel assumption, it is assumed that the body and shock are concentric hemispheres, equation (52) 
becomes 





1/2 Y V 2 a I 

+ 4 — 4- + 0[y 3 ] 
3 8 48 


where y = [(8/R) + 2](8/R). For S/R « 1, the expression is 



1 


2 / 2 / 8 | 1/2 " 
3 \R / 


which reduces to the plane parallel case for 8/R -> 0. Thus it can be seen that for the case of 8/R = 0.04, equation (54) 
overestimates the radiant flux by approximately 15%. From this we see that the assumption of a plane-parallel gas cap is 
good for the transparent case, and is expected to be belter yet for reabsorption. 
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which is the same as equation (54). 


Qg-- 


E t 5 


2 


SHOCK 



Figure 7-3 - Isoenergetic transparent gas cap. 

Having formulated the problem and expressed quantities of interest to be derived from its solu- 
tion, it is instructive to mention some features of the method of solution. 


Solution of the Equations 

In this section we will first discuss some of the information needed to initiate the solutions. Then 
the sequence of steps used in the numerical integration will be described in a general way. Additional 
information pertaining to the evaluation of the radiative flux term in the energy equation will be pre- 
sented in an appendix. 

The differential equations (36) and (37)_subject to boundary conditions (40) and (41) can be inte- 
grated if the transport properties <p, cp/Pr, andK and the thermodynamic properties pandT are known 
as functions of enthalpy at constant pressure. Except where noted, these are obtained from reference 22. 
In that reference, the transport and thermodynamic properties of Hahsen (ref. 19) and Moeckel and 
Weston (ref. 25) for equilibrium air to 15,000 K were represented by correlation formulas. The Planck 
mean mass absorption coefficient for equilibrium air obtained from reference 23 was also expressed by 
correlation formulas in reference 22. In passing, it is noted that since H ~ 1 everywhere, <p = p[i/p s |i s . 
However, the derivative of cp is also used and must contain the derivative of Jt; thus, 



Obviously, the derivative of cp/Pr must also contain the derivative of H. 

The functions Ei(£) and E2(C) are also needed in the solution of equation (37) and in the radiant 
heat-transfer equations (49) and (50). These can be obtained from the tables of reference 26 in fairly 
coarse intervals of £ or (as in the present work) can be generated from the exponential integral tables of 
reference 27 in very fine intervals noting that 

Ei(C) = -Ei(-0 (56) 
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where the E, designates the exponential integral and is not to be confused with Ei . . .E n defined by 
equation (11). The function E 2 (£) and subsequently E 3 (Q) is generated from equation (56) by the 
recurrence formula for n > 1 


E^O-jrrte-C-CEn-l©] 


(57) 


The procedure for solving the differential equations begins with the assumption of profiles of <p, 
<p/Pr, p- 1 , and Q as functions of T|, where Q is the entire right side of equation (37). Equation (36) is 
then integrated numerically using the Adams-Moulton predictor-corrector variable step integration 
scheme (ref. 28). The solution is obtained by iteration to find the value of f '(0) that makes the resulting 
solution satisfy the first two boundary conditions (41). The solution of equation (36) is used to solve 
equation (37) directly without iteration by evaluating the integrated form of equation (37) which is 
obtained as follows. The energy equation (37) can be rewritten: 





(58) 


Let 



(59) 


Equation (58) becomes 


drj (p 


(60) 


An integrating factor 0 is defined by 

i 


f(Pr/<p)ds f(Pr/cp)ds 

0=e“ , 0' = e 0 • f — = 0 

<P 


f*| 

9 


( 61 ) 


Multiply equation (60) by 0, 


0 •— +af— -0 
dri q> 


0-Q(q) = O 


(62) 


Then 
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and 


^- = 0Q(ri) 


(63) 


a0 = 


I 


0Q(ri)dri + const. 



n 

f(Pr/cp)ds 


ms f 

Q0l)e ( 


f(Pr/<p)dw 


ds + const. 


At T| = 0, g' = g' w yet unknown. Say 


const. = 


£ 

Pr 


± g 


w 


w 


Substitute equation (66) into (65) to obtain 


g’(Tl) = 


dg = 
dq 


J 

r 

f(Pr/tp)ds 


<p/Pr • e Jo 



We integrate again with, at r| = 0, g = g w : 


(64) 


(65) 


■ (66) 


(67) 
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gw + 


gCn) = 


T 

Pr 


w 


gw 


J. 


(cp/Pr) • e 


•s 

Jo 


f(<p/Pr)dco 



which is the solution of the transformed energy equation. To evaluate g' w , let q =q s in equation (68) 
and use the boundary condition that g(T| s ) = 1 and solve for g' w , which is 


g w =- 




(cp/Pr) 


w 


dq 


(<p/Pr)e‘ 


•n 

f 

Jo 


f(Pr/<p)ds 




X 


< 1 - gw - 


1 


r o j f(Pr/<p)ds 

((p/Pr)e 0 


r 

Qe 


f(Pr/<|>)dco 


dq 


(69) 


If a foreign species is considered, equation (43) is used to integrate equation (42) with the result 


cfCH ) = c fw + f w (l - c fw ) I — e Jo 


- (Sc/cp)fdco 


ds 


(70) 



which is evaluated as before without iteration so that cf w is selected to make cf(T) s ) equal to zero in 
accord with equation (44), thus 


C fw=i 




1-1/ 


*w 


Sc 

<P 


. Jo 


(Sc/q>)f ds 


doo 


( 71 ) 


_ The resulting profiles of g(n), g'Ol). and Cf(r|) are used to calculate new profiles of cp, <p/Pr, 
p _1 ,T, and K. The corresponding profiles of Q are calculated as described in the appendix. Using the 
new profiles, the next major iteration is begun by seeking a new solution of equation (36) and so on. 
The procedure is repeated until f '(0) does not change from one major iteration to the next, at which 
point 5, qr, and q^ are calculated from the solution of the differential equations using equations (45), 
(49), and (50). The method of evaluating equation (49) is presented in the appendix. 


RESULTS OF THE ANALYSIS 


The solutions of the integro-differential equations, including the resulting fluid flow and property 
profiles across the shock layer, are presented for a variety of flight conditions. Then heat-transfer results 
and shock standoff distance are shown and discussed. Finally, results with mass transfer are presented 
and their significance with regard to heat transfer is discussed. 


Solutions 

Nonionized regime- In presenting solutions of the differential equations, results for flight in the 
nonionized regime (but near the ionization boundary in fig. 7-1) are examined first in order to evaluate 
the present analysis in the light of known phenomena. It is pointed out that for the low-speed case pre- 
sented (fig. 7-4 only) two specializations have been introduced in the differential equations. These are 
(1) that the fluid in the gas cap emits, but is transparent to radiation, and (2) conduction and diffusion 
effectors are separated (and thus the frozen Prandtl number is assumed to be 0.72, and the Lewis number 
is unity in keeping with footnote 1). Both specializations are valid where they are used. The conse- 
quence of the first one is that the integral radiation term in the energy equation (37) is replaced by 


2R 

Et 

(j + 1)U 3 _ 

P 


( 72 ) 


where E t (p,T) is obtained from reference 23. 
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Figure 7-4- Flow-field solutions for U = 32,000 ft/sec at 1 10,000 ft altitude; p s = 10 atm, R = 5 ft. 

A solution is shown in figure 7-4 for flight at 32,000 ft/sec at an altitude of 1 10,000 ft 
(p s = 10 atm) for a body having a 5-ft nose radius. It was obtained by setting <p = 1 in the flow equa- 
tions. When cp was allowed to vary in accord with the computation of p|i by the correlation formulas 
of reference 17 for dissociated nonionized air, the resulting profiles could hardly be distinguished from 
those shown for 9=1. The standoff distance was the same in either case (5/R = 0.0409) and was about 
3% below the value estimated by Hayes' method (ref. 24). Although the radiative heat transfer was 
unaffected by the assumption 9 = 1» the convective heat transfer (which is presented subsequently) was 
approximately 15% lower for the 9=1 case, an effect that has been observed and discussed in refer- 
ences 20 and 29, and was as expected. 

The structure of the flow field indicated by the f and f profiles is that of a momentum boundary 
layer with constant external vorticity. The thermal structure indicated by the g profile is somewhat like 
that of a thickened thermal boundary layer — much like that of reference 30. 

The above comments apply generally to solutions obtained in the nonionized regime. At lower 
shock-layer pressure levels, the standoff distance was still closer to Hayes' estimate and the flow-field 
structure was even more like that of a boundary layer joined to an isoenergetic shock layer. 

Thus, in every respect, the present results, for which the flow field is not separated into a bound- 
ary layer and inviscid flow field, agree with results obtained when the flow field was so divided. The 
use of equation (15) instead of equation (5) and, further, the agreement of the present results with those 
obtained previously justify the use of equation (15) instead of equation (5). We thus conclude that the y 
momentum equation (5) is not important in the stagnation region flow field, for the high Reynolds num- 
ber regime. Having established the validity of the present method, we now apply it at higher flight 
speeds. 


Ionized regime- In the ionized regime, results are presented with and without radiative transfer 
to show its effect on the solution. Properties, computed by the use of reference 22, include the effects of 
ionization. 
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Flow-field profiles corresponding to flight at 50,000 ft/sec at 190,000 ft altitude (p s = 1 atm) for 
R = 1 ft without radiative transport effects are shown in figure 7-5. Here the boundary-layer structure is 
less distinct than was exhibited at lower flight speed. Far from the wall, however (y/5 > 0.3), the flow 
field is isoenergetic and has almost constant vorticity as indicated by the nearly constant slope of the 
u/u s curve. The transport property parameter ( <p/Pr )/(cp/Pr) w is shown in the figure. Using the criterion 
of reference 22, it is noted that the region between the body and the peak in this parameter (essentially 
3% of the standoff distance) is essentially not ionized. 



Figure 7-5.- Flow field with ionization but without radiation effects; U = 50,000 ft/sec at 190,000 ft 

altitude, R = 1 ft, p s = 1 atm. 

A solution at the same flight condition, but with radiative transport expressed by the integral 
form included in the energy equation (37), is shown in figure 7-6. The enthalpy level is substantially 
diminished from that of the previous Figure because of the transport of energy by radiation. Here the 
thermal structure of the flow field is that of a very thick thermal boundary layer joined to a nonisoener- 
getic shock layer, and the question is raised as to the applicability of first-order boundary-layer theory at 
this flight condition. Some speculation on this question will be made subsequently. 

The effect of increasing nose radius can be seen by comparing figure 7-7 with figure 7-6. In fig- 
ure 7-7, the flow field is far from isoenergetic everywhere. This behavior will, of course, affect the 
convective heat transfer as will be shown subsequently. 
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Figure 7-6- Flow field with both ionization and radiation effects; U = 50,000 ft/sec at 190,000 ft 

altitude, R = 1 ft, p s = 1 atm. 



Figure 7-7.- Flow field for large nose radius; U = 50,000 ft/sec at 190,000 ft altitude, R = 5 ft, 

p s = 1 atm. 

Solutions at higher and lower pressure levels (10.0 and 0.1 atm) are shown in figures 7-8 and 7-9. 
Boundary-layer structure is, of course, more evident in the high-pressure profiles (fig. 7-8). 
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Figure 7-8- Flow field profiles for U = 40,000 ft/sec at 1 18,000 ft (p s = 10 atm, R = 1 ft). 



Figure 7-9 - Flow field profiles for U = 50,000 ft/sec at 248,000 ft; p s = 0.1 atm, R = 5 ft. 

A solution with mass addition at the surface is shown in figure 7-10. The flight condition is the 
same as that for figure 7-6. A comparison of the two figures shows that the profiles for the injection case 
are shifted a small amount to the right for 0 < y/8 < 0.5, but coincide with the no-injection results for 
y/8 > 0.5. 
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Figure 7-10.- Flow field profiles with air injection at the wall; U = 50,000 ft/sec at 190,000 ft altitude, 

R = 1 ft, p s = 1 atm, f w = -0.4. 

Although figure 7- 10 is for air injection at the wall, the profiles (except for c f ) resemble those for 
a foreign gas, with a = 50, injected at the same mass-flow rate. The latter result is shown in figure 7-11. 
In general, the flow field is slightly cooler than that for air injection. 



Figure 7-11.- Flow-field profiles with foreign gas injection at the wall; U = 50,000 ft/sec at 190,000 ft 

altitude, R = 1 ft, p s = 1 atm, f w = -0.4, a = 50. 

A solution was attempted for a larger injection rate (f w = -0.7), but it became unstable. Although 
the solution did not diverge indefinitely, the profiles oscillated, apparently within limits, with sufficient 
amplitude to render them invalid. Of course, for large injection rates (corresponding, for example, to 
high ablation rates caused by large heating rates), the flow field will become physically unstable. 
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Heat Transfer 


Attention is now directed to the stagnation-region heat-transfer rates obtained from the solutions. 
Convective heat-transfer rates are presented and compared with those of others. The effect of radiative 
transport on convective heat transfer is shown. The conflicting theoretical (and experimental) results of 
various investigators for convective heat transfer mentioned earlier are discussed and an explanation is 
advanced for the differences. Radiant heat-transfer results are presented and compared with those cal- 
culated by use of the isoenergetic shock-layer assumption. 

Convection- The results of convective heating analysis are presented in the form of reference 6 
in figure 7- 12(a). Results of references 2, 3, and 6 are shown for comparison. It is seen in figure 7- 12(a) 
that if radiative transport is neglected, the present result agrees well with the results of references 2 and 
3. However, it (like the results of refs. 2, 3, 4, and 5) differs greatly from the results of reference 6 at 
high speed. The result of reference 29 is also shown for comparison, although ionization effects were 
not included in that analysis. 

The computed results of Howe and Viegas end of Scala are compared with experiments as shown 
in figure 7- 12(b). The experiments of Warren (ref. 31) agree with the theory of Scala, while the experi- 
ments of references 2 and 31-34 are more in agreement with Howe and Viegas. Although the experi- 
mental discrepancy was never resolved (to my knowledge), the theoretical results were, as will be dis- 
cussed subsequently. 

When radiative transport effects are included in the differential equations, the convective heat 
transfer is diminished, as shown in figure 7- 12(c), the effect being larger for larger nose radii and larger 
pressure levels. For R = 5 ft, p s = 1 atm, and U = 50,000 ft/sec, the convective heating is less than 50% 
of what it would be if radiative transport were ignored 7 and is only 16% of the value predicted by the 
theory of reference 6. 

Insofar as convective heat transfer is concerned, two questions are suggested. First, what is 
responsible for the difference between the results of reference 6 and those of other analyses when radia- 
tive transport is ignored? Second, is it possible to use the modified boundary-layer edge properties with 
existing boundary- layer results to predict the effects of radiative coupling on convective heat transfer? 

Of course, some of the differences between the results of reference 6 and those of other analyses 
can be attributed to the assumed chemical state of the gas as well as to the method of formulating and 
solving the problem (comparison of several analyses is made in ref. 35). However, it is likely that the 
principal reason that reference 6 obtains much higher convective heating results than do other analyses 
lies in the transport properties. 

Significantly, references 2-4 and the present analysis all use the transport properties of Hansen 
(ref. 19), but reference 6 does not. A comparison of the thermodynamic and transport properties of 
reference 6 with those of reference 19 yields the following result. The properties of viscosity, frozen 
specific heat, frozen thermal conductivity, and frozen Prandtl number generally agree within a few 


7 It is noted that although the dashed line labeled "present result without radiation coupling" in figure 7- 12(b) is shown to be 
for p s = 1 atm and R = 1 ft, it is not very sensitive to changes in pressure level and nose radius and can be used as a basis 
for comparison for all R and p s . 
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Figure 7-12.- Comparison of convective heat-transfer rates, (a) Radiative transport neglected in flow 
equations, (b) Convective heating — comparison with experiment, (c) Effects of radiative transport 
in flow equations. 
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percent except at high temperatures, where viscosity differs by approximately a factor of 2, which also 
shows up in the frozen Prandtl number. 8 However, when the binary diffusion coefficients of reference 6 
are put into the form of an equivalent reaction conductivity and added to the thermal conductivity, the 
result is very different from Hansen's total conductivity at high temperatures as shown in figure 7-13. At 
5000 K, the value deduced from reference 6 agrees well with Hansen’s value. At 10,000 K, the former is 
greater than the latter by a factor of about 25; and at 13,000 K, by a factor of about 20. 


REF. 19, HANSEN, AIR 



Figure 7-13.- Total thermal conductivity comparison; p = 1 atm. 

Thus, it is likely that the difference between the convective heat-transfer results of reference 6 
and those using Hansen's properties can be attributed to the reaction conductivity (or to the equivalent 
diffusion coefficients or their corresponding cross sections). Some experimental values of total thermal 
conductivity in nitrogen obtained from reference 37 are shown in figure 7-13. They are in substantial 
agreement with both the curve obtained from reference 38 for N 2 and Hansen's curve (ref. 19) for air, but 
do not agree with the values deduced from reference 6. However, more experimental work on conduc- 
tivity, diffusion coefficients, or cross sections is required before the problem can be completely settled. 


Turning to the second question, it is noted that convective heat transfer has been correlated by 
means of the driving enthalpy across the boundary layer and other outer edge boundary-layer properties 
(refs. 2, 3, and 29). The question is, can the reduction in convective heat transfer due to radiative cou- 
pling be accounted for by simply making the driving enthalpy (He - H w ) less than (H s - H w ) by an 
amount corresponding to the radiant energy emitted from an isoenergetic shock layer. The result is not 
encouraging; even for a small nose radius where coupling is expected to be relatively small, the 


According to a private communication from Scala, ihe viscosity in reference 6 is off by a factor of 2 because of an error in 
reference 36, but when viscosity is corrected, there is no appreciable change in the convective heat-transfer result of 
reference 6. 
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approximation gives an erroneous result. For example, the driving enthalpy (and thus the estimated 
convective heat transfer) for a body with a 1-ft nose radius flying at 50,000 ft/sec at 190,000 ft altitude 
would be diminished by about 8% because of this radiant emission, but the actual convective heat trans- 
fer calculated by the method of this paper is diminished by about 15% from the case where radiative 
transport is neglected (as seen in figure 7- 12(b)). Furthermore, the comparison is still worse if the 
change in boundary-layer edge properties corresponding to the diminished driving enthalpy is taken into 
account. Then the estimated convective heat transfer would increase by 5% in contrast to the actual 
decrease of 15% because p e |ie increases as He decreases. 9 The situation is much worse for larger nose 
radii. 


It is speculated that this failure of simple boundary- layer theory may be attributed to failure to 
match gradients of the flow properties at the outer edge of the boundary layer. 

Radiation- Studies of radiant heat transfer from the isoenergetic transparent shock layer have 
been made in references 7 and 1 1. Our purpose here is to indicate how the radiant heat transfer may be 
altered if transport processes are considered, and if the transparent gas assumption of those references is 
relaxed. 

The effect can be seen by comparing the solid line with the dashed line in figure 7-14. It is seen 
that for larger nose radii at 50,000-ft/sec flight speed with p s = 1 atm, the radiant heat transfer is actually 
about half what would be predicted for an isoenergetic transparent shock layer. It is noted that a similar 
effect was obtained for the convective heat transfer at this flight condition. Thus the coupling of trans- 
port processes reduces the convective heat transfer, and the fact that the shock layer is not isoenergetic 
reduces radiative heat transfer by large amounts at severe flight conditions. 


2.0 x io 4 



Figure 7-14.- Effect of nose radius on radiative heat transfer; U = 50,000 ft/sec at 190,000 ft altitude, 

p s = 1 atm. 


9 It is assumed for illustrative purposes that (after Cohen, ref. 3) q c <* (H e - H w )(p c p < .) (1 - 4 ^. 



The ratio of actual radiative heat transfer to that which would be predicted by an isoenergetic 
transparent shock layer assumption 10 is shown in figure 7-15 for several flight conditions. Interestingly 
enough, the ratio drops off more rapidly at 40,000 ft/sec than it does at 50,000 ft/sec (for p s = 1 atm). It 
can be shown by a detailed comparison of the two solutions for R = 1 ft (which is beyond the scope of 
this chapter) that the result can be attributed to real-gas properties. 



R, ft 


Figure 7-15.- Comparison of radiative heat- transfer results to those of isoenergetic shock layer 

(T w = 1500 K). 


Heat-transfer correlation- Attempts to correlate the convective heat-transfer results shown in 
figure 7- 12(c) in a simple way have not been successful. However, net radiative heat transfer (radiation 
flux incident on the wall less that emitted by the wall) has been correlated by the expression 


— -qr VR/ ggL_ - -0.0024 + 0.053 p s K s R — 

PooU(H s -H w )(U/U re 0 5/4 Ps 


(73) 


where R re f, U re f, and p re f are 1 ft, 10 4 ft/sec, and 21 17 lb/ft 2 , respectively. The correlation is shown in 
figure 7-16, and includes results for flight speeds between 30,000 and 50,000 ft/sec, nose radii between 
1 and 5 ft, and shock-layer pressure levels between 0.1 and 10.0 atm, at wall temperatures of 1000 to 
1500 K. 


10 Both q r and <f are the radiant heating rate to the wall less that emitted by the wall. Equation (49) is used to compute q r 


while equations (49) and (51) are used to obtain q* for an isoenergetic transparent flow field. The result is 

q; = <nf(Tt-25p s K s ). 
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Figure 7-16.- Correlation of radiative heat transfer. 

If it is noted that p s = poJT 2 and (H s - H w ) » U 2 /2, the expression for net radiant stagnation- 
region heating rate becomes 


-qr = pooU 


1/41 


u 4 


Cli = + C2^VR'lJ 2 

VR e 


(74) 


where Ci = -0.12x10-7 sec 5 / 4 fr 3 / 4 , C2 is 0.561xl0" 3 sec 5 / 4 lb fr 11 / 4 , and the units for q r are 
ft lb/ft 2 sec. 

Effect of mass addition on heat transfer- When air is injected from the surface into the flow 
field (corresponding to the solution shown in fig. 7-10), the convective heat transfer is only about 65% of 
that when there is no injection (solution shown in fig. 7-7), but the radiative heat transfer is about 106% 
of the no-injection value. The increase in radiation is caused by a 5% increase in shock standoff distance 
because of the mass addition into the flow field. 

If a foreign gas that is approximately 50 times as strong an absorber and emitter of radiation as 
air is injected into the flow field, the convective heat-transfer result for the foreign gas injection is about 
86% of that for air injection. The decrease in convective heat transfer is accompanied by an increase in 
radiative heat transfer of about 7%. Standoff distance for the foreign injection is approximately 10% less 
than that for air injection, possibly because the increased radiant emission for the former cools the shock 
layer and increases the density slightly. 

It is interesting to note that injection of a gas that is 50 times as strong an emitter as air increases 
the radiant heat transfer by only 14% over the no-injection case for this flight condition. The reason is, 
of course, that most of the mass of foreign gas is in the cold part of the flow field, as can be seen from 
the profiles of p/p w , Cf, andT in figure 7-11. 


• i 1 1 
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Shock Standoff Distance 


In the absence of mass transfer at the body surface, there appears to be very little effect of trans- 
port phenomena on shock standoff distance at flight speeds up to 50,000 ft/sec for the cases considered. 
For example, at a speed of 40,000 ft/sec for p s = 1 atm, 5/R varies from 0.0467 to 0.0468 for all nose 
radii, with or without radiative transport in the flow equations. These values are very near the 0.0472 
predicted by Hayes' method (ref. 24). 

The standoff distance for a flight speed of 50,000 ft/sec without mass addition is shown in fig- 
ure 7-17. The open symbol is the present result without radiative transport considered and is about 7% 
above Hayes' estimates. With radiative transport, the present result lies from 3 % about to 1 1 % below the 
estimate of Hayes. 



R, ft 

Figure 7-17 - Shock standoff distance; U = 50,000 ft/sec, p s = 1 atm. 

It is expected that energy transport will diminish the standoff distance because it tends to cool the 
flow field and raise the density. The present results indicate that this is the case. It appears in fig- 
ure 7-17 that as nose radius increases above 4 ft, and radiant transport becomes more severe, the standoff 
distance tends to grow less rapidly with R. Although results are not available, it may be that there is a 
strong transport phenomena effect on standoff distance for very large radii (R > 5 ft) at severe flight 
conditions. 

Finally, for the cases of moderate rates of air injection 11 at the surface, the standoff distance 
increased approximately 5 % above the no-injection case, while for a strongly absorbing gas (a = 50) 
injected at the same rate, the standoff distance was about 6% below the no-injection value. 


11 A moderate injection rate is one low enough not to upset the stability of the flow field. The limiting rate corresponds 
roughly to that which would just "blow off' the laminar boundary layer (refs. 8 and 39). 
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SUMMARY 


In this chapter the derivation has been used for solving the entire flow field in the stagnation 
region of bluff bodies where the flow field cannot readily be separated into a shock layer and a boundary 
layer. The method has been tested wherever possible to establish its validity. At flight speeds of about 
30,000 ft/sec it has reproduced known results in terms of heat transfer, standoff distance, viscous effects, 
and flow-field structure. At higher speeds (in the ionized flow regime), it has reproduced known heat- 
transfer and standoff results when radiative transport effects were neglected. The method has been 
applied to study the effects of coupling between radiative and convective transport as well as the effects 
of mass addition on energy transport. A number of general conclusions may be drawn. For high-speed 
flight in the equilibrium regime (flight speeds of the order of escape speed and higher for shock-layer 
pressure levels of about 1 atm and higher) strong coupling among the various modes of energy transport 
exists. The prediction of convective heat transfer without including radiative coupling effects overesti- 
mates convective heating by a factor of about 2 for flight at 50,000 ft/sec for R = 5 ft. Similarly, the 
prediction of radiative heat transfer by the isoenergetic approximation overestimates radiative heating by 
a factor of about 2 for the same flight condition. 

It is interesting to note that coupling between convective and radiative heat transfer generally 
reduces the convective heating rate more than can be accounted for by considering the fact that the 
driving enthalpy for convective heating has been diminished by the emission of radiant energy from the 
shock layer. The conclusion is that the combination of all modes of transport in the entire flow field is 
important in determining convective heat transfer at the wall. 

Without mass addition at the body surface, only about 7% to 3% of the flow field is free from 
ionization at flight speeds from 40,000 ft/sec to 50,000 ft/sec for p s = 1 atm. This, of course, is the 
region near the cold wall. Its extent appears to be about doubled when a gas that does not ionize more 
readily than air is injected (at moderate rates). 

Injection of a foreign gas that is 50 times as strong an emitter of radiant energy as air increases 
the radiant heat transfer by only 14% while diminishing the convective heat transfer by about 45% com- 
pared with the no-injection case (for U = 50,000 ft/sec, p s = 1 atm, R = 1 ft). The total heat transfer was 
diminished by about 14%. The radiant heating increase is small because the bulk of the foreign gas 
remains near the cold wall. Very little foreign gas gets into that part of the flow field where the air is 
ionized. 

Finally, these results of real-gas effects were obtained by use of what the author believes are cur- 
rently the most reasonable values of real-gas properties. The results are considered to be quantitative to 
the extent that the gas properties are quantitative. In any event, the phenomenological effects shown in 
the results are considered to be qualitatively correct. 
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APPENDIX 


EVALUATION OF RADIATION INTEGRALS 


The term Q in the energy equation (37) that accounts for the emission and absorption of radiant 
energy is 


Q = 


Jr r 

(j + DU 3 Jo 


T^COEKIt - xl)dt - 2T 4 + T^x) 


j 


(Al) 


The expression must be evaluated at every point in the flow field between the body and the shock. 

The terms in the brackets can be expressed as a finite series if the curve T^t) is replaced by M 
straight-line segments in equal intervals (At = x^/M) of x and the result is integrated. The brackets then 
become 


f 

Jo 


T 4 (t)Ei(lt - xl)dt - ZI 4 + T^E 2 (x) 


M-l 

= - ( E 2 (t s - X) + j- X CT?+1 - TV)[E3(lti+l - Xjl) - E 3 (lti - Tjl)] 
At i=0 


(A2) 


The series is evaluated at a constant x (starting with x = 0), with the dummy variable t ranging 
from t = 0 to t s . This leads to one value of Q corresponding to that Xj (or to T], which is related to x 
by eq. (38)). Then x; is increased by Ax and the entire process is repeated until x = x s . The method 
can be illustrated in detail by use of table Al. 

At the end of each integration of the^ energy equation (37), new profiles of the various thermody- 
namic and transport properties (including T and K) and the optical depth x are obtained as functions of 
T] from the new profiles of g. The resulting optical depth at the shock x s = t s (= 0. 10 for purposes of 
illustration in the table) is divided into M (= 10 in the table) equal intervals in x^anjJ t sojjiat 
xo = O...Xj = jA...XM = x s = MA and to = 0. . .ti = iAt. . .tM = t s = MAt. A table oF, T i+1 -Tj and, for 
each value of t, a table of E 3 (lti - xl) corresponding to ti is formed. Now it is easy to see the appli- 
cation of the table for evaluating the right side of equation (A2). The term E 2 (t s - Xj) is evaluated by 
use of the table and the recurrence formula (57). The series is evaluated at each Xj (for example, 0.02) 
using column 3 and the column corresponding to Xj (column 7). 

Because of the symmetry of the table E 3 (ltj - xl) about its diagonal, the table can be reduced to 
only two columns of E 3 values. It is noted that the column corresponding to Xj differs from that of 
Xj_i only in that the former is shifted downward one step and the top value in the Xj column is the ith 
value in the Xj = 0 column. Thus, using the Xj = 0 column as a master column plus one working Xj 
column makes it possible to generate the next working column. In that way, only 2(M +1) values of E 3 
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need be retained in the digital computer memory instead of (M + 1)2 values. A very large saving in 
memory space is achieved for large values of M. 


In similar fashion, the integral in the heat-transfer equation (50) can be replaced by a finite series 



'I 4 (t)E 2 (t)dt-Tt = -2 



(A3) 


Here, t — 0 and the series can be evaluated by a one column table of F.^(q), 
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CHAPTER 8 


HYPERVELOCITY FLOW FIELDS WHICH HIGHLIGHT THE EFFECTS OF MASS 

ADDITION 

INTRODUCTION 


In this chapter, many solutions of the viscous shock layer equations with mass addition are 
obtained. As before, flow-field equations include the effects of heat conduction, diffusion of reacting 
species, and emission and absorption of gaseous radiation for dissociated and partially ionized air in 
chemical equilibrium. Convective and radiative heating rates with mass addition are obtained from the 
solutions. Algebraic equations are derived for predicting the nose radius that minimizes total heating 
rates at a given flight speed and shock-layer pressure level. Values for the corresponding natural abla- 
tion rate, both the intrinsic and the effective heat of ablation, the ratio of radiative to convective heating 
rate, surface shear stress, and shock-wave standoff distance are given. The effects of ablated gases that 
radiate more strongly than air are examined. Rules for scaling flow-field structure with mass addition 
are discussed. Solutions without mass addition at low Reynolds numbers where external vorticity, 
energy depletion, and flow energy limiting are important are compared with existing theory and experi- 
ments. This chapter is largely the text of a report originally published as reference 1. 

The study of mass addition in the stagnation region is of interest because any object, blunt or 
pointed, which enters the atmosphere at high speed will generally have a blunted stagnation region as the 
material suffers thermal erosion. Interest is further enhanced because the aerodynamic heating rate is 
likely to be a maximum in the stagnation region. 

The hot, thin gas cap over the forward surface of an object entering a planetary atmosphere is the 
host to a myriad of interrelated physical phenomena. The study of the gas cap is especially cumbersome 
if the flow field is a mixture of air and foreign species which were added to the flow because the surface 
is ablating. The knowledge of mass addition effects at speeds below which ionization and gaseous radi- 
ation effects may be neglected is highly developed from both the flow-field and materials points of view. 
(A small part of the extensive literature on the subject will be brought into the discussion subsequently 
where appropriate.) On the other hand, mass addition at speeds greater than 30,000 ft/sec for which the 
gas cap is both ionized and radiating has received comparatively little attention. 

The present purpose is to examine the effects of mass addition at flight speeds greater than 
30,000 ft/sec as it influences and is influenced by some of the many other phenomena, parameters, and 
physical properties of the gas cap. To this end we consider mass addition in the general sense, transpi- 
ration; and in the special sense, ablation. 

For mass addition in general, we particularly want to know its influence on convective heating. 

Is heat blockage as effective at the higher speeds as it is at the lower speeds? Do the existing correlation 
formulas obtained for lower speeds apply at hypervelocity? 

With respect to ablation, we are especially interested in finding the conditions for which the total 
heating rate at a given flight condition and given material is minimized, for two reasons. First, minimum 
heating of itself is intrinsically advantageous. Second, it gives one ideal situation in terms of nose radius 
and ablation rate at each flight condition for which we can examine some of the other questions of 
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interest, thus affording some economy in the range in which other parameters need to be varied We 
obtain an appreciation of what may be achieved under advantageous conditions. 

Other questions have to do with the effect of mass addition on related quantities. For example, 
we may expect that mass addition will alter the flow-field structure and change the standoff distance At 
once the question arises, What is the effect of altered standoff distance on radiative hearing? Moreover, 
the species added to the flow field may be expected to radiate differently from air. The question is, How 
important might this effect be? Some ablating materials suffer from lack of physical strength, so we are 
interested in the magnitude of viscous shear stresses at the surface and how they compare with pressure 
stresses. . We are, of course, concerned with the relative importance of radiative to convective heating 
because it bears on the type of surface useful for heat-shielding, and it indicates where improvements in 
our knowledge are more important — in gaseous radiation emission properties or in total thermal 
conductivity of the gas. Ablation rates and effective heat of ablation are important to determine at 
speeds greater than 30,000 ft/sec. 


Because much of the experimental work on mass addition is performed in ground-based facilities 
for which both the entry object and the environment must be modeled, we wish to examine the problem 
of scaling mass addition effects. 

Finally, because of interest in pointed or very slightly blunted entry bodies, low Reynolds num- 
ber (based on nose radius) effects or external vorticity effects have become important. This problem will 
be examined briefly. 


GENERAL FEATURES OF SOLUTIONS 


In this study we are concerned with both the structure of the stagnation-region flow field (which 
will be obtained from solutions of the flow-field equations) and with quantities derived from flow-field 
solutions. The details of the method of solution of the flow-field equations are contained in Chapter 7 
and in references 1 and 2. 

Briefly, the conservation equations for mass, momentum, and energy are solved in the stagnation 
region of blunt bodies from the body surface to the shock wave. Momentum transport by viscosity and 
energy transport by conduction, diffusion of reacting species, and emission and absorption of radiation 
are included in the integro-differential equations. 

The thermodynamic and transport properties (radiative transfer sometimes excepted as discussed 
subsequently) of the mixture of air and injected gases are assumed to be those of partially dissociated 
and ionized air in chemical equilibrium (ref. 2). This assumption should lead to reasonable results even 
if the injected gases become ionized, if the ablation products are nitrogen, oxygen, and carbon com- 
pounds, the reason being the similarities among these species and their compounds. That is, the atomic 
weights of monatomic species are nearly alike, as are the molecular weights of diatomic species More- 
over, the dissociation energy of C0 2 is close to that of O 2 ; that of CN like that of NO; and CO like N 2 ; 
while CO, N, and O all have nearly the same first ionization potential (which is not radically different ' 
from that of C). Actually, the bulk of the injected species will be near the wall of the body where the 
temperature and degree of ionization are suppressed, and the argument of Lees (refs. 3 and 4) may be 
employed. For nonionized gas mixtures, Lees has shown that it is not necessary to understand the 
extraordinarily complicated details of the chemical interaction between the atmospheric gas component 
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and the vaporized surface material insofar as convective energy transport is concerned as long as the 
molecular weights and collision cross sections of the injected gas and air are of comparable size, or that 
the ratio of mass diffusivity to thermal diffusivity (Lewis number) of the mixture is near unity. 

On the other hand, radiative transport by the mixture of injected vapor and air may be signifi- 
cantly different from that of air alone. This behavior is allowed for in the governing equations by 
including a quantity proportional to the injected species concentration in the Planck mean mass absorp- 
tion coefficient expression and a diffusion equation to account for this species. In the results to be pre- 
sented, radiation from injected gases is specified to be like that of air except where stated to the contrary. 

The results of the analysis have been tested against those by others where possible in order to 
establish the validity of the method. In reference 2 it was shown that the calculated structure of the flow 
field exhibits both an isoenergetic shock layer and a boundary layer in the low speed (such that energy 
depletion by radiation is negligible), high Reynolds number regime as it should. The analysis repro- 
duced the well-known effects on convective heating of the assumption p|i = const when p(i was 
artificially set constant. The shock standoff distances predicted by this method agree with those pre- 
dicted by other methods (e.g., ref. 5). When radiation coupling is negligible, it produces convective 
heating results at high speeds (up to 50,000 ft/sec) that agree with the boundary-layer results of refer- 
ences 6 and 7 which use the same transport properties. It will be seen subsequently that in the low 
Reynolds number regime the method leads to flow fields which exhibit the expected shock-layer vortic- 
ity structure and the corresponding enhanced surface shear stress and convective heating rates. One 
additional test of the method is shown in figure 8-1. The solid lines are enthalpy profiles across the flow 
field for the nose radii and flight conditions noted, and the symbols are the results of K. K. Yoshikawa 
(ref. 8), corresponding to the one-dimensional flow of radiating air behind the shock wave. It is seen 
that both analyses show that the shock layer is nonisoenergetic and both give the same results for that 
half of the flow field nearest the shock. The present analysis shows lower enthalpy and larger enthalpy 



Figure 8-1- Comparison of enthalpy profiles with those of reference 8. 
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gradients as a result of convective transport in that half of the flow field near the body, because neither a 
body nor energy transport by conduction was included in the one-dimensional analysis of reference 8. 

Before mass addition effects are discussed, a brief comment should be made regarding terminol- 
ogy. Throughout the rest of the paper, comments relating to mass addition and injected gases, in gen- 
eral, apply to both forced mass addition (transpiration) and natural mass addition (thermal erosion or 
ablation). When our comments are specialized to thermal erosion only, the word ablation will be used. 


EFFECTS OF MASS ADDITION ON HEATING RATES 


Radiative Heating 

Mass addition can affect radiative heating in two ways: by altering the temperature and structure 
of the flow field and by adding species which radiate differently from air. However, results of numerous 
flow-field solutions 1 with mass addition show that (except for a combination of low Reynolds number 
and strong injection of gases which radiate more strongly than air, discussed subsequently) radiative 
heating is much less affected by mass addition than is convective heating. 

Mass addition of gases like air tends to thicken the shock layer, thus tending to enhance radiative 
heating, but it also tends to cool the flow field — inhibiting radiative heating. In almost all of the solu- 
tions, the net effect of injection of a species which radiates like air was to increase radiative heating 
modestly. For example, for a flight speed of 50,000 ft/sec, a nose radius of 0.25 ft, a shock-layer pres- 
sure level of 1 atm, and a surface mass flux of 13% of the free-stream mass flux (f w = -1.5), the radiative 
heating was enhanced about 24%. 

Moreover, it was shown in figure 9 of reference 2 that for high Reynolds number, the bulk of the 
injected species remains close to the vehicle surface where the temperature is low compared with that 
behind the shock wave. Thus the radiant flux at the surface is enhanced only 7% by a gas which emits 
50 times as strongly as air injected at the same rate (f w = - 0.4, or mass addition rate 2% of the free- 
stream mass flux; same flight condition as above with R = 1 ft). 

On the other hand, it will be shown subsequently (in the discussion of fig. 18) that if the 
Reynolds number is low, injected species will also be present in the hot part of the flow field and need 
not radiate much more strongly than air in order to have an appreciable effect on radiative heating. 


Convective Heating 


Convective heating is very strongly influenced by mass addition. In the flight regime for which 
air is dissociated but not ionized, many studies (ref. 9) have employed a linear approximation relating \j/ 
(the ratio of convective heating with mass addition to that without) to the product of mass addition rate 
and driving enthalpy divided by the convective heating without mass addition. The linear \|/ was based 
upon an empirical correlation of experimental and theoretical transpiration results from references 10-12. 
Swann (ref. 13) and Swann and Pittman (ref. 14) obtained a quadratic expression for \j / in terms of the 


1 A11 of the results of this section of the paper conrespond to wall temperatures between 1500 and 3000 K. 
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above variables by an empirical fit of the results of ideal gas boundary-layer solutions of Beckwith 
(ref. 15). 


In the present analysis, the influence of mass addition on convective heating was obtained from 
the results of the ionized radiating flow-field solutions (both q c and q Co , which are convective heating 
with mass addition and that without, respectively, were obtained from the solutions). The parameter y 
for speeds of 40,000 and 50,000 ft/sec is shown in figure 8-2 as a function of f w , the dimensionless 
stream function at the wall. The quantity f w is proportional to the mass addition rate by the relation 
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Figure 8-2 - Separate correlations of results of blowing on convection at TT = 40,000 and 50,000 ft/sec. 

For the moment, attention is directed to the solid curves of the figure. (The dashed lines are an applica- 
tion of the linear y approximation and will be discussed in a subsequent section of this chapter.) Each 
solid line correlates results for one flight speed, various nose radii, and various shock-layer pressure lev- 
els, as can be seen from a comparison with the plotted symbols. Each of these curves is represented by 
the exponential 2 


2 Some experiments by Vojvodich, Pope, and Dickey of Ames Research Center at conditions corresponding to subsatellite 
speed indicate that y for strong ablation of some materials may approach an asymptote different from zero, possibly of the 
order of 10 -1 (this effect was remarked upon in ref. 16). An appropriate form of y for that case would be 
y = a + (1 - a)e~ b ( - & v > n , where a is the value of the asymptote. The asymptote does not appear to be caused by wall 
temperature effects alone. In the examples of figures 2 and 3, no asymptote other than zero could be distinguished even 
though wall temperature was changed from 1500 to 3000 K. For example, at 40,000 ft/sec flight speed, 1-atm shock-layer 
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( 2 ) 


\j/ = e -b(-fw) n 


where b is a function of flight velocity alone and n is 3/2. Similar curves can be constructed for other 
flight velocities. Note that a given f w is less effective in retarding convective heating at higher speeds. 

We can correlate results over a wide range of conditions by expressing b as a function of H , 
where II is flight velocity in units of 10,000 ft/sec. That is, results of 36 mass addition solutions for 
nose radii 3 ranging from 0.01 to 5 ft, flight speeds from 30,000 to 50,000 ft/sec, and shock-layer pres- 
sure levels from 0.1 to 10 atm are correlated in figure 8-3 by use of equation (2) with 



Figure 8-3.- General correlation of result of blowing on convection. (The points are identified in 

table I.) 


pressure level, 1-ft nose radius, and f w = -1.0, both q c and q Co changed as T w was changed from 1500 to 3000 K, but their 
ratio v remained the same to four decimal places, 0.0626. It is conceivable that molecular weight of surface vapors may 
have something to do with the asymptote. For example, the vaporization temperature of Teflon is low enough that surface 
vapors may have a molecular weight of 100 rather than that between 16 and 30 for the air injection (or for that matter, 15 for 
vaporizing phenolic nylon) under consideration. This large a disparity in molecular weight may be significant and would 
tend to raise y for the higher molecular weight gas. The presence of an asymptote for ablation is expected for other reasons; 
that is, if y — » 0, q c -> 0. For test models that are small, qp is also negligible. Thus in the extreme, the heating which causes 
ablation vanishes and some asymptotic value of \\i, or finite value of q c is necessary to initiate ablation. 

3 Two points should be mentioned in connection with the small nose radii. First, the correlation holds for examples for which 
there is strong vorticity in the entire flow field as long as q C(J also includes the external vorlicity effect. Second, the 
chemical equilibrium assumption is somewhat in doubt for the small nose radii, a point which will be discussed subsequently. 
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TABLE L- CONDITIONS FOR POINTS SHOWN IN FIGURES 8-3 AND 8-4 


Group 

Point 

Ps, atm 

u 

R, ft 

n 

Comments 

1 

1 

1.0 

m 

0.01 

0.1 



2 

1.0 

Eg 

1.0 

.i 



3 

1.0 


1.0 

.2 


2 

4 

msm 

mm 


A 



5 

iEl 

10 


.45 



6 


5.0 

.1 

.5 




■89 

5.0 

.25 

.5 




mm 

5.0 

.5 

.5 




Ktn 

5.0 

msm 

.5 



BSIb 

1.0 

5.0 

Wm 

.5 



■Dl 

.1 

5.0 

EH 

.5 


3 

12 

10.0 


.i 

ma 



13 

10.0 

Ks 

.25 

■9 



14 

1.0 

BE 

1.0 

si 



15 

1.0 

EE 

1.0 

■9 



16 

1.0 

Bi 

1.0 

.75 


4 

17 

1.0 

EH 

.1 

|{!| 

17 and 18 appear as one point on figures 8-3 


18 

1.0 

Wm 

.25 

09 

and 8-4 


19 

1.0 

5.0 

.5 

HI 



20 

1.0 

5.0 

1.0 

1.0 



21 

.1 

5.0 

5.0 

1.0 


5 

22 

1.0 

5.0 

1.0 

1.25 



23 

1.0 

4.1 

.01 

1.0 



24 

1.0 

4.1 

1.0 

1.0 



25 

10.0 

4.0 

.1 

1.0 

25 and 26 appear as one point on figures 8-3 


26 

10.0 

4.0 

.25 

1.0 

and 8-4 


27 

1.0 

4.0 

1.0 

1.0 


6 

28 

msm 

3.0 


msm 



29 

189 

5.0 

.1 

■ 

29 and 30 appear as one point on figures 8-3 


30 

mm 

5.0 

.25 

m 

and 8-4 


31 

in 

5.0 

.5 

1.5 

31 and 32 appear as one point on figures 8-3 


32 

msm 

5.0 

ill 

1.5 

and 8-4 


33 

.i 

5.0 

E89 

1.5 


7 

MM 

10.0 

4.0 

.25 

lii 

34 and 35 appear as one point on figures 8-3 


Km 

1.0 

4.0 

1.0 

mWM 

and 8-4 

8 

36 

1.0 

5.0 

1.0 

2.0 
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(3) 


in n 6 ) 

and n = 3/2. This formula is useful for extrapolation to flight speeds above 50,000 ft/sec, but cannot be 
used for speeds much below 30,000 ft/sec because it changes sign at about 26,000 ft/sec. The same 
mass addition solution results are correlated by letting: 

b = 0.706+ 1.6 n-0.2811 2 (4) 

and n = 3/2. The result is shown in figure 8-4. This formula can be used to extrapolate to speeds below 

30.000 ft/sec, but cannot be used for speeds much above 50,000 ft/sec, because it changes sign at about 

61.000 ft/sec. ° ° 



Figure 8-4 - General correlation of result of blowing on convection. (The points are identified in 

table I.) 

The flight condition, nose radius, and value of f w for each point shown in figures 8-3 and 8-4 
are listed in table I. Groups of points are numbered consecutively from left to right on the figure. Points 
within a group are numbered consecutively from top to bottom. 

These same results are also compared in figure 8-5 with the linear \j/ approximation of reference 
9 and quadratic V approximation of reference 14. The points on the figure correspond to our solutions 
listed in table II, where now the points are simply numbered from left to right in figure 8-5. In the fig- 
ure, if the constant 0.49 in the linear approximation corresponding to Teflon (which would be 0.5 for 
phenolic nylon (ref. 9)) is changed to 0.6 (mentioned by ref. 9 as obtained from refs. 10-12), the fit is 
improved for the initial part of the data out to about H s m/(-qc 0 ) equal to unity. Clearly, however, it 
cannot be made to fit the results beyond unity. 
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Figure 8-5 - Present results compared with other correlation formulas. (The points are identified in 

table II.) 

The quadratic expression of reference 14 fits the present result well for H s m/(-q Co ) equal to unity 
also, but not beyond. That approximation is set to zero when H s iiV(-qco) is 2.5. However, three points 
are shown to the right of 2.5 for which \}/ ^ 0. Subsequently, during the examination of special condi- 
tions which minimize total heating rate, it will be important to have a simple correlation which must dif- 
fer from zero at the high mass addition rates, for if y were zero, the radius which minimizes total heat- 
ing would be that which minimizes radiative heating, namely, zero. The exponential correlation (eq. (2)) 
will be especially useful in that regard. 

The two approximations were obtained from relatively low-speed results for which the air was 
either dissociated or inert, and it is not surprising to see that they do not fit the present results which 
include (among other differences) thermodynamic and transport properties of partially ionized air and 
are extended to higher values of H s m/(- q Co ). 

The simple expression, equation (2), describing the effects of mass addition on convective heat- 
ing at high speeds can now be applied with other information to the special case of mass addition by 
ablation. 
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TABLE H.- CONDITIONS FOR POINTS 
SHOWN IN FIGURE 8-5 
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ABLATION 


General Relations for Ablation Quantities 


Subliming ablators are attractive for heat protection because of heat absorption due to vaporiza- 
tion and the heat blockage effect in the boundary layer. Moreover, high-temperature subliming ablators 
have still another asset — the rejection of heat by reradiation. 


On the other hand, at conditions corresponding to subsatellite speeds, subliming ablators have at 
least one liability from the heat-rejection point of view. A theoretical study of Scala (ref. 17) shows that 
reactions between air and a graphite surface impose a significant heat load on the vehicle. Theoretical 
studies by Hartnett and Eckert (ref. 18) and Cohen, Bromberg, and Lipkis (ref. 19) show enhanced con- 
vective heating caused by gas-phase reactions between air and ablated vapors. Experimental results of 
Vojvodich and Pope (ref. 20) confirm that both heterogeneous and homogeneous combustion between 
air and charring ablators impose a heat load comparable to the net convective heating (q c w) at low lev- 
els of shock-layer pressure (10-3 to 10~ 2 atm) and driving enthalpies up to about 8000 Btu/lb m . How- 
ever, they show that the relative importance of combustion diminishes with increasing driving enthalpy 
and increasing pressure level. Very likely, the reason for diminished importance of surface reactions 
between air and ablation material is that higher injection rates prevent air from reaching the surface at 
the more extreme conditions. An analogous phenomenon was studied theoretically by Chung (ref. 21) in 
which he showed that heterogeneous recombination reactions are inhibited by air transpiration at a cold 
wall, preventing dissociated shock-layer air from reaching the surface. Further, the relative importance 
of energy release by gas phase reactions between air and injected species is diminished, probably 
because of the increased energy release by recombination reactions of air components themselves at the 
more severe conditions. 

Thus we assume that at the higher levels of shock-layer pressure (1(H to 10 atm) and higher 
enthalpies (20,000 to 50,000 Btu/lb m ) with which the present study is concerned, one need not sort out 
gas phase combustion reactions from other recombination reactions. So, we appropriate Lees' argument 
mentioned previously and neglect the details of the combustion reactions, but consider for practical pur- 
poses that their effects are included implicitly in q^, the convective heating results. 

. high flight speeds, radiative heating must be included with convective heating as causing 
ablation. The ablation rate for a given vaporizing material is related to the actual total heating flux, q T , 


v 


X (piVi)w 

i=l 



C Pi dT + h ivj 


Qj 9c ~ Qr 


(5) 


where Vj is the absolute velocity of species i and the minus sign arises from the convention that posi- 
tive flux is outward from the surface (that is, q means flux; it will be represented by a negative number 
the wall is receiving heat), and the summation is over the ablation products appearing at the surface. 
Tiie lower limit on the integral is T m , the material temperature before heating began. The term h; is 
the enthalpy of vaporization of species i. Now we define the intrinsic heat of ablation as 
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( 6 ) 



(Pi^Dw 


i=l 


The definition of the mass averaged velocity applied at the wall is 


PwVw — (pairV air)w + | ^ Pi^i 
, i=l 


(7) 


In order for air not to penetrate the surface, (pairVairXv = 0. Then combining equations (5)-(7) and (1) 
yields 


- q T = mh a (8) 

In equation (5), q r is the net radiative heat flux accepted by the wall. It is a combination of the 
accepted incident radiative flux from the gas and the reradiated flux from the wall; thus 

% = qr g + (9) 

where e w is the surface emissivity, and a is the Stefan-Boltzmann constant. 

To express qp in a simple way, assume that (1) there is no coupling between radiative and con- 
vective heating rates (a very reasonable assumption for the application to moderate sized bodies made 
subsequently) and (2) q rg is not a strong function of -f w or mass addition rate (which will be verified 
for conditions of special interest in the next section) and that (3) 

qr g = -<x w BR (10) 

where B is a constant for a given flight condition and a w is the surface absorption coefficient. 
Numerical solutions indicate that, for purposes of the estimates in this section, equation (10) is a good 
approximation even though the flow field is nonisoenergetic, as long as B is obtained from the non- 
isoenergetic solutions at a given flight condition. Because no coupling is assumed (1) we can use the 
simple, no-blowing, convective-heating correlation (which, by the way, excludes external vorticity 
effects) of Hoshizaki (ref. 6), which can be put in the form 

^ = q Co = -4.03xlO-5 (2 e)^n° 19 (H s -H w ) '\J~^ (11) 

where the units on 4.03x10-5 are (lbf) 1 ^ 2 sec/ft 3 ^ 2 and the units on q Co are lbf/ft 2 sec. Combining equa- 
tions (2), (5), and (9)-(l 1) leads to an expression for total heating rate 
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(12) 


-q T = a w BR - e w T^ + 4.03xl0- 5 (2e) 1/4 n° 19 (H w - H w ) e - b <- f w> n 
It is convenient to define 


A =- 


a w B 
haV PsM-sU (j + 1) 


4.03xlQ- 5 (2e) 1 ( 4 II 019 (H s -H w ) Vp7 ; 
h a Vp s |i s U(j + 1) 


D = 


ae 


W 


h a VpsH s U(j + 1) 


(13) 


In terms of these quantities, equation (12) becomes 


— q-p — h a V psHsU(j 1) 


AR+^e^-^-Df, 

fR 




(14) 


Now if equations (1), (5), and (8)-(l 1) are combined, a transcendental expression is obtained for 
the natural blowing rate parameter f w in terms of R, the flight conditions, and material properties: 

-f w = AR 3/2 + Ce-^w) 11 - DR 1/2 I^ ( 1 5) 

Another quantity of interest is the effective heat of ablation (where the subscript o means with- 
out mass addition). 


h e ff = 


-QTo 


m 




m 


(16) 


If we assume 

qr 0 = qr ( 17 ) 

which is approximately true for the moderate ablation rates which will concern us in this and the next 
section (and will be demonstrated at the end of the next section), it is simple to show that 

Jt = i_[i_e-K-fw)']^2. (18) 

q-To I- J q T 0 

Combining equations (8), (13), (16), and (18) leads to 


131 



heff=ha {l + ^[ 1 - e - bW "' ) | 09) 

where -f w is obtained from equation (15). 

convectiononly! inSm,C ' iVe ‘° Specialize some of ,hese ablado " relationships still further-to the case of 

Ablation Due to Convection Only 

im fl l gbt speeds g r ® at e r than satellite entry speed, both radiative and convective heating are 

abIanon rates ’ even for sma11 nose radii (which will be shown subsequently) 
However, if for the moment convection is assumed to be the only heating mode, and both radiation and 

We ^ '? S0mc ' nlere sting comparisons. pfr these conditionsZth eqna 

tions (15) and (19) reduce to simpler forms, thus 4 

-f w e b (~ f w) n = c , 9m 


^=e b (- f w) n =l 


respectively. So we see that for the convective case, both the ablation rate in terms of f w (eo (20)) and 
effective heat of ablation (eq. (21)) for a given material are independent of nose radius and amS den- 
sity (or pressure level) but depend only on flight speed (because b is a function of H alone (eq (4)) 
and C is essentially a function of XT because in equation (13) 64 ' 

(H s - H w ) Vps/psPs * (U 3 /^) Vpoc/PsPs 

Adams S 3 Vef 7 WCak [ UnCti0 " ° f a T bient densit y)- Lees (ref. 22) and Bethe and 

Adams (ref. 23) reached a similar conclusion for melting and glassy ablators at subsatellite speeds- that 
is, velocity is the important parameter. p eus ’ mat 

respect ? *5^?™ C °7i re thi l fw andbeff result obtained from V ‘hat is exponential with 

respect to (-f w ) (eq. (2)) with those obtained from y that is linear with respect to (-f w ) of reference 9 

The \|/ of that paper written in terms of f w by use of equation (1) is 

v = i - ft ,,/ iwio + 1) 

q Co V r w l > 

ptp'eO air " ^ COrreSP ° ndin8 f » < fOT * ha ‘ 
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(23) 


4.03xlQ-5(2e) 1 / 4 n O19 (H s - H w ) / p s 

[h a + P(H S - H w )] V psH s U(j + 1) 


or 


l + (P/ha)(H s -H w ) ^ 

while h e ff for the linear \| / would be 

!gl= 1+ i(H s -H w , (25) 

The comparison of the natural f w for Teflon as predicted by the two results (eqs. (20) or (22) in 
(24)) is shown in figure 8-6 (using p = 0.49 and h a = 2.38xl0 7 ft 2 /sec 2 from reference 9 for Teflon). 4 It 
is seen that the two results are in close agreement and that both depend strongly on velocity, but not on 
pressure. Similarly, the results of h c ff/h a as a function of velocity are in close agreement for convection 
only (using eqs. (21) and (25)) as shown in figure 8-7. 



Figure 8-6.— Wall blowing parameter for natural ablation of Teflon at stagnation region — convection 

only. 


In this application and throughout the rest of this chapter, the values of b and n for a given flight condition were obtained 
by passing equation (2) through results from the two flow-field solutions having the highest values of — f w . Numerically, b 
and n are slightly different from equations (3) and (4) for b and 3/2 for n obtained for the general correlation. 
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u 

Figure 8-7.- Effective heat of ablation of Teflon (convection only). 

The agreement between the two methods might indicate that e b T f w) T of equations (20) and (21) 
equals 1 + (P/ha)(H s — H w ). At first glance one might suspect that these natural blowing rates are small 
enough to be on the initial part of the \j f curve where the exponential can be replaced by the first two 
terms of its series. It turns out that this is not the case, as can be seen in figure 8-2. The solid curves and 
dashed lines correspond to the exponential and linear \|r (to eqs. (2) and (22)), respectively, for the flight 
speeds shown. Although agreement between the solid curve and dashed line is better at the lower 
speeds, in neither case is the solid curve well represented by the dashed line except at the one point 
where they intersect. That intersection just happens to occur very near the natural value of f w in each 
case. Thus the apparent agreement in the results does not imply any general agreement in the function 
\ji but, rather, is considered to be fortuitous. 

Now we turn again to the case of ablation caused by both radiative and convective heating to 
examine some special conditions. 


Conditions at Minimum Heating Rate With Ablation 

Because of the many combinations of variables, parameters, and phenomena associated with the 
ablation problem, it is convenient to seek an optimum condition, in terms of minimizing total heating 
rate or total mass loss rate of a given material at each flight condition, and then present some of the other 
quantities of interest correspondingly. 

Nose radius and ablation rate for minimum heating rate - Equation (14) shows that at a given 
flight condition, convective heating rate becomes large with small R whereas radiative heating rate 
becomes large at large R. Thus there is an intermediate value of R for which total heating rate is a 
minimum. This is illustrated graphically for one flight condition in figure 8-8. The family of solid light 
curves is calculated by use of equation (12) for specified values of f w , and represents approximately the 
flow-field solution results. The family of dashed curves is calculated for Teflon by use of equations (8) 
and (1), also for specified values of f w (the physical properties used for Teflon and phenolic nylon are 
listed in table III). The intersection of a line of each family corresponding to one value of f w denotes a 
combination of qj, R, and f w that satisfies both the flow-field results and the material behavior. The 
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heavy solid curve is the locus of such intersections. Its minimum gives the optimum nose radius R*, 
which minimizes total heating rate (and total ablation rate) for Teflon for this flight condition. Thus R* 
is 0.109 ft and the ablation rate is given by -f w = 1.94 in this example. The right branch of the heavy 
curve shows an interesting result. That is, for a nose radius half again as large as R*, convective heating 
has been essentially eliminated by strong ablation, and therefore q T increases linearly with R and is 
independent of ablation rate (-f w = 3 line coincides with -f w = ~), as it should be for radiative heating 
only, according to the approximation of equation (10). Correspondingly, ablation is caused by radiative 
heating alone and its rate must be increasing in proportion to qj, and thus R, in accord with equation 
(8). This points out the potential importance of reflecting ablative surfaces for heat protection for radii 
greater than R*. 

FLOW FIELD 

RESULTS, EQ. (12) 



Figure 8-8.- Optimum nose radius for Teflon at U = 50,000 ft/sec, p s = 1 atm (a w = e w = 1). 

Analytically, R* and f w for a specific material and flight condition are obtained as follows. The 
partial derivative of qj with respect to R (obtained by differentiation of eq. (14) noting that f w is a 
function of R by eq. (15)) is set to zero. After some algebra, the expression 

rv-b(-f w ) n A , n 

fw + [bn(-f w ) n + 3] - DT^R 1 / 2 = 0 (26) 

is obtained. The simultaneous solution of equations (15) and (26) yields the optimum nose radius and 
natural ablation rate (f w ) for a given material and flight condition. 

It can be noted that for no reradiation (DT^R 1 / 2 negligible), equation (26) is uncoupled from 
equation (15) and the former can be solved directly for f w after which the latter can be used to obtain 
the optimum R. Thus we obtain the result that f w is independent of the radiative properties of the gas 
expressed by A in equation (13). Although f w depends only on the convective heating properties of 
the gas (expressed by C in eq. (13) which appears in eq. (26)), it differs from the f w for convection 
only (eq. (20)) for the following reasons. If equation (15) is satisfied by AR^ * 0 by use of the f w 
obtained from equation (26), it cannot be satisfied for AR^/ 2 = 0 by the same f w . Moreover, if radia- 
tion is zero, we cannot use equation (26) to calculate f w , but must revert to equation (20). 
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TABLE III.- PROPERTIES OF ABLATING SURFACES 3 


Material 

h a , ft 2 /sec 2 

T W ,K 

£&) 

a w 

Teflon b 

2.38x107 (9) 

1000 

1 and 0.5 

1 and 0.5 

Phenolic nylon 

3.755x107 (9) 

3700 (50) 

0.6 (50) 

0.6 (50) 


a Values in the table are estimates made from information in the 
references shown in parentheses. 

b The temperature shown for Teflon is higher than the 800°F (or 
700 K) given by reference 9. However, it is of no consequence 
because reradiation from the higher temperature is still negligible. It 
should be noted that in applying our \|/ results to Teflon, we have not 
made any allowance for the presence of a finite asymptote (see foot- 
note 2). If indeed a finite asymptote does exist at hypervelocity, we 
must consider our Teflon results to apply instead to a material which 
has the properties shown above but which does not have a finite \|r 
asymptote. 


For radiation different from zero but negligible reradiation, heff is independent of the radiative 
flux because f w used in equation (19) is independent of radiative properties. The same is true for \|/ 
because of equation (2). On the other hand the optimum R obtained from equation (15) depends on the 
radiative properties. Of course, m depends on the radiative properties because it depends on R 
(eq. (1)), and the same is true for both q c (eq. (11)) and q r (eqs. (9) and (10)). 

The optimum nose radius with convection, radiation, and reradiation obtained from equations 
(15) and (26) is shown as a function of flight speed in figure 8-9 for Teflon and phenolic nylon. It can 
be seen that the Optimum nose radius diminishes with increasing shock-layer pressure at a given flight 
speed. On the other hand, the optimum nose radius increases as flight speed is diminished, shock-layer 
pressure level being constant. These trends can be related to actual entry trajectories, by noting (ref, 2, 
fig. 1) that typical trajectories consist of essentially a path of increasing shock-layer pressure at constant 
velocity followed by a path of diminishing velocity at constant pressure. Now if the major heating 
occurs at constant velocity (the case of plunging probes), total heating rates could be minimized only by 
artificially tailoring the nose shape (for example, by pushing concentric rods of progressively smaller 
radii out the front of the vehicle in a programmed sequence— a refinement of a suggestion by H. J. Allen 
(ref. 24)). On the other hand, if the major heating occurs at constant shock-layer pressure (typical of 
g-limited entry), the problem of minimizing total heating rates is simplified because R grows by 
ablation naturally in the direction of the growing optimum. 



D °°l TEFLON 

PHENOLIC NYLON 

p s (atm), a w = e w 



Figure 8-9 - Optimum nose radius for ablation. 

The wall absorptivity, a w , is important to both charring and noncharring ablators because it 
influences the amount of radiant heat accepted from the gas cap and thus R*. In figure 8-9, if a w is 
diminished by 50% for Teflon (noncharring), R* is increased by about 50%. 

The wall emissivity, e w , is important to high-temperature charring ablators for which reradiation 
is an important heat-rejection mechanism. For phenolic nylon, reradiation is partly responsible for a 
larger optimum nose radius than that of Teflon (for which reradiation is negligible because of its low 
vaporization temperature) — by about a factor of 3 in the speed range 40,000-50,000 ft/sec as can be seen 
in the figure. The phenolic nylon calculation for p s = 1 atm was not extended to lower speed because 
of the uncertain wall temperature; that is, the wall temperature becomes a function of the heating rate at 
less severe heating conditions. 

The ablation rate, in terms of f w , corresponding to optimum heating rate conditions is shown as 
a function of flight speed for various pressure levels and absorptivities for Teflon and phenolic nylon in 
figure 8-10. The heavy lines and symbols represent minimized total heating-rate conditions, while the 
light lines represent the result for convection only (shown previously in fig. 8-6 for Teflon). It is impor- 
tant to point out that the difference between the total-heating heavy lines and convective-heating light 
lines does not represent the contribution of radiative heating. When minimum total heating rate is con- 
sidered, there is a complete rearrangement of convective and radiative contributions. The result is that 
the radiative contribution is usually larger (sometimes much larger) than the convective contribution 
(which will be demonstrated subsequently in fig. 8-15). The result of the sum of the two readjusted 
heating components on f w is shown by the heavy lines in figure 8-10. 
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Figure 8-10.- Blasius wall-blowing parameter corresponding to ablation with optimum nose radius. 


It is especially interesting to note that for Teflon, optimum f w is a very weak function of both 
absorptmty and pressure level. This will have consequences in figure 8-12, where it will be remarked 


The ratio of the surface mass flux to free-stream mass flux is related to fw 
number. Thus 


by the Reynolds 


f w -/S± 
>ooU V eR 


1 ) 

eRe 


(27) 


is Reynolds nunter P- UR /^s- It should be mentioned that e 1 ^ is almost invariant 
with flight speed (between 30,000 and 50,000 ft/sec) at a given level of shock-layer pressure. It only 
vanes at the worst from 0.233 to 0.278 as pressure level changes from 0.1 atm to 10.0 atm. 

The mass flux ratio of Teflon is presented in figure 8-11 for optimum conditions. The location of 
e point corresponding to 0.1 atm indicates that relationships at optimum conditions are not systematic 

Z & (SI? P a cnLJJ 6 / SUrf o Ce maSS fluX VarieS fr0m 15 t0 25% of the free-stream mass flux between 
JO.UOO and 50,000 ft/sec. Some approximate calculations show that it is reasonable to extrapolate this 

result to 70,000 ft/sec. The result is that the mass flux ratio is still less than 0.5 for optimum conditions, 
this is m sharp contrast to the values in excess of unity for entry of some meteors for which the nose 

radius is very different from R* (Pribram meteor, for example, which likely has a radius of the order of 
a meter, ret. 25). 
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Figure 8-11- Mass addition ratio for optimum conditions (Teflon, a w = 1). 

In fact, it is apparent from the present optimum results for U = 50,000 ft/sec and p s = 1 atm 
that if nose radius is increased from the optimum of 0.109 to only 0.5 ft, the mass loss ratio will exceed 
unity as a result of radiative heating alone. 

Now that we have the nose radius which minimizes total heating rate, and the corresponding 
ablation rate, we can examine several other interesting quantities corresponding to these conditions. 

Effective heat of ablation- The ratio of the effective heat of ablation to the intrinsic heat of abla- 
tion is calculated by use of equation (19) and is presented in figure 8-12. The ratio depends strongly on 
velocity, but is relatively insensitive to pressure level (altitude) and, for the case of Teflon, surface 
absorptivity. The last is a consequence of the insensitivity of f w to a w shown in figure 8-10, and the 
fact that h e ff is very strongly dependent on f w in equation (19). The ratio h e ff/h a subject to combined 
radiative and convective heating is less than that obtained from convection alone, which can be seen by 
comparing figure 8-12 with figure 8-7. 

In spite of the fact that h e ff/h a in figure 8-12 is smaller for phenolic nylon than for Teflon, the 
actual h e ff for phenolic nylon is larger. At 50,000 ft/sec, the ratio of effective heat of ablation of phe- 
nolic nylon to that of Teflon is 1.13, which includes the reradiation effects. (It should also be remem- 
bered that the phenolic nylon is not ablating as rapidly as Teflon for these optimum conditions.) 

Effect of mass addition on standoff distance - The ratio of standoff distance with optimum mass 
addition (and nose radius) to that without mass addition, 5 0 , was obtained by solving the flow-field 
equations using values of R* and f w obtained from figures 8-9 and 8-10. The result is shown in fig- 
ure 8-13 for Teflon. For speeds below 30,000 ft/sec, where mass addition rates are low (3% of free- 
stream mass flux, fig. 8-11), the standoff distance with blowing is actually less than that without blowing 
(ratio is about 0.94). This interesting effect occurs at very low mass addition rates over a broad flight 
range and may be attributable to a cooling effect that increases the flow-field density (and diminishes 
standoff distance) more than enough to overcome the space required for the increase in mass flow in the 
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Figure 8-12 - Effective heat of ablation corresponding to natural ablation for optimum conditions. 


p $ , atm 



Figure 8-13 — Ratio of standoff distance with blowing to that without blowing for optimum conditions 

(Teflon, a w = 1). 

flow field. At higher speeds, and thus higher ablation rates, the standoff distance is enhanced by abla- 
tion. It is seen that for optimum conditions, the shock layer is thickened by about 50% at a speed of 
50,000 ft/sec. 

Effect of mass ad dition on radiative heating - Mass addition can influence gaseous radiative 
heating flux in two ways: (1) by altering the stand-off distance and temperature distribution in the flow 
field (it should be mentioned that gaseous radiation reabsorption is negligible in the regime being con- 
sidered), and (2) by adding chemical species to the flow field which radiate differently from air. To 
examine the first effect, due to alteration of flow-field structure, the foreign species was assumed to 
radiate like air. The resulting incident gaseous radiant flux at the wall as obtained from flow-field solu- 
tions corresponding to the optimum R* and f w of Teflon is compared with the no ablation value in 
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figure 8-14. In the flight regime studied, the maximum effect of the ablation on radiant flux is an 
increase of about 17% for these optimum conditions. This is in sharp contrast with the very large effects 
on convective heating presented earlier and supports the approximations of equations (10) and (17). For 
example, for a flight condition in which radiative flux was changed 17%, ablation diminished the 
convective flux by two orders of magnitude (to 0.8% of its nonablation value). 



u 

Figure 8-14.- Ratio of incident gaseous radiation flux with ablation to that without ablation at optimum 

conditions (Teflon, a w = 1). 

To examine the second effect, that of introducing species which radiate differently from air, a 
(the rate of absorption in the foreign species relative to that of air (chap. 7)) was specified to be different 
from unity. Briefly the result is that if a foreign gas that radiates twice as strongly as air (a = 2) is intro- 
duced at the same rate, the radiative flux is enhanced at the most by only 5% over the air value. If the 
foreign gas radiates 10 times as strongly as air (a = 10), the radiative flux is enhanced at the most by 
about 50% over its airlike value. Finally, the influence of these radiative properties that differ from 
those used in estimating optimum conditions on the optimum conditions themselves is as discussed pre- 
viously. That is, to the extent that reradiation is negligible, there is no effect on f w and h e ff. However, 
the effect on R* is to reduce it by a factor of 1.5~ 2 ^ (or 0.76) and the effect on fn is to increase it by a 
factor of 1.5 1 / 3 (or 1.14) at the worst by virtue of equations (15) and (1), respectively. 

Comparison of ra diative and convective heating - The ratio of the radiative to convective heating 
rate was obtained from solutions of the flow-field equations in which R* and f w corresponded to mini- 
mum heating conditions for Teflon. The result is that radiation exceeds convection by a factor of from 
2 to 9 as shown in figure 8-15. The importance of convection increases as the level of shock-layer pres- 
sure is increased. For phenolic nylon, the results are almost the same without reradiation. That is, q r Jq c 
is greater than unity. However, because of reradiation, the ratio q r /q c is less than unity except at high 
pressures (p s = 10 atm). In any event, the gaseous radiant flux incident on the wall is considerably larger 
than the convective flux for either Teflon or phenolic nylon at optimum conditions. 



Figure 8-15.- Comparison of radiative and convective heating at optimum conditions (Teflon, a w = 1). 
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Effect of mass addi tion on surface shear stress - In figure 8-16, the surface shear stress divided by 
distance from the stagnation point is shown as a function of flight speed for various pressure levels for 
Teflon at optimum conditions. It was obtained from the flow-field solutions by use of the relationship 



V(j + l)psHs 


(28) 


which can be derived by use of the transforms in chapter 7. Generally, wall shear stress increases with 
velocity and shock-layer pressure level. For reasonable values of x, the shear stress is not excessive 
(even at p s = 10 atm, the surface shear stress per foot is of the order of standard atmospheric pressure 
per foot for optimum conditions, and the optimum size is considerably less than a foot). The shear stress 
for phenolic nylon would be slightly larger than that for Teflon because the mass addition rate for the 
former is lower at a given flight condition. 



Figure 8-16.- Stagnation region shear stress at optimum conditions (Teflon, a w = 1.0). 


Of course, without mass addition, the shear stress would be considerably higher. The ratio of 
wall shear stress with mass addition to that without mass addition is shown in figure 8-17 for Teflon 
under optimum conditions. It is noted that the ratio is only 0.035 at TJ = 50,000 ft/sec and p s = 1 atm. 



u 


Figure 8-17.- Ratio of stagnation region shear stress with ablation to that without ablation at optimu> w 

conditions (Teflon, a w = 1). 
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CQmrngnt$ Qn validity .— Our optimum-condition considerations have led to results of small nose 
radii for which the assumption of a flow field in chemical equilibrium may be doubtful. In spite of the 
fact that the chemistry and some of the thermodynamic and transport properties may be grossly in error 
for nonequilibrium flow fields, the basic structure of the flow field (velocity and enthalpy profiles) is not 
expected to be seriously in error. The reason is twofold: First, the analysis in the appendix of refer- 
ence 26 shows that in the absence of transport phenomena, the enthalpy profile is almost unaffected by 
large departures from chemical equilibrium. Second, the velocity and enthalpy profile results of refer- 
ence 27, which include both chemical nonequilibrium and transport phenomena, do not differ in a sig- 
nificant way from those of reference 2 or this chapter for chemical equilibrium (excluding low Reynolds 
number results). 

Moreover, the convective heating is not likely to be very much in error because the wall is 
expected to be catalytic, and thus nonequilibrium convective heating would be essentially the same as 
that for equilibrium since recombination would occur at or near the wall in either case. 

We have noted that at the minimum heating condition, radiative heating dominates. Our radia- 
tive heating estimates could be in error for several reasons. Although the neglect of nonequilibrium 
radiation is a potential source of error for the small optimum nose radii, the present estimate is that it is 
not an important effect as gaged by the results of reference 28. We have already noted that there is some 
change in radiative heating caused by flow-field distortion and the presence of foreign species which 
radiate differently from air. Moreover, there is still considerable uncertainty in the radiative properties 
of air itself (the radiative properties of air given by ref. 29 may be high by a factor of 2 according to 
ref. 30) and in the absorptivities and emissivities of the surface material. 

The effect of underestimating the radiative properties on optimum conditions has been noted in 
the section on effect of mass addition on radiative heating. Briefly, we now examine the effect of 
overestimating radiative properties by a factor of 2. As before, f w , h ef f, and \\r are unchanged (for a 
material which does not reradiate importantly). But R* is enhanced by about 60%, so m and q c are 
diminished by a factor 1 (or 0.79) by virtue of equations (1) and (1 1). So both q r and q c change, 
but in such a way that f w is constant. The error in radiative properties by a factor of 2 has a large effect 
on R* (60%); only a 20% effect on q c and m; and no effect on f w , h e ff, and \jr. 

Thus, in spite of the many uncertainties in the flow-field chemistry and in our knowledge of gas 
properties, the results of the minimized heating-rate study remain meaningful. 


MASS ADDITION AND SCALING 


To a large extent, experimental studies of high-speed planetary entry problems consist of expos- 
ing models to a simulated entry environment in a test facility, such as an arc-heated wind tunnel, a shock 
tube, or a ballistic range. Experimental results are then scaled to the actual flight conditions by one 
means or another. 


Strictly speaking, we cannot expect to scale flow-field profiles at all because thermodynamic and 
transport properties used in the flow equations do not scale. However, in this portion of this chapter, we 
will examine briefly how to scale dimensionless foreign species profiles approximately for forced mass 
addition (transpiration) and then specialize the result to natural mass addition (ablation). 
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Scaling With Arbitrary Mass Addition 


It is well known that in order to scale stagnation-region flow fields in general, the Reynolds 
number should be fixed. In order to scale foreign- species concentration profiles, the mass addition rate_ 
must also be fixed. Or, because of the Reynolds number factor in the expression relating f w andm/pooU 
(and in view of the small variation m e 1/2 noted earlier), we may simply say that in order to scale mass 
addition effects, both f w andm/pooU must be fixed. We will illustrate this by use of flow-field 
solutions. 

The main points of the demonstration are briefly as follows. In figures 18-21, solution scorre- 
sponding to the Conditions shown in table IV are presented in which either, but not both, f w orm/pooU is 
the same between pairs of examples (assuming that thesequantities can be varied at will). These results 
can be compared with those for which both f w andm/pooU (or Re) are constant, shown in figures 8-22 
through 8-25. 



Figure 8-18.- Flow-field profiles; U = 41,000 ft/sec, p s = 1 atm, R = 0.01 ft, f w = -1. 


144 


I'l 1 1 



Figure 8-19.- Flow-field profiles; U = 41,000 ft/sec, p s = 1 atm, R = 0.01 ft, f w = -0.1. 



Figure 8-20.- Flow-field profiles; U = 41,000 ft/sec, p s = 1 atm, R = 1 ft, f w = — 1. 
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Figure 8-21.- Row-field profiles; U = 41,000 ft/sec, p s = 1 atm, R = 1 ft, f w = -0.1. 


TABLE IV.- SCALING FOREIGN-SPECIES EXAMPLES 


Figure 

Ux KT 4 , ft/sec 

Ps, atm 

D 

m/pooU 

R, ft 

Re 

18 

4.1 

1 

1.0 

0.5 

0.01 

1.14xl0 2 

19 

4.1 

1 

.1 

.05 

.01 

1.14xl0 2 

20 

4.1 

1 

1.0 

.05 

1.0 

1. 14x10* 

21 

4.1 

1 

.1 

.005 

1.0 

1. 14x10* 

22 

3 

1 

.3 

.05 

.063 

1.19xl0 3 

23 

4.1 

10 

.3 

.05 

.011 

1.05xl0 3 

24 

4.1 

1 

.3 

.05 

.1 


25 

5 

1 

.3 

.05 

.074 

1.26xl0 3 
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Figure 8-24.- Flow-field profiles; U = 41,000 ft/sec, p s = 1 atm, R = 0.1 ft, f w = -0.3. 



Figure 8-25 - Flow-field profiles; U = 50,000 ft/sec, p s = 1 atm, R = 0.074 ft, f w = -0.3. 

The foreign-species profiles corresponding to these two groups of figures are summarized in fig- 
ures 8-26(a) and (b), respectively. Obviously, the concentration profiles resulting from fixing only one 
of the parameters do not scale, as can be seen in figure 8-26(a). The figure shows that f w essentially 
controls the foreign-species concentration at the wall, while Reynolds number determines its penetration 
into the flow field. Thus, if both f w and Re (or m/p<»U) are fixed, scaling should be much improved, as 
is evidenced by the Cf profiles summarized in figure 8-26(b). Moreover, since fixing f w andnVp«,U 
essentially preserves the Reynolds number (the slight variations in Re in table IV are caused by varia- 
tions in e 1 ^ in eq. (27)), the velocity profiles of the second set (figs. 8-22 through 8-25) are scaled; 
where in the first set (figs. 8-18 through 8-21) they are not. Finally, the enthalpy profiles of the first set 
exhibit less similarity than those of the second set. The latter are summarized in figure 8-27. 
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Figure 8-26 - Foreign species concentration profiles, (a) Either f w orm/pooU fixed. 

(b) Both f w (= -0.3) and m/p<x>U (= 0.05) fixed. 

There are additional features of some of the solutions mentioned above which ^although sec- 
ondary to the argument, are worth comment. In figure 8-18, for a flight condition of U = 41,000 ft/sec 
and p s = 1 atm, a nose radius of 0.01 ft, and a blowing rate f w = -1.0, both momentum (associated with 
u/u s profile) and thermal (associated with H/H s profile) boundary layers are conspicuously absent. 
Indeed, the vorticity (slope of u/u s curve) is approximately constant throughout the flow field. Inter- 
estingly, the flow field is far from isoenergetic everywhere; not because of enthalpy depletion by radia- 
tion but, rather, by convection. The mass addition rate at the surface is half the free-stream mass flux 
(table IV), and standoff distance is 90% higher than the no-blowing value (this is the most exoreme result 
we have in both regards). The foreign species completely permeates the flow field, which is important 
from the point of view of the relative importance of air and ablation species radiation. In this regard, 
results show that the incident gaseous radiant-heating flux at the wall is enhanced by 49% over the 
no-blowing value if air is the injected gas (because of the thickened shock layer), but is enhanced by 
93% over the no-blowing value if a species that is three times as strong an absorber and emitter as air 
(a = 3) is injected at the same rate. If the mass addition rate is diminished by a factor 1/10, the principal 
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Figure 8-27.- Comparison of enthalpy profiles for both f w (= -0.3) and m/p«,U (= 0.05) fixed. 

effect is to greatly diminish the foreign-species concentration (fig. 19). However, the foreign species 
still permeates the entire flow field because Re is moderately low (~10 2 ). 

Now if we increase body size (going from figs. 8-19 and 8-20), the result is that we regain the 
staicture of both a momentum and a thermal boundary layer (i.e., there are large changes in u/u s and 
H/Hj near the wall). The foreign species vanishes at only a third of the distance from the wall to the 
shock in spite of the fact that its concentration at the wall is an order of magnitude larger than that of 
figure 8-19 and the mass addition rate is the same. These are all the effects of going to a larger Reynolds 
number (-10* in table IV), and they underscore the importance of preserving (at least approximately) the 
Reynolds number in scaling mass addition effects. 

Now we specialize the scaling discussion to the ablation case. 


Scaling With Ablation 

Although it is generally not possible (and sometimes not desirable) to simulate or scale all of the 
pertinent parameters in the laboratory, it is nevertheless worth examining the extent to which mass addi- 
tion can be scaled in the presence of conduction, gas and surface radiation, and ablation. 

. Conduction, radiation, and reradiation - The expression fnr the material f w as obtained from 
equations (15) and (13) can be written 
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w 


= j -\[^ 

ha V e(j + 1) 


a w B(U,p s ) R 


PsU 


+ 2.015xl0- 5 (2e) 1/4 n ol9 U 2 A / — e“ b ^ f w) n - 

V p s R 


ae w T?, 


W 


PsU 


(29) 


We confine our attention to one material and assume that h a , a w , e w , and T w are constant and, more- 
over, that e and iJO.19 are almost constant. Then, in order to scale the flow-field concentration, we 
require (from the preceding section), that both f w and m/pooU be fixed; or, alternatively, we require that 
both f w and Re be fixed. From equation (29) these can be fixed if the brackets in equation (29) are 
fixed. Conceivably, one could find a range of flight conditions and nose radii for which the bracket and 
f w within the bracket are fixed. The constancy of the bracket would then constitute a somewhat unap- 
pealing scaling law. It could be specialized to a set of simple laws, namely that 


RB(U,p s ) 

PsU 


U _Kf_f \n 1 

e , and — =■ 

Vp s R p s U 


are individually constant, where these pertain to radiative, convective, and surface reradiative transport, 
respectively. If convective heating is negligible, we can scale by keeping p s U and R fixed between 
model and prototype. 


Note that, in general, we cannot have constant-velocity scaling with this set of simple laws; that 
is, p s and (because b = b(£J)) R must be individually constant so that there is no hope of scaling either 
flow field or model to achieve foreign-species scaling. Even if we neglect either radiative or convective 
heating, we cannot have constant-velocity scaling because of the surface reradiative term. 

Radiation only- If this is the only energy transfer mode, scaling can be accomplished if f w and 
Re are constant, which requires that BR/p s U be constant. Or the last can be replaced by constant 
BR 2 /p s because the Reynolds number is fixed. For scaling at constant velocity, BR/p s must be 
constant. 


Convection Qnly- This situation is especially applicable to ground-based tests with noncharring 
ablation models in arc-heated wind tunnels and ballistic ranges. If this heating mode is assumed to pre- 
vail also for the prototype at actual flight conditi ons (a s was the case for the mass addition calculation of 
ref. 31), the scaling law requires that either UWp s R]e~ b ( _fw ) oij^^/Vj^Je-^-fw) 1 be constant (as 
well as constant Re and thus constant fw)- For constant- velocity scaling, tne requirement is simply that 
•vPs be constant (along_with Re). Since -v Ps is a weak function of altitude for constant velocity, the 
conclusion is that m/pooU is constant and the foreign-species profiles are scaled simply if Re is pre- 
served^ T. N. Canning and G. Chapman of Ames Research Center have advanced the former conclusion 
(nVpooU constant for constant Re) based on phenomenological arguments under the same constraints 
(simulated shape, Reynolds number, velocity, and noncharring material (no reradiation)), and for 
convection only. 


In short, then, ablation scaling (by simple rules) in terms of ratio of mass flux at the wall to that 
in the free-stream and foreign-species concentration profile can best be achieved for very special heating 
conditions. Constant-velocity scaling cannot be achieved for materials which reradiate importantly. 
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Low Reynolds Number Effects Without Mass Addition 

It is common to study mass addition effects on convective heat transfer by use of the effective 
heat of ablation and a heating rate corresponding to no mass addition. For large Reynolds number (no 
shock-layer vorticity) and high speeds (up to 50,000 ft/sec), the convective heating rate without mass 
addition is quite well known both by experiment (refs. 6 and 32-35) and theory refs. 2, 6, 36, and 37). 
Moreover, a number of investigators have studied the regime of Mach numbers up to 8, where Reynolds 
numbers are low enough that shock-layer vorticity affects the convective heating rate for no mass addi- 
tion (refs. 38-44). Van Dyke (refs. 45 and 46) has studied the vorticity effect up to infinite Mach num- 
ber without real gas effects. Hoshizaki (ref. 47) has examined the vorticity effect for the incompressible 
shock layer. Our intent is to examine briefly the external vorticity (or low Reynolds number) effect on 
wall shear stress and convective heat transfer at high speed using real gas properties. 5 

Surface shear stress - The present flow-field analysis is a single-layer analysis in which the equa- 
tions are solved from the body to the shock. A very simple comparison of the single-layer result with 
the no-vorticity, two-layer (boundary layer plus inviscid shock layer) result for surface shear stress can 
be made as follows. 

The single-layer shear stress has been expressed by equation (28). It is simple to show for the 
two-layer, no-vorticity analysis for a cold wall (ref. 2, for example) that 

t»l _(uf/2 (2 e)W VpwMj + 1) f; 

x (no von W in ' ' 

where FJJy is the value at the wall corresponding to the solution of the Blasius equation for no mass 
addition, and is given by reference 48. The ratio of t w /t Wno „ obtained by combining equations (28) 
and (30) is 


% _ 2 lf2 (cpw)^ p 

(Tw)no vort (2 e)^ p" 

The result is shown in figure 8-28 in which the ratio is plotted as a function of Reynolds number of the 
form used by reference 39. As would be expected, shear stress increases over its no-vorticity value as 
Reynolds number decreases. At a given Reynolds number, the effect is enhanced by increasing speed 
and decreasing pressure. 


5 Because of a conflict which arises in the outer boundary conditions for very small Reynolds numbers, the results of this 
section and those of figures 8-18 and 8-19 are considered to be first approximations. The conflict is discussed and evaluated 
in an approximate way in the appendix. 
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Figure 8-28 - Effect of Reynolds number on surface shear stress with no mass addition. 

Convective heat-transfer rate - The convective heat-transfer results for the same examples are 
shown in figure 8-29. For present purposes, q Cno vort was obtained by use of equation (1 1), which is the 
result of Hoshizaki (ref. 6). It can be noted that the convective heating result parallels the shear-stress 
result as would be expected. 


■ 5.7 AND 8 
T s = 1170 AND 1280 K 
e = 0.167 
(REF. 38) 

(REF. 39) 



HICKMAN AND GIEDT 
= 2-6, T s = 300 K 
0.282 (REF. 43) 


(U(ft/$ec), p s (atm)) 
PRESENT RESULTS 



Figure 8-29.- Effect of Reynolds number on convective heat transfer with no mass addition. 
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The theoretical and experimental results of others (as obtained from refs. 39, 44, and 47) are 
shown for comparison. In each separate study the ratio of convective heating with vorticity to that with- 
out is enhanced by an increase in speed or Mach number (M) or (for the case of ref. 47) e at a given 
Reynolds number. Among the various studies, the ratio does not always increase with Mach number and' 
there is some disagreement as to the reason (discussed in refs. 38, 44, and 46). The present results lie 
considerably higher than those of references 39, 41, 43, 45, and 46 as would be expected because of our 
comparatively high speed (or Mach number) and/or low e. 

Our results are close to the viscous layer results of Probstein (ref. 40) and Hoshizaki (ref. 47). 

The results of these three studies are all derived from flow-field analyses which employed the Navier- 
Stokes equations from the body to the shock wave. Both Probstein and Hoshizaki assumed constant 
density and Prandtl number. Our results lie above theirs for a given Reynolds number and e, possibly 
because of our variable Prandtl number corresponding to higher speeds for which ionization occurs, and 
possibly because of compressibility effects near the cold wall. The slopes of our lines are alike and are 
much like those of Hoshizaki at lower Reynolds numbers. The vorticity results show that convective 
heating may be as much _as 60% higher than the no-vorticity value and that the ratio of the two is 
enhanced by increased U or diminished p s (or e) at constant Reynolds number. 

The ratio of convective heat-transfer rate to total free-stream energy flux for these same examples 
is shown as a function of Reynolds number in figure 8-30. The results of references 38 and 39 at Mach 
5.7 and 8 are also shown. In each case the slanted line corresponds to the no-vorticity result. At a given 
Reynolds number, Ch, the ratio of convective heat transfer rate to the free-stream energy flux without 
vorticity, increases with increasing speed. The symbols represent the present result and are attached to 
the appropriate no-vorticity line by a vertical line. Conservation of energy requires that Ch not be 
more than unity. Thus Ch of unity represents the flow energy limit. For the examples represented by 
the symbols, the convective heating was less than half the flow-energy limit at the most. 


FERRI, ZAKKAY, M = 5.7, 
T s - 1170 K (REF. 39) 



Figure 8-30 - Effect of Reynolds number on heat-transfer coefficient for no mass addition. 

Finally, flow-field solutions for some of the nonablating small-body points of figures 8-28 
through 8-30 are shown in figures 8-31 through 8-33. The trend toward increased vorticity near the wall 
relative to that near the shock, and toward a more nearly isoenergetic flow field behind the shock for 
higher pressure can be seen from comparison of figures 8-31 and 8-32. 
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Figure 8-31- Nonisoenergetic flow field with shock-layer vorticity for no mass addition 
(U = 41,000 ft/sec, p s = 1 atm, R = 0.0104 ft). 



Figure 8-32.- Nonisoenergetic flow field with shock-layer vorticity for no mass addition 
(U = 41,000 ft/sec, p s = 10 atm, R = 0.01 ft). 
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Figure 8-33. Nonisoenergetic flow field with shock- layer vorticity for no mass addition 

(U = 50,000 ft/sec, p s = 1 atm, R = 0.02 ft). 


CONCLUDING REMARKS 


The flow equations in the stagnation region of the shock layer of blunt bodies (including mass 
momentum, and energy transport phenomena) have been solved for numerous examples at flight speeds 

up to 50,000 ft/sec in air. The thermodynamic and transport properties of dissociating ionizing air were 
used m the analysis. to 

Many results with mass addition (by transpiration or ablation) were obtained. It was shown that 
convective heating was more strongly affected by mass addition than was radiative heating for mass 
addition rates up to half the free- stream mass flux (excluding effects of radiation from ablation prod- 
. Convective heating results were correlated by a simple relation which shows that mass addition 
diminishes convective heating exponentially, where the argument of the exponential is a simple function 
of flight speed and Blasius-type wall stream function to the 3/2 power, (-f w ) 3/2 . Results with mass 
addition for body nose radii between 0.01 and 5.0 ft, flight speeds from 30,000 and 50,000 ft/sec wall 
temperatures from 1500 to 3000 K, shock-layer pressure levels from 0.1 to 10.0 atm, and surface mass 
addition rates up to half the free-stream mass flux were correlated by the simple expression. 

Previous correlation formulas obtained from subsatellite speed results do not correlate the present 
higher speed convective heating results corresponding to high mass addition rates. 

The results with mass addition were used to study ablation at hypervelocity for which convec- 
tion, gaseous radiation, and surface reradiation were taken into account. At specified flight conditions 
(ranging in speed between 30,000 and 50,000 ft/sec and between 0.1- and 10.0-atm shock-layer pressure 
level) for a given ablator, the nose radius which minimizes total heating rate was determined. For this 
optimum ' nose radius, the following results were obtained. 


TT"Tn 
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1. The ablation rate in terms of the stream function at the wall depends only on the convective 
heating properties of the gas and is independent of the gaseous radiation properties if reradiation from 
the surface itself is negligible. 

2. The ablation rate in terms of mass flux at the surface is not more than one-fourth the free- 
stream mass flux for the flight regime cited above and is not more than one-half the free-stream mass 
flux at a speed of 70,000 ft/sec. 

3. The influence of mass addition on standoff distance is moderate, the distance being enhanced 
by not more than 50% over its no-ablation value. 

4. The gaseous radiation flux incident on the surface is larger than the conduction flux by a 
factor of 2 to 9. 

5. If the ablated vapors radiate like air, mass addition enhances the incident radiant flux at the 
surface by less than 20% over the no-ablation value. 

6. If the ablated vapors radiate 10 times as strongly as air, mass addition enhances the incident 
radiant flux at the surface by less than 50% over the airlike value. 

7. The surface viscous shear stress is low, generally not more than 5% of the surface pressure. 

In order to scale mass addition effects in terms of dimensionless concentration profiles of the 
foreign species, it is necessary to match both Reynolds number and either the stream function at the wall 
or the ratio of the mass flux at the wall to that of the free stream. 

Ablation rate and foreign-species concentration profiles can be scaled conveniendy for very spe- 
cial heating conditions only. Constant-velocity scaling cannot be achieved for materials which reradiate 
importantly. 

Finally, at low Reynolds numbers, the shock-layer vorticity enhances both shear stress and con- 
vective heating over the no-vorticity values, the effect being greater at higher speed and lower pressure. 
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APPENDIX 


BOUNDARY CONDITIONS BEHIND BOW SHOCK WAVE FOR VERY SMALL NOSE RADn 

OR REYNOLDS NUMBER 


The outer boundary condition employed for the energy equation is 

H s = U 2 /2 


(Al) 


However, for nose radii small enough that a conductive heat flux exists just behind the shock, the 
boundary condition may be quite different from equation (Al) as follows. 6 

For simplicity, we consider only the normal portion of the shock wave and equate energy flux on 
each side of the shock thus, 


pooU hoo ^ 


U2\ 


, , k dh \ I, v s 

+ 1 — — = -p s vJ h s + -f- 

c p dyL 1 s 2 


k dh 


I V-puy/oo \ *■] lcpdy/ s 

where y and v are positive outward from the body. We employ the strong shock approximation 


(A 2) 


hoo ^ ' 


ui 


(A3) 


and assume 


From mass continuity, 


By use of equations (A2)-(A5) 



-v s = eU 


(A4) 


(A5) 


H = U^ 1 /j<_dM 

S ~ 2 Poo u|c p dy| s 


(A 6) 


So we see that H s is always less than 11^/2 by an amount that cannot be determined until the solution 
of the flow field behind the shock has been obtained and (dh/dy) s is known. The problem at once 
becomes one of iterating not only on the simultaneous differential equations, but on their boundary con- 
ditions as well. 


6 A similar difficulty exists with the outer boundary conditions on the momentum equation because of vorticity behind the 
shock, as noted in reference 47. One way out of these difficulties is to integrate through the shock wave and use the free 
stream as boundary conditions, as presented in reference 49. Even at that, many uncertainties in the shock wave (e.g., 
thermal and chemical relaxation phenomena) remain at these high velocities. 
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We can estimate how much H s differs from U /2 by assuming 


dh _ h$ — h w ~ h s 
dy ~ 5 J 


Then H s becomes 


^U 2 (1 +e 2 y) 
s 2 (1 +y) 


(A7) 


(A 8) 


where 


1 

I— 

Pr s 

P ooU(5/R)R 

i c P< 

s Re(8/R) 


(A9) 


Since 5/R is not a strong function of R,U, or p s in the examples considered, y is inversely propor- 
tional to the Reynolds number and directly proportional to the Prandtl number behind the shock wave. If 
y« 1, we have the usual boundary condition H s =U 2 /2. However, as Re gets smaller, that is. 


if 

y - 1 

Y » 1 

e 2 Y » 1 


then 
H - 02 

Hs ~T 

Hc»S(- + e 2 


H s -?-e 2 


In the limit the approximation tells us that H s 


— 9 

is much less than U /2. 


(A10) 


For the moment, we adopt the point of view that H s was specified, and that the result of the 
corresponding solution might apply directly to a hig her flight speecLT hat is, the true flight speed is 
larger than the assumed flight speed by the factor 7(1+ Y)/(l + e 2 y). For the results in question, 
shown in figure 8-29, the factor varies from 1.0 to 1.2. 

However, we cannot make direct application to the higher speed because we then_violate another 
boundary condition in which the smaller value of U was employed; namely, p s v s = -p<x,U. For these 
reasons, the present low Reynolds number results are considered to be first approximations. 
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CHAPTER 9 


VERY SEVERE FLIGHT ENVIRONMENTS, THE JUPITER PROBE 


Thus far we have learned something about hypervelocity flight in Earth's atmosphere. We have 
considered stagnation-region heating in considerable detail. The flow has been in chemical and thermo- 
dynamic equilibrium in the illustrations presented. We can now say that the apparent sterile formalism 
of kinetic theory and radiative transfer have been applied with some fidelity to real gas flow. 

Flight in other planetary atmospheres has been studied: Mars, including return to Earth (ref. 1); 
Venus, Uranus, and Saturn (and its largest moon, Titan (ref. 2)); and the most severe entry environment 
that we have studied in detail — Jupiter (ref. 3). The last of these is the subject of this chapter. 

Both convective and radiative heating rates for various hypervelocity flights (circles) and some 
facilities that were used to simulate those flight environments (squares) are shown in figure 9- 1 of the 
AIAA Survey Paper presented in this chapter. Thus we see Apollo, Pioneer Venus, Space Shuttle, and 
the Jupiter Probe. The Jupiter entry is two orders of magnitude more severe in terms of heating than the 
Apollo flight. Correspondingly, the Giant Planet Probe Facility (GPPF) at NASA Ames Research 
Center is shown, which partially simulated the Jupiter Probe entry for a small body. 

The Jupiter Probe is part of what is called the Galileo Mission. Although the probe has been 
designed and built, its launch has been delayed for lack of a launch vehicle. Nevertheless, it is instruc- 
tive to examine the phenomenological events involved in the design of the probe destined to enter an 
environment that is an inferno. 

It may be remarked at the outset that the heat-shield material for the probe was selected to be 
carbon phenolic. Although that choice was partially attributed to the writer, I preferred a heat shield that 
back-scattered radiation in the depths of the material (even though the surface was ablating) in this 
intense radiative environment. Indeed, that technology was well advanced, but it was not given serious 
consideration. 

A descriptive survey of the effort constitutes the balance of this chapter, which was originally 
presented as AIAA Survey Paper 81-1068. References starting with reference 4 are those of that paper. 
The work is that of many people at NASA Ames Research Center, Langley Research Center, The Jet 
Propulsion Laboratory, the General Electric Co., the Acurex Corporation, and other organizations and 
people. It was a major technology effort. There is no reason to highlight the features, uncertainties, and 
difficulties of the supporting research and technology effort in advance, so we move directly into the text 
of the paper. 
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(Text of the) 

SURVEY OF THE SUPPORTING RESEARCH AND TECHNOLOGY FOR THE THERMAL 

PROTECTION OF THE GALILEO PROBE' 


Abstract 

The Galileo Probe, which is scheduled to be launched in 1985 and to enter the hydrogen-helium 
atmosphere of Jupiter up to 1475 days later, presents thermal protection problems that are far more diffi- 
cult than those experienced in previous planetary entry missions. The high entry speed of the Probe will 
cause forebody heating rates orders of magnitude greater than those encountered in the Apollo and 
Pioneer Venus missions, severe afterbody heating from base-flow radiation, and thermo-chemical abla- 
tion rates for carbon phenolic that rival the free-stream mass flux. This paper presents a comprehensive 
survey of the experimental work and computational research that provide technological support for the 
Probe's heat-shield design effort. The survey includes atmospheric modeling; both approximate and 
first-principle computations of flow fields and heat-shield material response; base heating; turbulence 
modeling; new computational techniques; experimental heating and materials studies; code validation 
efforts; and a set of "consensus" first-principle flow-field solutions through the entry maneuver, with 
predictions of the corresponding thermal-protection requirements. 


Introduction 

The Galileo Probe is scheduled to be launched sometime in 1985 and to enter Jupiter's atmo- 
sphere as many as 1475 days later. After entry, a parachute will be deployed, both the forebody and 
afterbody heat shields will be jettisoned, and scientific measurements will be made as the payload 
descends through the atmosphere. Because this entry severely tests our thermal-protection design capa- 
bility, a significant research base — both experimental and analytical — has been formed in support of the 
probe heat-shield design effort. Since about July 1975 (ref. 4) the supporting research and technology 
base has involved Ames Research Center, Langley Research Center, and several contractors. 

The high-speed entry (up to 48.2 km/sec relative to the hydrogen-helium atmosphere of Jupiter) 
is expected to lead to probe forebody heating rates that are hundreds of times greater than those of 
Apollo and 10 times the Pioneer Venus rates (fig. 9-1). The corresponding thermo-chemical ablation 
rate of the carbon phenolic heat shield is expected to rival the free-stream mass flux because of the 
intense radiative and convective heating from the hot turbulent gas cap that will envelop the probe. 
Moreover, for the cool dense model of the Jovian atmosphere, it is likely that radiative heating rates near 
the stagnation point will be so severe that the heat-shield material will undergo spallation as well as 
thermo-chemical ablation. It is calculated that the severe environment will cause the probe forebody to 
lose about one-third of its mass during the entry heating pulse. Not only is forebody heating a severe 
problem for the Jovian probe, but the afterbody will incur significant radiative heating from the wake 
region of the flow. As a result, it is necessary that the flow field and the heating rates about the entire 
probe, including the wake, be understood as thoroughly as possible in order for the design of the thermal 
protection systems to be as effective as practical. 


'Written by J. T. Howe, W. C. Pitts, and J. H. Lundell, Ames Research Center, as AIAA Survey Paper 81-1068. 
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Figure 9-1- Heating environments. 


The flow environment and material response about the probe is beyond our previous capability to 
simulate experimentally, or to compute realistically. In recent years, new experimental facilities have 
been devised, constructed, and used to simulate some but not all of the features of the entry environ- 
ment; indeed they cannot simulate the flow field in all respects. The Ames Giant Planet Probe Facility 
(GPPF) can produce heating rates that simulate the heating levels on the flank of the probe, and the 
Ames Gas Dynamic Laser Facility can simulate heating rates to the stagnation region of the probe and 
produce heating levels that cause spallation of the heat-shield material. These facilities will be discussed 
subsequently. Both the forebody and afterbody flow fields and radiative transfer have been studied in 
shock-tube facilities, and by firing free-flight models into a noble gas mixture to simulate Mach number 
and Reynolds number in order to observe forebody and wake flow-field structure. Moreover, various 
computational codes have been developed to predict the effects of a great many coupled physical phe- 
nomena in the flow field and in the heat-shield material. 


A number of additional problems have been addressed to assess their effects on the probe heat 
shield. These include the composition and thermal structure of the atmosphere, probe shape, and tar- 
geting (initial entry) conditions. But some phenomena (such as spallation) that cannot readily be 
assessed by computational physical analysis must be investigated experimentally. The final objective of 
the comprehensive supporting research and technology effort is to provide an acceptable set of first- 
principle computational codes to be used as a standard for the more approximate design codes of others. 
To enhance confidence in the results of these first-principle, or benchmark, codes it is necessary to vali- 
date the codes experimentally, where possible, and to provide experimental corrections to account for 
phenomena that cannot be modeled in the codes. 


Targeting, Atmospheric Modeling, and Trajectory Computations 

Targeting for the entry probe has been at an altitude of 450 km at various latitudes near the 
Jovian equator. Effects of variations in the azimuthal angle (near equatorial), of the entry angle of the 
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probe, with respect to the planet, and of probe mass and configuration on thermal-protection require- 
ments have been examined. That target altitude was fixed early in the study, because it corresponded to 
what was thought to be the "top" of the model atmosphere, according to reference 5. Subsequently, 
however, there were a number of other model atmospheres which differed both in thermal structure and 
in relative abundance of hydrogen and helium. The thermal structures of these different nominal atmo- 
spheres are shown in figure 9-2, which was adapted from an Acurex Corporation report. The original 
atmospheric model, which corresponded to that reported in reference 5, was composed of 85% hydrogen 
and 15% helium by volume. That was followed by the atmosphere shown by the solid line in figure 9-2; 
it was referred to as the Hunten model. That model had the thermal structure given by reference 6, but a 
composition of 89% hydrogen and 1 1% helium in accord with reference 7. This was followed by a 
nominal atmosphere model developed by G. Orton of the Jet Propulsion Laboratory which has the ther- 
mal structure shown in figure 9-2, and a composition of 89% hydrogen and 1 1% helium, as noted above. 



Figure 9-2.- Thermal structure for nominal model atmospheres of Jupiter. 

The atmosphere was assumed to be in hydrostatic equilibrium. Hence, an ordinary differential 
equation could be solved for each model atmosphere composition and thermal structure to compute 
density-altitude profiles for each model. These computations were performed by M. J. Green, W. C. 
Davy, and D. Kirk of Ames Research Center. (Mr. Kirk and his co-workers also contributed to the 
Galileo SRT effort in the area of flight mechanics, which is not included in the scope of this paper.) A 
comparison of two sets of these atmospheres is shown in figure 9-3, in which the corresponding atmo- 
spheric composition is also shown. Thus for each atmosphere defined by a thermal structure there is a 
warm, nominal, and cool version, depending on the composition. The Orton cool and nominal model 
atmospheres are seen to be significantly higher than those of reference 5. (Subsequently, for the heavy 
(310-kg) probe study, a revised version of the Orton atmosphere was used.) Nevertheless, targeting was 
still specified at 450 km, and trajectories were calculated such that the probe would arrive at that altitude 
with the specified entry angle. Trajectories were calculated to include the effects of latitude, azimuthal 
angle, entry angle, and ballistic coefficient for each candidate probe shape, and to accommodate the fact 
that Jupiter is a rotating oblate spheroid. 


ATMOSPHERIC 

MODEL 

COMPOSITION H 2 /H e (VOL) | 

SP 8069 

ORTON 

WARM 

96/4 

100/0 

NOMINAL 

85/15 

89/11 

COOL 

70/30 

78/22 



Figure 9-3.- Atmospheric models of Jupiter. 

Entry was usually at quite shallow entry angles (5 to 11° inertial entry angle or 6.25 to 13.72° 
relative to the rotating atmosphere (ref. 8)) with flight in the posigrade direction to take advantage of the 
velocity of the rotating atmosphere (about 12 km/sec). In this way the velocity of the probe relative to 
the atmosphere was significandy diminished (from about 60 km/sec inertial to 48 km/sec relative) in 
order to reduce the heating rates and the weight of the thermal-protection system. Probe shapes were 
studied that varied from 35°-half-angle hyperboloid forebodies (ref. 9) to 60°-half-angle sphere cones 
(ref. 10) with a bluntness ratio of about 1/2. The baseline shape was a 45°-half-angle sphere cone with a 
nose radius of 0.31 m and a base radius of 0.62 m. The shape was modified slightly late in the study. 
The probe weight was increased during the course of the study as greater demands were placed upon the 
probe. Much of the early work was for 242-kg probe; later work was for 290-kg and 310-kg probes, as 
will be shown. 


Flow About the Forebody and Heat-Shield Response 

During the heating pulse the flow is determined over the forebody at a number of times through 
the trajectory. To provide a reference set of benchmark solutions, the physical phenomena illustrated in 
figure 9-4 (adapted from ref. 11) must be included in the analysis. Thus the flow equations include the 
transport of mass, momentum, energy, and species. The gases are hot and reacting, and emit thermal 
radiation which is incident on the heat shield. The radiative and convective heating cause the heat-shield 
surface to ablate at a rate that is sometimes comparable to the free-stream mass flux. The ablated gases 
mix with and react with the atmospheric gases, and enter into the radiative transfer in the shock layer. 
The gas mixture makes a transition from laminar to turbulent flow as it is swept around the body. To 
obtain the benchmark flow-field solutions, a quasi-steady-state ablation boundary condition is applied 
such that the ablation rate is simply related to the net heating to the wall. Subsequently, the heat shield 
is designed by applying the heating rates obtained from the flow field to a transient material-response 
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analysis of a material of finite thickness; the latter material is receding and losing mass by thermo- 
chemical ablation of the surface char, and by the outgassing of the pyrolysis gases caused by the internal 
conduction of heat. The thickness of the heat shield is determined by a specified allowable temperature 
at its back surface. 
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Figure 9-4.- Heat-shield and flow phenomenology. 

A number of flow-field studies were performed early in the research effort, using either inviscid 
codes or various codes that employ correlations and physical estimates to assess the existing state-of- 
entty environment- prediction methods (ref. 1 0); and the extent and importance of atmospheric structures 
on radiative heating (ref. 8), radiative heating distributions (ref. 12), problems concerned with an atmo- 
spheric reconstruction experiment (ref. 13), and configuration effects on radiative heating (ref. 14). 

Moreover, an extensive effort was made to acquire and refine the radiative properties of carbo- 
naceous ablation gases both experimentally and theoretically (refs. 15-25). The detailed effect of these 
carbonaceous species is illustrated in figure 9-6, which was adapted from reference 26. Figure 9-5(a) 
shows the spectral radiative flux incident on the stagnation region of the entry probe near peak heating, 
and figure 9-5(b) shows the spectral flux incident on the surface. The difference between the two 
spectra represents absorption by spectral bands of the ablation gases, as shown in the figure. The 
spectral-property studies were used to provide input data primarily to the detailed benchmark flow-field 
codes, as mentioned previously. 

Three primary benchwork codes, which were developed over an extended period of time, solve a 
coupled set of partial differential equations for the mass, momentum, energy, and elemental species 
transport between the forebody and shock wave. This set of equations was complicated by a radiative 
flux divergence term, which is a double integral over space and wavelength. The three codes— RASLE, 
HYVIS, and COLTS— were developed primarily by W. E. Nicolet, J. N. Moss, and A. J. Kumar, 
respectively. Briefly, RASLE (ref. 11) solves the equations simultaneously by a Newton-Raphson 
numerical technique; HYVIS solves the equations sequentially by iteration; and COLTS is a time- 
dependent version of HYVIS. The development of these codes cannot be described in detail here; in 
general, however, the laminar stagnation region shock layer was first solved by RASLE (ref. 27) and 
was used to provide ablation blocking corrections for more approximate codes used in a parametric 
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Figure 9-5.- Spectral flux attenuation by ablation species: Shock temperature = 16,090 K, pressure = 

5.1 atm, 89% H2, 1 1% He atmosphere, carbon phenolic heat shields, (a) Spectral flux incident on 
ablation area, (b) Spectral flux incident on wall. 

study (ref. 28), in which comparisons with other parametric studies were made. These results were used 
in reference 29 to assess heat-shield requirements and were compared with the pioneering work of 
Tauber (ref. 30). 

The HYVIS code was a pioneering code in several respects. The code produced laminar solu- 
tions about hyperboloids of revolution. The solutions were over the entire carbon phenolic forebody 
heat shield, and included coupled ablation as a quasi-steady-state boundary layer (ref. 31). Subse- 
quently, the code also produced turbulent-flow solutions about both sphere-cones and hyperboloids with 
coupled ablation for carbon phenolic heat shields (ref. 32). 

The RASLE code also computed turbulent flow over a sphere-cone heat shield at various points 
in an entry trajectory (ref. 33), and the results were compared with those of reference 29. The turbulence 
models for the two codes differed somewhat. HYVIS used a Prandtl mixing length near the wall and a 
Clauser-Klebanoff outer eddy viscosity; RASLE used an ordinary differential equation for the mixing 
length near the wall, selected partly on the basis of an experiment of reference 34 for a highly cooled 
wall, and a wake-like mixing length differential equation away from the wall. Surface radiative condi- 
tions for the two codes also differed, so that there were significant differences in the integrated mass loss 
obtained from results of the two codes through the trajectory, as shown in figure 9-6. The inertial entry 
angle was -9° in the nominal atmosphere. Results from the COLTS code (ref. 35) are also shown; in 
general, they agree with the results of HYVIS. It may be noted that HYVIS was able to obtain solutions 
at the beginning of the heat pulse and after peak heating, but not in between. RASLE solutions were 
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obtained throughout the heating pulse. The COLTS code seems to have overcome some instabilities 
experienced by HYVIS and to have obtained two additional solutions on the upward part of the heating 
pulse. 


242-kg PROBE 
X RASLE, a = e = 0.9 



TIME FROM 1000 km ALTITUDE, sec, 
ORTON NOMINAL ATMOSPHERE 


Figure 9-6 - Probe mass-loss-rate comparison for quasi-steady ablation of carbon phenolic. 

A very detailed comparison of the results of these benchmark codes was performed at two times 
in the trajectory (ref. 36): 1 10.2 sec and 1 1 1.3 sec. Differences in physics, methods, boundary condi- 
tions, and "convergence" criteria were presented, and the effects were assessed. An attempt was made to 
reconcile results of the codes by making adjustments to correct for some differences. The work reported 
in reference 36 is a remarkable effort to sort out many coupled interrelated effects to determine those 
which are most significant. Global mass and energy conservation checks were made. Among three 
major contributors to differences in the results, the turbulence models are considered to be the most 
important, because of their effect on radiative transfer to the wall. Although the various turbulence 
models employed do not affect state profiles near the wall significantly, profiles remote from the wall 
are affected. Thus, temperature profiles away from the wall differ because of the turbulence model. 
Correspondingly, species profiles differ as well. Thus the species that emit and absorb radiation are dis- 
tributed differently in differing thermal environments by the different turbulence models; the radiation 
emitted differs correspondingly. Because radiation is the most important phenomenon causing heat- 
shield ablation, the effect of turbulence modeling on radiative transfer is of primary importance. Results 
of this important paper (ref. 36) have affected subsequent benchmark solutions and contributed to a 
'consensus" set of benchmark solutions shown in part subsequently in this paper. 

Turbulence modeling of a very hot, chemically reacting flow over a relative cool surface that is 
ablating at a rate comparable to the free-stream mass flux is very poorly understood. A comparison of 
several candidate turbulent models (including those cited above) has been made (ref. 37) in the context 
of a given flow-field code to isolate the effect of the model alone, and to rank the various models 
according to severity. It is clear that this is an important problem. Experiments must be devised that 
will make it possible to construct realistic turbulence models and to provide input data for those models. 

As noted previously, the flow-field results shown in figure 9-6 are obtained by use of a quasi- 
steady-state ablation condition. The heat shield designed by the CMA code — using the net heating out- 
put from the RASLE code as input to the in-depth materials-response code — is shown in figure 9-7. The 
CMA code computes the transient thermal response of a charring ablating material of finite variable 
thickness, and includes the effects of temperature-dependent properties and in-depth pyrolysis. The 
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criteria that "design" the heat shield are the maximum allowable backface temperature of 700 K and the 
safety margins (shown in fig. 9-7). These margins are percentages of the ablated material that are added 
for safety purposes. The figure shows the initial and final shapes of this probe designed for entry into 
the nominal atmosphere. The mass loss (47 kg) is less than that shown in figure 9-6 (53 kg). The 
difference is affected by the heat storage in the transient, in-depth solution in only a minor way; this 
point will be discussed subsequently. When the heat shield with the safety margins defined by fig- 
ure 9-7 is allowed to enter the cool dense atmosphere at an inertial entry angle of-1 1°, an additional 
thermal protection mass (22%) is required for survival, with no safety margin (ref. 38). These prelimi- 
nary results are for a 242-kg probe and are illustrative only. Subsequently, we shall cite research that 
reconciles some of the differences among the benchmark codes and present results for a more recent 
heavier probe (310 kg). 


242-kg PROBE 
NOMINAL ATMOSPHERE 
50% MARGIN AT STAGNATION POINT 
34% MARGIN ON CONE 

a = e ~ 0.9 

MASSOF SHIELD 73.4 kg 
MASS LOSS 47 kg 


ORIGINAL SURFACE 
ABLATED SURFACE 



BACKFACE TEMPERATURE 
CRITERION 


HEAT SHIELD REMAINING 
AFTER ENTRY 


Figure 9-7.- Transient heat-shield response: carbon phenolic. 

The codes that produced the results of figures 9-6 and 9-7 are based on first principles and are 
considered to be benchmark flow-field and material response codes. They might even be labeled "first- 
generation" benchmark codes, because new computational techniques are being devised and applied to 
this problem, as will be shown subsequently. However, even these first-generation benchmark codes are 
very costly to operate and are difficult and time-consuming to run — especially for broad parametric 
studies. Thus for engineering purposes, more flexible, more approximate codes are useful for economic 
reasons. Zoby et al. (ref. 39) have developed an approximate code that simulates inviscid radiating 
flow-field analysis by the use of analytic shock shapes and a 58-step radiative transfer model. Moreover 
a code has been developed to predict laminar and turbulent convective heating of reactive and nonreac- 
tive gases about blunt reentry configurations for hypersonic flight (ref. 40). These approximate codes 
produce results that are in good agreement with more complicated codes, and are used extensively for 
parametric studies. 


Flow about the Afterbody and Base Heating 

Although the convective heating of the base regions is small (ref. 41), radiative heating is not. 
The significance of the base radiative heating was noted by Stephenson (ref. 42) in 1965 by measure- 
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merits made of radiation from the wake of an ablating blunt body launched into air, the measurements 
were made as the body flew past the radiation sensors. 

Park (ref. 43) developed a base flow-field model and derived the base radiative heating from the 
model. Observations from both ablating and nonablating free-flight tests were used to devise the flow- 
field model shown in figure 9-8. The main features of the flow field were an expansion about the shoul- 
der; a recirculation region behind the body; and a wake which recompressed and formed a neck about 
1.5 body diameters aft of the frustum, with a neck diameter between one-half and two-thirds of the body 
diameter. Method-of-characteristics solutions about the shoulder were matched to the recirculation 
region by a choice of base pressure that gave the observed turning angle. Physical reasoning led to a 
theoretical model of the recirculation, recompression, and neck-region flows. In the ablation case of 
interest, chemical and radiative modeling predicted that the most severe radiation would occur immedi- 
ately behind the frustum, because of radiation from hydrogen and carbon atoms that were not in chemi- 
cal equilibrium in the expansion region. At the base stagnation point, radiation is received from both the 
neck and recompression region, and the recirculation region — as shown in figure 9-9. These regions are 
probably in chemical equilibrium, and detailed spectral radiation computations from the hydrogen- 
carbon mixture produce the upper line corresponding to the ablation case. (It is interesting that the non- 
ablating case gives almost comparable results — probably because the temperature is higher.) It is 
believed that the high base pressure obtained by the model is the most significant physical property that 
leads to the high base radiative heating. The effect of base pressure was examined and the radiative 
heating was found to vary as the square of the ratio of base-to-forward-stagnation-point pressure. 

SHOULDER- 





Figure 9-8.— Schematic of base-region flow field, (a) Shoulder region, (b) Velocity profiles. 

(c) Enthalpy profiles. 




ANGLE d, rad 

Figure 9-9.- Radiation at base stagnation point. 

Experiments performed by instrumented nonablating models in shock tubes (ref. 44) showed that 
radiative heat fluxes are indeed significant, and the results were generally in agreement with those of 
reference 43. Comparison of the experimental and theoretical results was made with the early work 
reported in reference 42; agreement was generally good. It was determined that the ratio of base-to- 
front-stagnation radiative heat fluxes and pressure ratios for a Jovian probe are significant, and are in 
agreement with the theory of reference 43. 

A summary of the work of references 43 and 44 is presented in reference 45. The work 
described in references 43-45 is for laminar flow. Work in progress in free-flight facilities addresses the 
turbulent-base-flow problem. 


Second-Generation Flow-Field Codes 

Although benchmark flow-field codes were being used extensively to produce forebody solutions 
with ever-increasing complexity and physical detail, it was apparent that the method was being strained 
to the utmost, and accuracy and convergence problems were becoming more severe (ref. 36). Moreover, 
afterbody solutions were becoming more important, and they were intimately related to forebody events 
(ref. 43). Thus efforts were initiated to develop a new method, one in which new coordinate systems 
and advanced numerical techniques would be used. A new code, CAGI (ref. 46), that was tailored to 
capture the major flow-field features of both the forebody and the afterbody, was initiated; it is presently 
in an advanced state of development. The code solves the unsteady compressible Navier-Stokes equa- 
tions for two-dimensional axisymmetric three-dimensional flow. A finite-volume formulation was 
devised, rather than the usual finite-difference approximations of the differential equations. The finite- 
volume formulation maintains global conservation of mass, momentum, and energy over the computa- 
tional volume in accord with specified fluxes at the volume boundaries. A mesh is constructed over the 
entire probe configuration; the mesh is flow-aligned over four principal regions — forebody shock-layer, 
base recirculation, outer inviscid wake, and inner inviscid wake (ref. 47). There is a singular coordinate 
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mesh topology in the recompression region. The computational code employed is an outgrowth of the 
basic factored implicit algorithm of references 48-52. 

The first results of the code were presented in reference 46, where real-gas properties were mod- 
eled rigorously by a variable "gamma," which reduces to the ratio of specific heats for the case of a per- 
fect gas. A "homogeneous sweep" algorithm was used in that study instead of the full-factored implicit 
algorithm. Results were presented in reference 46 for the case of flight at Mach 50; some results for 
mass addition along the forebody surface were included. Of particular interest is the flow-field structure 
about the entire bluff body and the corresponding computational mesh shown in figure 9-10 (actually 
adapted from ref. 47). The computed pressure distribution is shown in figure 9-11, and a comparison of 
the computed velocity vector field with an experimental shadowgraph by C. Park is shown in 
figure 9-12. 
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Figure 9-10 - Probe flow-field features and the computational grid. 



176 




3.5 X 10 3 



Figure 9-11- Computed pressure distributions. 
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Figure 9-12.- Computed velocity field compared with free-flight shadowgraph. 
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A chemically reactive gas module of the code was introduced (ref. 53) and solutions were 
obtained. The pressure distribution, and the temperature, enthalpy, compressibility, and species profiles 
were computed and compared with results of reference 9. There was generally good agreement between 
the results of the two studies. 


Experimental Program 

In addition to the analytical program, there is also an extensive experimental program on both the 
forebody and afterbody heat-shield materials. The objectives of the program are (1) to provide ground- 
based facilities that simulate the expected Jovian entry heating conditions as closely as possible and 
(2) to evaluate candidate materials in those facilities. The magnitude of the simulation problem is 
illustrated in figure 9-13, which shows the heating distribution over the forebody, at peak heating, for 
one of the early study versions of the entry probe. Note that the nonablating peak heating rate is greater 
than 40 kW/cm 2 at the stagnation point, but that it decreases to 8 kW/cm 2 at the rear end of the conical 
frustum. The convective heating is much more uniform over the forebody, decreasing from about 
12 kW/cm 2 at the stagnation point to 9 kW/cm 2 on the frustum. 


GALILEO AT PEAK HEATING FACILITIES 



0 .5 1.0 1.5 2.0 2.5 ARC JET LASER 


S/R n 4 cm diam 1 X 1 cm 

Figure 9-13.- Galileo peak-heating facilities. 

The heating rate on the conical frustum, where 85% of the forebody heat-shield weight is 
located, is well simulated by the GPPF; the GPPF is illustrated in figure 9-14 and described in detail in 
reference 54. In the arc heater, an equimolar mixture of hydrogen and helium is injected uniformly 
along the 4.32-m-long constrictor and heated to about 15,000 K by an electric arc, which terminates on 
external graphite: electrodes. For the Initial heat-shield test program, the facility was operated under the 
following conditions: power = 57.7 MW, current = 5.33 kA, flow rate = 0.12 kg/sec, bulk enthalpy = 
185 MJ/kg, and model stagnation pressure = 2.2 atm. The combined convective and radiative heating 
rate measured by a 4-cm-diam. thin-shell tungsten calorimeter, with a radius of curvature of 20 cm, was 
14.6 kW/cm 2 . A radiometer of similar geometry was used to measure the radiative heating component, 
and the measured value was 6.0 kW/cm 2 . By combining these two numbers and correcting for the 
reflectivity of the tungsten-calorimeter cap, the applied convective heating rate was determined to be 
12.0 kW/cm 2 . These heating rates pertain to the model test position, which is 7 cm downstream of the 
7-cm-diam. nozzle exit. The facility is equipped with four swing arms, so that a calorimeter, a pressure 
probe, and two ablation models can be inserted into the free jet stream during a given run. The ablation 
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models are 4-cm-diam. flat-faced cylinders consisting of a concentrically wound carbon phenolic shroud 
and a 2-cm-diam. core on which recession and weight measurements are made. 

ATMOSPHERIC 

NOZZLE, M = 1.7 VENT 



SUPPLY 


Figure 9-14 - Giant Planet Probe Facility. 

Heat-shield materials studied to date include six versions of the carbon phenolic forebody 
materials, carbon-carbon composites, graphites, hyperpure silica reflective materials, and the nylon phe- 
nolic afterbody material. Only the work on one version of carbon phenolic and graphite will be dis- 
cussed. In all the test programs, ATJS graphite is used as a control material against which the other 
materials are compared. It also serves as a check on the arc-jet repeatability. Results of some of the 
earliest tests on carbon phenolic and graphite are compared with theoretical calculations in the next sec- 
tion. In later work, all materials were tested at four different exposure times so that the steady-state 
ablation rates could be established. A comparison of the steady-state rates for chop-molded carbon 
phenolic (the nose-cap material) with those of the control material is shown in figure 9-15. For ATJS 
graphite, the steady-state mass-loss rate is 0.41 g/sec and the recession rate is 0.070 cm/sec. Corre- 
sponding values for carbon phenolic are 0.528 g/sec and 0.104 cm/sec. Thus, the carbon phenolic mate- 
rial has a mass-loss rate that is 28.8% higher than that of ATJS, and the recession rate is 48.6% higher 
than that of ATJS. The greater disparity in the recession rate is accounted for, in part, by the fact that the 
density of carbon phenolic is 1.45 g/cm 3 compared with 1.83 g/cm 3 for the graphite. 

Of course, ablation performance is not the sole criterion for a heat-shield material. In addition to 
accommodating entry heating by ablation, the heat shield must protect the bond line between the shield 
and the structure from excessive temperatures. Thus, in addition to its ablative properties, the insulative 
properties of a heat-shield material are important. In flight, the superior insulation properties of carbon 
phenolic will tend to compensate for its somewhat poorer ablation performance. 

Since the GPPF is not capable of producing radiative heating rates comparable to the peak values 
expected in flight, a gasdynamic laser (GDL) must be used to evaluate materials under these severe 
heating conditions. The facility, which is described in detail in reference 55, is a conventional 
combustion-driven laser which produces radiation at 10.6 pm by burning CO to produce CO 2 as the 
lasing medium. Although the lasing wavelength is significantly different from the UV and visible radia- 
tion expected in flight, the difference is not expected to be important, because the heat-shield materials 
are graybody absorbers. When the internal mirrors are freshly polished and aligned, the laser will pro- 
duce an output power of 33 kW when room-temperature nitrogen is injected along with the CO and up 
to 45 kW when the nitrogen is heated to 1400 K by a gas-fired heat exchanger. Thus, when the output 
beam is focused to a 1.0-cm 2 spot, the peak radiative heating rate expected in flight (fig. 9-14) can be 
simulated. 
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Figure 9-15 - Ablation results at 57.7 MW. (a) Carbon phenolic, (b) ATJS graphite. 


Early in the laser testing (ref. 56), it was discovered that carbon phenolic tends to spall (Fig. 9-16) 
when exposed to severe radiative heating. Although this phenomenon may be unique to laser testing, it 
must be assumed, to be conservative, that it will also occur during peak heating of the Galileo probe. In 
order to account for spallation in the heat-shield design, the phenomenon must be quantified; such tests 
were performed. 



Figure 9-16 — Carbon phenolic spallation. 

The experimental setup for these tests is illustrated in figure 9-17 and described in reference 5~\ 
The output beam of the GDL is collimated into an 8-cm diameter and then transformed into a spatially 
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uniform 1.0-cm by 1.0-cm spot by means of a segmented mirror. This device consists of an eight-by- 
eight array of 1-cm by 1-cm polished molybdenum segments. Although the segments are optically flat, 
they are mounted on a curved, water-cooled substrate which gives the device a 1.0-m focal length. The 
mirror divides the incident beam into many 1-cm by 1-cm beams and superimposes them at the focal 
point. The focal image is then reimaged one-to-one, by a 1.5-m-focal-length mirror, into a canister in 
which the test specimen is mounted. The canister (fig. 9-18) is designed with suitable baffling so that 
the laser beam can enter, but the spalled particles cannot escape. An air jet in front of the specimen 
deflects the laser plume upward and allows the beam to reach the surface without significant absorption 
in the plume. 

SEGMENTED 



O 


COLLIMATED 

BEAM 


Figure 9-17 - Experimental setup for spallation tests. 
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Typical spallation results are shown in figure 9-19. The total mass loss is determined by weigh- 
ing the sample (3 cm square by 5.1 cm long) before and after exposure, and the spallation mass loss is 
determined by collecting and weighing the particles in the bottom of the canister. At least three expo- 
sure times are used so that the steady-state total mass loss and spallation rates can be obtained. Note that 
for the data shown, 18.3% of the mass loss is by spallation at an intensity of 23 kW/cm 2 . By repeating 
the test at a variety of intensities, a correlation of spallation fraction as a function of intensity will be 
developed. The calculated radiation history in flight will then be used in conjunction with the correla- 
tion to determine the total mass loss by spallation over the complete heating pulse and over the entire 
forebody heat shield. 



Figure 9-19 - Spallation results for carbon phenolic. 

Because of the power limitations of the laser, the peak flight radiative rate cannot be simulated 
on a spot larger than 1 cm by 1 cm. Use of such a small spot size results in a hole being bored in the end 
of a square sample 3 cm by 3 cm in size, and may produce results that are unique to the test setup. To 
check for scale effects, future work will involve tests with a square segmented mirror 2 cm on a side, 
which will produce a square spot 2 cm on a side on the end of a square sample 2 cm on a side. Thus the 
test sample will be fully enveloped in the beam in a better simulation of the flight situation. Because of 
the power limitation, these tests will be limited to intensities of 10 kW/cm 2 or less. The results will be 
compared with the 10-kW/cm 2 results with the 1-cm segmented mirror to see whether there are signifi- 
cant differences between the two test techniques. Note that the intensity of 10 kW/cm 2 , which can be 
achieved with the 2-cm mirror, is a good simulation of the maximum radiative rate expected on the 
conical frustum of the entry probe. 
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Code Validation 


It has been noted that the computational studies that have been discussed push the limits of theo- 
retical analysis, and that they were used to calculate material responses to environments that cannot be 
simulated fully in any existing experimental facility. Therefore, it is very difficult to validate experi- 
mentally these theoretical methods. The best efforts to do this to date were reported in references 58 and 
59. We present here a very brief summary of some of the results presented in reference 59. 

The data used for this code validation were obtained using the GPPF (ref. 54), wherein the equi- 
molar mixture of hydrogen and helium gas was injected uniformly along the constricted arc (fig. 9-14), 
and was heated to about 15,000 K, as noted previously. The facility operating conditions for the data 
discussed below were as follows: arc power, 55 MW; arc current, 5330 A; mass flow rate, 0.1 18 kg/sec; 
pressure, 5 atm. 

The models used for the test were 4-cm-diam. flat-faced cylinders. The flat face was 7 cm from 
the nozzle exit, and the exit was 7 cm in diameter. The heat flux to the models was measured using a 
thin-shell calorimeter; the calorimeter was also 4 cm in diameter, but it had a 20-cm radius of curvature 
on the front face. The models and calorimeter were mounted on separate rotating arms so that they 
could be swung alternately to the same position of the arc stream, as described previously. 

Absorbed heat flux data from this test are compared with heat fluxes computed by the RASLE 
code in figure 9-20. The calculated values generally fall within the estimated ±10-MW/m7 error bands 
of the calorimeter data. The RASLE solutions were obtained using a nonablating boundary condition for 
the tungsten surface of the calorimeter. The application of the RASLE code to these experimental con- 
ditions was not direct, primarily because the RASLE code is based on hypersonic approximations and 
the test Mach number was only 1.2. These approximations had to be modified for the subject analysis 
without significantly modifying the structure of the RASLE code. It was found convenient to obtain 
RASLE solutions for an equivalent sphere rather than the cylindrical models. From a parametric study 
(ref. 59), it was found that a 4-cm-diam. cylinder is equivalent, in terms of stagnation-point heat- transfer 
rate, to a spherical model with an 8-cm radius. 

The response of the ablation models to the GPPF arc flow is shown in figure 9-21 for two 
materials. For the computed curves RASLE was used for the flow-field analysis and CMA for the 
material response. The agreement is good for the ATJS graphite, but the predicted recession is twice the 
measured recession for the carbon phenolic. As for the Jupiter entry conditions, the predicted convective 
heating was reduced to a negligible value by ablation products. The reason for the much better agree- 
ment for ATJS graphite than for carbon phenolic is not clear. One reason may be that the material prop- 
erties were better known for the ATJS graphite than for the carbon phenolic. The question of material 
properties for carbon phenolic is under investigation. 
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Figure 9-20.- Comparison of computed and measured absorbed calorimeter heating. 



Figure 9-21.- Comparison of predicted and measured ablation tests in GPPF. (a) ATJS graphite. 

(b) Carbon phenolic. 
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Consensus Benchmark Flow-Field Solutions 


To provide a set of solutions as a standard for calibrating other codes, a matrix of flight condi- 
tions for entry into the nominal Jovian atmosphere was specified, and forebody flow-field solutions were 
obtained by the COLTS and RASLE codes. The probe was a 44.25°-half-angle sphere-cone; it is de- 
scribed in reference 60. Two probe weights were studied: 290 kg and 310 kg. The results of the study 
are shown for the 310-kg probe in terms of mass loss rate at various times in the entry trajectory by the 
symbols in figure 9-22. It may be noted that there is consistently good agreement between COLTS and 
RASLE, except near the peak mass-loss rates at 50.3 and 51.5 sec. Examination of these results by J. N. 
Moss, W. E. Nicolet, A. Balakrishnan, W. C. Davy, M. J. Green, and J. T. Howe failed to resolve the 
differences near peak mass loss. Therefore, a "consensus" curve was agreed on; it is shown by the solid 
line in figure 9-22. The mass loss integrated over the trajectory is 101 kg. If the curve was drawn 
through the RASLE results (the square symbols), the result would be an additional 3 kg (approximate) of 
heat-shield mass. The integrated heat-shield mass loss is presented and discussed in references 60 and 
61 for both the nominal and cool dense atmospheres. 



TIME FROM 450 km ALTITUDE, sec 

Figure 9-22.- 310-kg probe benchmark solution results: nominal Jovian atmosphere, nominal entry 

angle. 

The distribution of mass loss rate about the body at 47, 50.3, and 54.1 sec is shown in fig- 
ure 9-23, which was adopted from references 60 and 62. The mass loss rate distribution is shown 
nondimensionally, where local mass-loss rate was normalized by the product of instantaneous flight 
velocity and atmospheric density. At 47 sec, both the integrated mass loss rate (fig. 9-22) and the mass- 
loss distribution about the body (fig. 9-23) are in good agreement. At 50.3 and 54.1 sec, both figures 
show a generally higher mass-loss rate by RASLE — especially on the conical frustum where most of the 
probe mass resides. 
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Figure 9-23.— Comparison of coupled carbon phenolic mass loss distributions about 310-kg probe: 

nominal atmosphere. 

These consensus solutions had a quasi-steady-state ablation boundary condition. Moreover, by 
prior agreement, the radiative boundary condition was that the surface absorbed all of the flux incident 
upon it — a very severe condition. 

In reference 61, the more severe surface heating-rate history derived from the RASLE flow-field 
results (ref. 60) were input into the CMA code, which is a transient material-response code for a pyro- 
lyzing ablating material of finite thickness wherein material properties are temperature-dependent. An 
allowable bond-line temperature between the aeroshell and the heat shield was specified, and the 
required heat-shield thickness distribution was determined. For the nominal atmosphere, safety margins 
were added to account for uncertainties (such as spallation effects). Thus a heat-shield "design" from a 
benchmark material code was produced (fig. 9-24). Similar results were obtained for the cool dense 
atmosphere (without safety margins). The question of survivability of the nominal benchmark "design" 
in the cool dense entry environment was assessed in reference 61. Results of that paper suggest that for 
a thermochemically ablating heat shield with physical uncertainties that fall within the prescribed margin 
allowances, the 310-kg probe would survive the nominal Jovian atmosphere. Moreover, there appears to 
be a reasonable probability that the probe would survive off-nominal atmospheres without margin. 
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Figure 9-24 - Transient heat-shield response: carbon phenolic. 


Conclusions 

A supporting research and technology effort related to the entry of the Galileo probe into the 
Jupiter atmosphere has been reported. Entry at various initial conditions into numerous model atmo- 
spheres by probes of various size and mass has been studied experimentally and computationally to 
determine both forebody and afterbody heating environments and forebody heat-shield requirements. 
Two kinds of codes were used: (1) codes that use engineering correlations and (2) benchmark codes that 
use basic modeled physical phenomena. Detailed comparisons of the benchmark flow-field codes have 
been performed, and they show that the turbulence modeling of a hot gas about a relatively cool probe 
that is massively ablating has important consequences and needs continued development. Ranking of 
various turbulence models according to severity has been performed. 

For the first time in entry technology experience, both forebody and afterbody heating are severe. 
For that reason, a second-generation benchmark code that computes the flow over the entire probe and 
into the wake is being developed. The code uses very advanced computational fluid dynamic concepts 
and techniques. Initial results have been obtained for flow about the entire Galileo probe configuration, 
and some results have been obtained for flow of a reacting gas. 

Experiments were performed in the GPPF to simulate probe flank heating levels, and with the 
gasdynamic laser to simulate stagnation-region heating for several materials. Spallation was studied and 
partially quantified in the gasdynamic laser. Efforts to validate the computer codes by tests in the GPPF 
were performed for several materials. 

For the "heavy" Galileo probe, a set of consensus benchmark flow-field solutions has been 
obtained. Results that contributed to these solutions have been used in a benchmark material-response 
code to assess the prospects for survivability of the probe during entry into the atmosphere of Jupiter. 
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The overall effort has revealed areas for future entry technology research that will lead to under- 
standing of phenomena that are currently uncertainly known and as a result are carried in the weight 
margin. These include turbulence modeling and mechanical spallation in particular. 
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CHAPTER 10 


HYPERVELOCITY FLIGHT AT LOW DENSITY, NONEQUILIBRIUM FLOW 


In figure 1-1 of chapter 1, near the center of the Figure, a short double horizontal line at about 
50-km altitude was labeled "thermochemical nonequilibrium." It was shown to illustrate that at high 
altitudes where gaseous molecules, atoms, etc., are less dense, collisional processes occur at diminished 
frequency. The practical significance is that the gas cap over a hypervelocity vehicle approaches 
thermodynamic and chemical equilibrium at a finite rate — it does not occur instantaneously. For early 
hypervelocity flight (e.g., Apollo), it was not necessary to understand this nonequilibrium regime in any 
detail- it was only necessary to be sure that a vehicle could safely pass through that high-altitude 
domain to lower altitudes where most of the deceleration and aerodynamic heating occurs. Neverthe- 
less, we thought we understood the main features of the high-altitude flight — it was out of chemical 
equilibrium, and ways were devised to study flow fields that were relaxing chemically at a finite rate. 
However, in this decade. Park (ref. 1) and others deduced that the flow field is out of equilibrium ther- 
modynamically as well. That is, the internal states of species are not equilibrated, and the gas may need 
to be characterized by more than one temperature locally (a translational temperature), as well as vibra- 
tional, rotational, and electronic temperatures. These temperatures equilibrate at finite rates that may 
need to be determined. Currently, there is intense activity aimed at achieving an understanding of ther- 
mochemical nonequilibrium because a new class of vehicles, hypervelocity aeromaneuvering vehicles, 
will function in that high-altitude regime. That regime is depicted by the "u"~shaped area centered 
about 9 km/sec in figure 10-1. An attempt to quantify the phenomena is shown by the solid lines labeled 
1 cm, 10 cm, and 1 m. These lines were inferred from figure 10 of reference 2. They represent experi- 
mental correlations of the distance behind a normal shock wave required for the flow to equilibrate 
thermochemically. For example, for the blunt-drag-configured, aeroassist flight-experiment vehicle 



km/sec 


Figure 10-1.- Flight regime highlighting thermochemical nonequilibrium. 
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being planned, the flow behind the bow shock tends to equilibrate in a distance of about 10 cm near peak 
heating (the total shock standoff distance is about 20 cm (ref. 3). For a lifting configuration, the flow 
field is more complicated. Nevertheless, an extensive literature on aeroassisted space transfer vehicles 
has been developed by the use of sophisticated approximate analysis. Such vehicles replace a costly 
near-Earth propulsive maneuver by an aerodynamic maneuver in the upper atmosphere. They are the 
key to economic utilization of Earth-moon space, the exploration of Mars, and other future space mis- 
sions. Nonequilibrium flow phenomena affect forces, moments, and heat transfer to vehicles. Some 
discussion of these vehicle concepts will be made in the next chapter. 

But in this chapter let us examine chemical nonequilibrium in a flowing gas in a simple way. 

The coupling of chemical rate equations to the flow-field equations, the solutions which result, and some 
useful notions will be illustrated. Importantly, all spatial variations in flow variables, thermodynamic 
quantities, and chemical composition occur because the chemical reactions proceed at a finite rate. 

The flow model chosen for the analysis is that of a normal shock moving at velocity U into qui- 
escent CO 2 at density p^. 1 It is like compressible flow in a constant- area duct (chap. 1) wherein a nor- 
mal shock wave alters the thermodynamic and chemical state of the gas — this time the latter changes at a 
finite rate. The corresponding flow field as seen by an observer traveling with the shock is shown in 
figure 10-2. 


co 2 

PocU 


CHEMICAL REACTIONS 

► — 

COUPLED TO FLOW 
P, U, Hj 

SHOCK 


Figure 10-2.- Flow of a compressible gas across a normal shock wave which leads to gaseous chemical 

reactions at a finite rate. 


FLOW-FIELD EQUATIONS AND BOUNDARY CONDITIONS 


The equations describing the flow field behind the shock are (neglecting transport phenomena): 


pu = PooU = constant 

(1) 


(2) 

u du + dh = 0 

(3) 


or, in the integrated form, 


'As of this writing, the current model of the Mars or Venus atmosphere is about 95% CO 2 and 5% N 2 (ref. 4). In 1962 the 
Mars atmosphere was considered to be less than 10% CO 2 while that for Venus was of the order of 10% CO 2 (ref. 5). 
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(4) 


u 

H = — + h = constant 

for the continuity, x-momentum, and energy equations, respectively, where x is the distance behind the 
shock. 


The enthalpy of the mixture of species is 


k 

h(p, p, ni, n2, . . nk) = £ njhi 

i=l 


(5) 


where nj is now the number of moles of species i per unit mass of fluid and hj is the enthalpy per 
mole of species i, or 



dT + hj 0 


( 6 ) 


where c pi and hi 0 are the specific heat at constant pressure and enthalpy of formation on a per-mole 
basis. 


The equation of state of the mixture is 


k 

P = PRT £ nj (7) 

i=l 

/^s 

where R is the universal gas constant per mole (1.98717x10-3 kcal/g mole K). 

The boundary conditions of the flow equations and the chemical-rate equation are specified 
immediately behind the shock and are at x = 0 


(for a strong shock (p^ « p s )), 


poo 

— = e 
Ps 

(8) 

u s = eU 2 

(9) 

pcoU^l-e) 

(10) 

— 9 2 

2 2 

(11) 
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ni = nu = T^~ , 112 = 113 = 114 = ns = n6 = n 7 = ng = 0 

s Mi 


(12) 


(that is, CO2 passes the shock without reacting), eight species are considered, and Mj is the molecular 
weight of species i. 

A A. 

For simple illustrative purposes, it is assumed that c p /R is constant, such that the combination 
of boundary conditions given in equations (9)-(12) and the equation of state, equation (7), yields 



■ 

A \ 


8 = 

2 

C P1 

i R i 

-1 





l-l 


(13) 


This completes the boundary conditions behind the shock wave. 


The Modified Flow-Field Equations 

It is convenient to rearrange the flow-field equations as follows. Equation (5) can be written in 
the differential form 


„ 5h , dh , v ^h . 

dh 'aF dp+ ^ dp+ ^5^ d ” i 

which combined with equations (l)-(3) and (5) expressed as equation (6) yields 

0h 


V on , 
2,;w dn ' 


du p0h/3p) fi 3n i 


1 — u^ 


Qh/dp)-(l/p) 


-(3hy9p) 


which can be written as 


du 

dx 


X "i[(c P /R) - !] 

{ PU 2 i=l 

P X 

i=l 


2 


k Cpidni 


i=l 


dx 


i=i R 


dx 


k Hn- ^ 
th ani 

^ hi ^ 

i=i 


\ 


k k r k Cn 

X«4 RTXni^ 

w i=i R i=i R 


(14) 


(15) 


(16) 


where dn;/dx will be expressed subsequently (eq. 28). 

From kinetic theory, the ratio of specific heat to the universal gas constant can be written 
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(17) 


Cp _ Cy i _ D + 2 

R R 2 

where D is the number of degrees of freedom of a molecule of species i. For monatomic gases D = 3 
so that Cr/R = 5/2; for diatomic gases D = 5 so that Cp/R = 7/2; for triatomic gases D = 7 so that 
ZpfR = 9/2. 

A XV 

From statistical mechanics, full vibrational excitation adds 1 to the ratio Cp/R for diatomic 
species (refs. 6 and 7). For illustrative purposes it is assumed that c p /R is fully excited vibrationally. 
Moreover, for the linear triatomic molecule its value is 7 according to reference 8. Electron excitation is 
not considered here. Its effect on the thermodynamics of the problem is discussed in reference 9, and 
will be mentioned subsequently. Its effect on reaction rates is currently being assessed, but is beyond the 
scope of our present considerations. Thus for full vibrational excitation (where the error in h; at 5000 K 
is from reference 8: 


Species, i CO 2 

CO 

O 2 

O 

C 

A 

Cp: 

7 

R 

9/2 

9/2 

5/2 

5/2 

0-8% 

at 5000 K 

3% 

1% 

2% 

— 


The set of differential equations with their boundary conditions describing the interrelated 
chemical and flow processes in that part of the flow field behind the shock that is not in chemical equi- 
librium is thus complete, and we turn our attention to the chemical reactions themselves and the rates at 
which they proceed. 


Chemical Reactions 

Behind the shock wave, the reactions assumed for the dissociation, ionization, and recombination 
of CO 2 and its components are (terminology is identified and discussed below) 
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r 


( 1 ) 


k fl 

1 CO2 + M<±C0 + 0 + M 

kbi 

( 2 ) 

2 C0 + M<±C + 0 + M 

(5) 

3 CO + M«±CO + + e- + M 

(3) (6) 

4 0 + M*±0 + + e _ + M 


(18) 


(19) 

( 20 ) 
( 21 ) 


(4) (7) (8) 

5 C + M<iC + + e" + M ^22' 

where species 1 through 8 are identified by (1) through (8); r is the reaction number, kf r and kb, are its 
forward and backward reaction rate coefficients, respectively; the ratio kf/kb,. is the equilibrium coeffi- 
cient Kc r ; a ri and ft-j are the stoichiometric coefficients for species i for the reactants and products of 
reaction r, respectively (unity everywhere here); and M is any collision partner. Although it may be 
expected that electrons are more efficient than the large particles in producing ionizing reactions, no 
attempt is made to distinguish between electrons and larger particles as collision partners 2 and 


n M = n (23) 

The rate at which reactions in equations (18)-(22) proceed to the right is characterized by the forward 
reaction rate coefficient expressed either in the form given by collision theory 


k fr = 


2PAd 2 


a(s - 1)! 


- V ' 2 

2kRT I 

(n \ s 
E fr 

M* 1 | 

RT 

/ 


:-l 


e -(E f /RT) 


(24) 


where M* is the reduced molecular weight corresponding to a collision partner M, and is 
(1/Mi + 1/Mm)' 1 ; P is the sterric factor; Ef r is the activation energy for the forward reaction; A is 
Avagadro s number; d is the average diameter of the colliding pair; s is the number of classical squared 
terms of energy contributing to the reaction; and a is a constant. The corresponding Arrhenius form of 
the reaction rate coefficient is 


For argon, it is generally agreed that electron collision partners are important to the rate at which equilibrium is approached 
behind shock waves, although there is little agreement as to the mechanism for initiating ionization (refs. 10-14). 

For air at shock speeds up to 10 km/scc, ionizing processes listed in order of decreasing importance are atom-atom 
collisions, photoionization, electron impact, atom -molecule, and molecule-molecule collisions according to reference 15. 
However, the importance of electron impact increases with shock speed. 

For carbon dioxide, although the ionization rate coefficients with electron collision partners are conceivably an order of 
magnitude larger than those with large particle collision partners, we do not know either within an order of magnitude and 
thus do not attempt to distinguish between them. 
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kf r = Bf r T c ^e*<* T (25) 

The backward reaction rate is included in the equilibrium coefficient 3 * * * * * expressed in either of two forms 


k fr 


K*-?-pII S = Bc,rt. e^T) 

^br 1 


i=l 


(26) 


where the bar refers to equilibrium values, E Cf refers to energy, and is the backward reaction rate 
coefficient. Symbolically for reaction r. 


X arjXi <± 

i=l 


8 


IWi 


(27) 


Thus equations (l)-(6) are combined with equation (14) to yield equation (15), which can be 
wntten in the form of equation (16), for which the chemical rate equation for the rth reaction is needed. 
By use of equation (27) this can be formed as (refs. 6 and 16) 


dnj 

dx 


— X(Pr,i-ar,i)k f] 

P u r=l 


"8 s 

II (Pni) 01 *-^-!! (P 11 !)^ 


i=I 


cr i=l 


(28) 


which is used to solve equation (16). 

If there are m types of atoms in addition to free electrons, (m + 1) of equation (28) can be 
replaced by m statements of conservation of atoms (in this case carbon and oxygen) of the form (where 
a Z j is the number of atoms z in species i) 


8 8 

X W = X Wh (29) 

i=l i=l 

plus one statement that the number of ions equals the number of electrons (only singly ionized species 
are considered). In applying equation (29), we consider an atomic ion to be equivalent to one atom of 
the appropriate species. It is noted that equation (29) can be differentiated to provide m values of 
dnj/dx for use in equation (16). 


3 The backward reaction rate is a misnomer for three body collisional backward reactions, especially in reaction 1, where the 

reaction does not proceed in the backward direction in a straightforward manner (ref. 17). Nevertheless, Fay (ref. 18) uses an 

equivalent reaction, 2CO + O 2 -> 2CC>2, while reference 19 uses the reaction CO + H 2 O -> CO 2 + H 2 to oxidize CO. But it 

appears that the reaction for oxidizing CO is a branched chain reaction which is not yet defined (ref. 17). Fortunately for our 

problem, the backward oxidation reaction for CO is not important, and our results are not affected significantly by our 

approach. 
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This completes the set of equations used to solve the nonequilibrium flow behind the shock 
wave. It is instructive to have the equilibrium conditions which the nonequilibrium properties eventually 
approach some distance behind the shock wave. 


Equilibrium Conditions 


The 13 equilibriurnproperties of interest (designated by a bar) include the flow velocity u; ther- 
modynamic properties p, p, and T; and chemical concentrations n, ni . . . ns. They are obtained from 
the simultaneous solution of the following 13 algebraic equations. By definition 

8 

n = X"i (30) 

i=l 

Statements of conservation of oxygen and carbon atoms are 

n3 = 2ni s -2ni -n2-n5-n6 (31) 

and 

h4 = ni s -ni -n2-n5-n7 (32) 


The number of moles of electrons equals the number of moles of ions and therefore 

ng = ns + n6 + nj (33) 


The equation of state is 


E = nRT 
P 


The strong shock relations give 


p = pooU(U -u) 
8 


ni s h? + y(U 2 -u 2 ) = X"i 


i=l 


/- \ 
c Pi 


l\ r / 


RT + hj 


_ pooU 

u = — — 

P 


(34) 


(35) 

(36) 

(37) 


The remaining five equations will not be listed, but are simply equation (26) written for the five reac- 
tions given in equations (18)-(22). 



Method of Solution 


Solutions for the nonequilibrium flow field were obtained by numerical integration of the 
differential equations, equations (2), (3), (16), (28), subject to boundary conditions, equations (9)-( 1 2), 
making use of equations (1), (7), and (13). The integration was performed numerically, making use of 
the Adams-Moulton predictor-corrector variable-step integration scheme (ref. 20). Values of dissocia- 
tion energy and some of the other physical constants employed in the analysis were either obtained or 
estimated from information in references 21 and 22. The rate and equilibrium coefficients will be pre- 
sented subsequently. 

The set of 13 equilibrium equations (eqs. (30)-(34) plus the five equations, eq. (26), in which r 
takes on the value 1-5), some of which are nonlinear and transcendental, were solved in a direct way by 
the Newton-Raphson method (ref. 23, p. 213). In addition, the results were verified by the method of 
tracing the locus of roots described in reference 24. 


Chemical Kinetics 

Reaction rate coefficients - The estimated values of the constants in the forward reaction rate 
coefficients were obtained with the aid of collision theory and are listed in table 10-1 in the Arrhenius 
form (eq. (25)). As in references 9, 25, and 26, it has been assumed for present purposes that P = s = 1 
in the forward reaction rate coefficient expression (eq. (24)), and that the activation energy equals the 
reaction energy. An equivalent assumption is that collision of pairs having total energy in a specific 
degree of freedom (such as translation) equal to or greater than the reaction energy will result in a reac- 
tion (ref. 16). The neglected effect that reaction rates tend to be lowered because not every such colli- 
sion results in a reaction (i.e., P is actually less than unity) is compensated to some extent by the 
neglected effect that reaction rates tend to be raised because more than one degree of freedom may par- 
ticipate in the reaction (i.e., s is actually greater than 1). For these reasons (which are discussed more 
fully in ref. 25) and for lack of experimental evidence, we estimate reaction rates in accord with the 
above assumption, and the result appears in figure 10-3. 


TABLE 10-1.- REACTION RATE AND EQUILIBRIUM COEFFICIENT DATA (ref. 9) 


r 

Reaction 

cm 3 

«fr 

Efr 

kcal 

»c r , 
g mol 

«cr 

Ec r 
k cal 
g mol 

g mol sec K 0 ^ 

g mol 

cm 3 K ttc i 

i 

co 2 *±co + o 

6.955xl0 12 

0.5 

125.75 

0.17141xl0l 3 

-2.6294 

135.73 

2 

CO C + O 

7.238xl0 12 

.5 

256.17 

.17713x107 

-1.1106 

268.57 

3 

CO CO + + e- 

7.238xl0 12 

.5 

323.18 

.41958X10- 4 

.72788 

333.44 

4 

0 0 + + e - 

7.344xl0 12 

.5 

314.05 

.42160xl0- 10 

1.9483 

307.73 

5 

C C + + e~ 

7.895x1012 

.5 

259.84 

.39181x10-6 

1.0415 

i 

264.04 
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Figure 10-3.- Reaction rate coefficient estimates (ref. 9). 

At a temperature of about 20,000 K, the rate coefficients of the various reactions exhibit rather 
large overall differences (~10 2 for example). Since the colliding pairs (electrons excepted) do not differ 
greatly in size or molecular weight, the large spread in the rate coefficients at low temperatures is 
attributable to differences in the activation energy Ef r which appears a^part of the exponent in equa- 
tion (25). At very high temperatures, however, the entire exponent Ef/RT is small so that differences 
in activation energy are of diminishing importance and kf r ~ T 0 ^ (~T^ in accord with the assumption 
s = 1). Thus at a very high temperature, the various reaction rate coefficient estimates do not differ 
greatly from one another. 

Equilibrium estimates - The equilibrium coefficients were obtained by empirically fitting data 
obtained from references 21 and 22 by use of the far right side of equation (26). The resulting values of 
the parameters in equation (26) are shown in table 10-1, and the equilibrium coefficients themselves are 
shown in figure 10-4. 


Chemically Relaxing Flow-Field Results 

Profiles of thermodynamic, kinematic, and chemical quantities in the nonequilibrium part of the 
flow field behind the shock were obtained by integration of equations (2), (3), (16), and (28) subject to 
their boundary conditions and with the use of the chemical rate coefficient estimates. The solution of the 
flow-field equations is fully coupled to the chemical relaxation equations. 
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Figure 10-4.- Equilibrium coefficients (ref. 9). 

An illustrative solution for a high shock speed and specified ambient density (15 km/sec and 
Poo/po = 10 -4 , respectively) is shown in figure 10-5. In part (a), all quantities have been made 
dimensionless with respect to their values immediately behind the shock wave. As chemical relaxation 
proceeds (increasing x), temperature diminishes monotonically. On the other hand, particle velocity 
first increases and then decreases while density varies in the opposite way with increasing distance 
behind the shock wave. This behavior has been discussed in reference 25. The variation of pressure and 
enthalpy in the flow field is very small, as would be expected from the appendix on flow-field sensitivity 
estimates which follows. 

The mole fraction profiles of the eight chemical species are shown in figure 10-5(b). It is seen 
that CO 2 vanishes a short distance behind the shock wave. Both CO and O rapidly increase in concen- 
tration for a short distance behind the shock wave and then diminish because of further dissociation and 
ionization processes. Electrons are the predominant species in the flow field for x greater than half a 
millimeter. 4 Although the nonequilibrium electron concentration has overshot the equilibrium value, it 
does approach that value when x is large. This behavior is observed in other examples and in refer- 
ences 9, 25, and 26 as well. It is interesting to see that the concentration of atomic oxygen remains high 
in much of the relaxation region. 


4 It should be noted that the electron concentration in the nonequilibrium profdes to be presented is probably slightly lower 
than it should be because doubly ionizing reactions have been neglected. However, for the ambient conditions considered in 
this analysis, doubly ionized species would tend to disappear as equilibrium is approached. 
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Figure 10-5.- Nonequilibrium profiles; U = 15 km/sec, poo/p 0 = 10 -4 (ref. 9). 

Profiles for a lower shock speed (9 km/sec) at the same low ambient density (poo/p 0 = 10 -4 ) are 
shown in figure 10-6. It is noted that the density is far from its equilibrium value 10 cm behind the 
shock. The dominant species in most of the relaxation region is atomic oxygen. Interestingly, it 
achieves an almost constant value only 3 mm behind the shock in spite of the fact that CO is still disso- 
ciating into C and O. Thus the ionization of atomic oxygen keeps pace with the production of atomic 
oxygen throughout most of the relaxation region. The concentration of electrons and ions is relatively 
small at this lower shock speed. This figure may be compared with figure 5 of reference 26, in which 
ionization reactions were neglected. The thermodynamic and flow-field structure is much the same in 
either case, although the overall concentration of the nonionized species is higher in that reference, as 
would be expected. 
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Figure 10-6.- Nonequilibrium profiles; U = 9 km/sec, Poo /p 0 = 1 (H (ref. 9). (a) Thermodynamic and 

flow properties, (b) Species concentration. 


(f ^ of / chemically relaxin S field at much lower velocity is shown in figure 10-7 

( rom ref. 26). The first part shows the velocity and thermodynamic variation, while the second part 
!„; S n J e B chemical species concentrations which correspond. The chemical relaxation distance for this 
, .f , c J eds °. cm ‘ Again, pressure and static enthalpy were almost invariant, while temperature, 
density, and flow velocity varied significantly. The variation of mole fraction of CO and O were identi- 
ca because there is no ionization, and the equilibrium value of CO 2 was essentially zero. The initial 

rev^H f th enSlt) I and veIocity b ^ h [ nd the shock tended to depart further from equilibrium before they 
reversed themselves to approach the equilibrium values. The conditions which govern the initial slopes 

° nd ^ vanables ^ denved in reference 25. Additional solutions are shown in both references 9 
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Figure 10-7 - Nonequilibrium flow-field profiles; U = 6 km/sec, p<,o/po = 10" 4 (ref. 25). (a) Velocity 
and thermodynamic properties, (b) Species concentration. 


Relaxation Distance 

Results of a large number of examples can be summarized in terms of a chemical relaxation dis- 
tance; the distance behind the shock wave at which the overall chemical relaxation has gone a given 
fraction (or percentage) toward completion. The fraction of reaction completion is defined in terms of 
the total change in the number of moles per unit mass occurring in the distance x divided by the total 
change in number of moles per unit mass required to achieve chemical equilibrium and is 
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In all of the examples that were studied, n, and thus rj, were monotonic throughout the relaxation region 
in spite of the fact that n, is not monotonic for several species. 

The relaxation distances for rj = 0.5, 0.8, and 0.9 were obtained. The distance for 50% reaction 
completion is quite short, not more than 10 -1 cm for the entire range of shock conditions studied (refs. 9 
and 25). The distance for 80% reaction completion is roughly a decade higher than for 50% completion. 
Finally, the relaxation distance for 90% reaction completion is quite large, ranging from IQ" 2 to 10 2 cm 
over the entire range of shock conditions, as shown in figure 10-8. 





Figure 10-8 - Relaxation distance for 90% reaction completion (ref. 9). 

It is worth mentioning that these relaxation distances are not very sensitive to the reaction rate 
assumptions. It was noted in reference 25 that a two to three order of magnitude change in reaction rates 
results in only a one order of magnitude change in relaxation distance. 


209 



CONCLUDING REMARKS 


The interdependent dissociative and ionizing chemical rate processes behind shock waves in car- 
bon dioxide have been examined for shock speeds from 6 to 20 km/sec at ambient densities of 10 -2 to 
10 -4 ratioed to standard Earth atmosphere density. Because of a lack of reaction rate data for carbon 
dioxide and its components, reaction rate coefficients have been estimated with the aid of some simpli- 
fying assumptions and collision theory. The effects of chemical relaxation on the flow field behind 
normal shock waves have been studied by use of the coupled fluid flow and chemical rate equations. 

Results have been presented for both the equilibrium and nonequilibrium regimes behind the 
shock wave. At the higher speeds (15 km/sec), electrons are the predominant species in most of the 
nonequilibrium chemical relaxation region. Atomic oxygen is the predominant species in the lower- 
speed examples and is also present in high concentration at the higher speeds. 

Equilibrium concentration of electrons is presented here in terms of mole fraction of the mixture 
and in terms of number of electrons per unit volume as in reference 9. The equilibrium concentration 
varies from about 5xl0 14 to 5xl0 18 electrons/cm 3 for a shock speed of 9 km/sec at a density ratio of 
10 -4 standard Earth atmosphere density to 20 km/sec at a density ratio of 10~ 2 . Results are also 
presented in the form of a chemical relaxation distance behind the shock wave. 

The nonequilibrium results obtained are considered to be preliminary estimates for a number of 
reasons. They are based on estimates of chemical reaction rate coefficients. Indeed, it is not certain that 
the reactions studied are the correct reactions. Not all of the conceivable reactions have been included in 
the analysis. The equilibrium results are in a somewhat better situation in that they are independent of 
any estimates of reaction rate coefficients. They are not influenced significantly by internal excitation 
phenomena for the high shock speeds considered. 

Interestingly enough, the study of reference 9 shows that the influence of the state of internal 
excitation on most of the chemically relaxing flow field is small insofar as this influence is manifested 
through the thermodynamics of the problem. However, it is possible that dissociation and ionization 
processes may proceed more readily from excited states, in which case the internal state of excitation 
would be important in the entire chemical relaxation region. 

Reference 9 shows that there is a small part of the nonequilibrium flow field just behind the 
shock wave in which the state of internal excitation of the various species has a significant influence on 
the flow properties. This region is likely to be important to nonequilibrium radiation problems; thus a 
knowledge of both the chemical and internal relaxation is important at high temperatures. 

Knowledge of the behavior of these common chemical substances at elevated temperature is 
meager in spite of the fact that they have long been important species in many combustion problems. 

The behavior of these substances has acquired new significance relating to planetary entry problems. 
There is much to be learned both from the macroscopic and the subatomic points of view. 


APPENDIX 


FLOW-FIELD SENSITIVITY ESTIMATES 


In the strict sense, all of the flow-field quantities p, u, p, h, and T and species concentration 
vary in the reacting flow field behind the shock. Some of these quantities, however, will vary more than 
others. The purposes of the following discussion is to show the relative variation of the quantities p, u, 
p, and h in the flow field behind the shock. 

At some distance x behind the shock, the ratio of local density to the density immediately 
behind the shock is defined as g. 



(Al) 


Thus the fractional change in p from p s is 

p - ps _ Ap • j 
Ps Ps 

The corresponding change in u is from equation (1) 

Au = J__ 1 
u s g 


(A 2) 


(A3) 


Similarly, from equation (2) the change in p (noting that p s = pooU 2 = p s u s U) is 



(A4) 


where e = p<Jp s . Finally, the change in h if we assume that hoo«h s andU 2 » u s and make use of the 
integrated form of equation (3) is approximately 


Ah 

hs 
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Assuming for exemplary purposes that g may be as high as 2, and that e * 10" 1 , the corre- 
sponding changes in density, velocity, pressure, and enthalpy (from eqs. (A2)-(A5)) are 100, 50, 5, and 
0.75%, respectively. Similarly, if g is as small as 1/2, we arrive at the same qualitative result. Thus we 
conclude that pressure and enthalpy do not vary greatly in the chemical relaxation region, although den- 
sity and particle velocity vary a great deal. 


It is useful to make these estimates because they highlight the main features of the results before 
we have them. But more importantly, they sort out the most sensitive variables; temperature, density, 
and fluid velocity which are affected significantly by chemical relaxation, but not static enthalpy and 
pressure. 'Diese notions are helpful in devising numerical schemes to compute chemically relaxing 
flows efficiently. Thus the insensitivity of pressure means that numerical methods need not be tailored 
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to allow for large pressure excursions. Making sensitivity estimates is a good habit to acquire. Such 
estimates in advance should influence the methodology, and may determine whether a computed solu- 
tion is obtained at all. 
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CHAPTER 11 


HYPERVELOCITY AEROMANEUVERING FLIGHT AT VERY HIGH ALTITUDES 


In recent years, interest has developed in the use of vehicles that can transport payloads from one 
orbit to another. Such vehicles could be space-based. The capability of such vehicles is greatly 
enhanced if propulsive maneuvers near Earth can be replaced in pan or almost entirely by a shallow dip 
into the upper atmosphere where the maneuver is performed aerodynamically. Such vehicles are called 
aeroassisted space transfer vehicles (ASTV), or by some, aeroassisted orbital transfer vehicles (AOTV). 


Thus a payload returned by these vehicles from a distant orbit (geosynchronous orbit (GEO)) for 
example, to a low Earth orbit (LEO) would experience an aerodynamic drag maneuver in the upper ’ 
atmosphere before skipping out to rendezvous with another object in a low Earth orbit. Or the ASTV 
could be returning at hypervelocity from an orbit considerably farther than geosynchronous; from a 
ubration center (a gravitational null point in the Earth/moon system) which is at lunar distance from 
barth or from a location outside the Earth/moon system altogether. Payloads may be enhanced greatly 
(doubled or even more) by the use of the drag aeromaneuver in the upper atmosphere rather than by a 
propulsive deceleration near Earth. 3 


„ A S t C °” d ^ ind of aeromaneuver would transfer a payload from one orbit to a second orbit at a 
different orbital-plane inclination angle. The corresponding maneuver would involve a lateral banked 

turn which is performed using aerodynamic lift. Or there may be a combined lift and drag maneuver in 
the upper atmosphere. 

In a primarily drag-type aeromaneuver, the vehicle configuration would be large and blunt like 
an umbrella. Because of the blunt-bow shock wave, the heating to the vehicle would be radiative to a 
large extent But the turning aeromaneuver requires a lifting capability, and would be more of a slender- 

^ C01 J eSp0ndm 5 ly ’ aerod y namic heating would tend to be high on the slender nose 

and on the leading edges of wings and control surfaces. It would be primarily convective heating. 

a . A /V h ;; d vehi f! e concept studied in recent years is ground-based, would fly into orbit (perhaps) 
and eventually would return to land on Earth. It, too, would perform hypervelocity aeromaneuver? 
which would have to be done at very high altitude because of the heating problem.^ The primary heating 
Site ° n thC aSCent ’ bUt thC technology of the liftin S ASTV would be relevant to this cSiceJut high 


This family of vehicle concepts requires extensive research to form the broad technology base of 
a new generation of aerospace vehicles. It is necessary to understand real gas flow fields, aerodynamic 
ating, improved reusable materials and structures, thermal protection systems, efficient aeromaneu- 
vering strategies, and guidance and control problems; and if the vehicle requires an air-breathing pro- 

ESS ££*££££!& S“ esradon prob,ems need ‘° “■"* 

Insofar as the ASTVs are concerned, some features of the flight environment can be assessed bv 

EJSH 7 cu . rrcn ^ f echn °logy base. Consider, for example, the atmospheric drag maneTve^y a 
blunt vehicle wherein radiative transfer is a dominant heating mechanism. We can infer the kinds of 

SS L fr * d *T d - Fig r,' M ^ ^ 1) depicts the “a 

vehicle from GEO, and from a distance of five times GEO (5xGEO)— which is about 40% of the 
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distance to the moon. The figure would be essentially the same if the lunar distance return were shown. 
The ordinate is a log scale of altitude in meters, and the abscissa is the inertial speed relative to a rotating 
atmosphere in kilometers per second. The vehicle leaves GEO (about 35,000-km altitude) and arrives at 
the "top" of the atmosphere at an inertial speed of about 10 km/sec, decelerates by atmospheric drag to 
low Earth orbital speed of about 7 km/sec, exits the atmosphere, and rendezvous with the Space Shuttle 
(for example). The perigee in the atmosphere is at an altitude of about 80 km. Entry from 5xGEO or 
lunar return would be about the same because Earth escape speed is about 12 km/sec. 

A Space Shuttle (STS) entry is shown for reference. Since the entry speed of an STS is about 
7 to 8 km/sec, we see that the ASTV entry is about twice as energetic as the Shuttle entry (energy varies 



Figure 1 1-1- Aeroassisted orbital transfer from distant orbits. 
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with the square of velocity). Moreover, since the ASTV decelerates at a much higher altitude than the 
Shuttle, the ASTV resides in a very energetic rarefied environment. 

The two horizontal lines called "quasi-continuum limit" pertain to altitudes below which contin- 
uum theory can be applied with minor modifications (ref. 2), especially to large blunt objects. For the 
10-m nose radius, that altitude is 1 12 km. Below an altitude of 103 km, continuum theory may be 
applied without modification. Most of the deceleration and heating occurs in the full continuum region 
because the air molecules need to be present to develop significant aerodynamic effects. 

Some of the real gas phenomena that should be considered are shown in the figure. An important 
family of lines of constant free-stream energy flux (p„U 3 /2) in watts per square centimeter are shown on 
the right. They vary inversely with altitude in an exponential manner. They serve as an index to the 
severity of aerodynamic heating and of thermal-protection requirements. Experience has shown (ref. 3) 
that maximum heating can be as much as 10% of the free-stream energy flux. Thus, if a thermal- 
protection material can withstand 50 W/cm 2 without the need to be refurbished, these lines indicate that 
flight below an 80-km altitude is not advisable unless there is a significant improvement in thermal- 
protection materials. Detailed computations of heating rates indicate that these estimates are about 
correct. 


Another important line on the figure is labeled the "ionization" threshold. Many mechanisms by 
which the hot gas cap can emit thermal radiation exist to the right of the line. Indeed, thermal radiation 
is the dominant heating mechanism for the blunt drag configuration. Moreover, the radiation is from air 
that is not in thermal or chemical equilibrium, as will be shown. Although vehicles have flown through 
this region routinely, the peak heating was at lower altitudes in a regime where the air is in equilibrium 
behind the bow shock wave. Thus the total contribution of nonequilibrium radiation to Apollo heating 
was about 1% and to the Space Shuttle about 5%; but for the drag-configured ASTV it is estimated to be 
about 50% (or more). Park (refs. 4 and 5) has explored ways to model the physical-chemical phenom- 
ena that produce nonequilibrium radiation. 

Aeroassisted space transfer vehicle problems have been studied by approximate but highly 
sophisticated methods. Configurations for both high drag (lift/drag ratios of 0.3) and high lift (LTD ~ 2) 
have been defined that are stable, controllable, and of minimum weight (refs. 6, 7, and 8). Two config- 
urations are shown schematically in figure 11-2. Each configuration can be placed into low Earth orbit 
in a single Shuttle payload. Note that the drag configuration is large (27 m in diameter), and is very 
blunt. It is not symmetrical, so that lateral aerodynamic forces can be generated to maneuver the vehicle 
aerodynamically to some extent, to make trajectory corrections. 

The lifting surface shown on the bottom of the figure is an early version (Davies and Park 
(ref. 8)). A payload would be placed in the dead-air region, along with some propulsive capability. 

The development of significant lift at high altitude is very difficult, so the lifting surface must be 
very carefully designed to accomplish that. Although simple Newtonian theory provided the basic shape 
shown, it was refined by applying rarefied viscous forces locally and integrating them over the entire 
surface. The configuration shown was the best of several configurations studied, which all had a com- 
mon dimension and payload. This shape was the only one which produced the necessary aerodynamic 
stabilizing moments. Thus both a high-drag and a high-lift configuration appear to be the best candi- 
dates for ASTVs. No real use for an intermediate configuration was found. A combination of moderate 
drag and moderate lift does not appear to do anything well. Admittedly, these configurations appear 
somewhat ugly, but in time they acquire an undefinable appeal. 
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(b) 

Figure 1 1-2.- ASTV configurations, (a) High-drag configuration, (b) High-lift surface. 

Both structural and thermal-protection weights have been estimated for these configurations 
(refs. 9 and 10). Moreover, clever flight strategies have been devised and mission capabilities were per- 
formed. A few of the round-trip mission payloads for the drag configuration are shown in figure 1 1-3. 
Although GEO is in the equatorial plane, other orbits are inclined to that, so that both an aerodynamic- 
drag maneuver and a propulsive plane-change maneuver (far from Earth) were involved in the analysis. 
Delivery (open bar), retrieval (shaded bar), and combined delivery-retrieval (solid bar) payloads for 
round-trip missions from low Earth orbit (Shuttle-compatible) were examined. It is noteworthy that the 
retrieval missions have a large payload capability which nonintuitively increases with distance in space. 
Although it is not shown, about 23 tons can be retrieved from a libration center (at lunar distance) by a 
vehicle which leaves the Space Shuttle weighing 30 tons (the Shuttle payload capacity). The retrieval 
time is less than 2 days. It may be advantageous to store space assets in space for retrieval in an emer- 
gency — rather than to store them on the ground where they cannot be deployed, even in bad weather. 
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MAXIMUM PAYLOAD CAPABILITY FOR DELIVERY. RETRIEVAL, 
AND COMBINED DELIVERY RETRIEVAL ROUND-TRIP SORTIE 
MISSIONS WITH 480-sec SPECIFIC IMPULSE ROCKET ENGINE 



Figure 1 1-3 - ASTV mission payload capability. 

The many mission analyses revealed a number of nonintuitive advantageous results. Reduction 
of aerodynamic loads and heating by use of multiple passes through the upper atmosphere has also been 
examined. Excursions to other planets have been studied as well (refs. 1 1 and 12). The simultaneous 
use of propulsion and aeromaneuvers has been considered (synergetic maneuvers). 

For the drag configuration, an Aeroassist Flight Experiment (AFE) is being planned whereby a 
highly instrumented one-third-scale model will be launched by the Space Shuttle, driven into the atmo- 
sphere propulsively, perform the drag aeromaneuver, and skip out of the atmosphere to be retrieved by 
the Space Shuttle. Figure 1 1-4 shows two candidate ASTV configurations for AFE. The configuration 
on the right is the NASA Johnson Space Center (JSC) shape partially modified at the extreme diameter. 
As shown at the bottom of the figure, models of these configurations were flown in the Ames Ballistic 
Range Facility (by P. Intrieri). The resulting shadowgraphs show three things. First, the models are 
stable (they fly forward without tumbling and tend to seek a fixed orientation). Second, the shock waves 
and flow-field structure are visible to some extent. Third, the windows of the facility need to be 
replaced. 

It appears likely that the modified JSC configuration will be chosen for the flight test. The 
radiometer experiment to be carried aboard the flight test vehicle is of special interest. The principle 
investigator for the experiment is W. C. Davy of NASA Ames Research Center (ref. 13). Radiative 
emission from the hot gas cap is especially interesting in the context of this monograph. Figure 11-5 
shows a diagram of an ASTV with some shock- wave structure. The left side of the figure shows sche- 
matic profiles of temperature, chemical species, and radiative emission along the forward stagnation 
streamline. Immediately behind the bow shock, the gas is out of thermal equilibrium, and three tem- 
peratures are shown to characterize the gas locally; the translational (or large-particle) temperature, a 
vibrational temperature for diatomic species, and an electron temperature. If there are many diatomic 
species (O 2 , N 2 , NO, and their ions), many vibrational temperatures may need to be used. Behind the 
bow shock, air composition is atmospheric N 2 and O 2 . These dissociate and ionize as flow proceeds 
toward the body until they (like temperature) approach a thermochemical equilibrium state. The 
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Figure 1 1-4 - AFE candidate configurations. 
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ENERGETIC FLIGHT AT HIGH ALTITUDE 
ASTV REAL GAS PHENOMENA 



BOW SHOCK BODY 



Figure 1 1-5.- Thermochemical nonequilibrium flow-field schematic. 

corresponding radiative emission "overshoots" because of excessive temperature, excited chemical 
species, and a nonequilibnum chemical composition. Again, as the flow proceeds toward the body the 
radiation tends toward a thermochemical equilibrium value. The extent of the thermochemical nonequi- 
libnum effect is estimated in figure 1 1-6. The figure highlights the region denoted by the short double 
horizontal lines at 50 km altitude that was shown in figure 1-1. The three curved lines labeled 1 cm 
0 cm, and 11 m indicate the distance behind a normal shock wave required for the thermochemical state 
to approach within 10% of the equilibrium value. The estimates were made by use of shock-tube 
experiments by Allen, Rose, and Camm (ref. 14) correlated in figure 10 of that reference. If these 
estimates are related to a large blunt, drag-configured ASTV having a shock standoff distance of about 
1 m about 10 cm of the flow behind the shock is not in thennochemical equilibrium at peak heating at 
an altitude of 80 km. Flight in this domain is not well understood. The AFE flight experiment is 
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Figure 11-6 - Comparison of vehicle flight regimes in Earth's atmosphere, highlighting thermochemical 

nonequilibrium. 


expected to increase understanding of hypervelocity flight at very high altitudes, and the forces and 
heating that occur because of real gas effects. 


Finally, with respect to hypervelocity lifting configurations, flight at very high altitudes reduces 
convective heating (to the leading edges of lifting and control surfaces because of rarefied gas effects), 
which leads to concepts that can result in vehicles that do not have to be refurbished after each flight. 
But although aerodynamic heating is diminished at high altitude, so is lifting capability. Figure 11-7 
(private communication, M. E. Tauber, Ames Research Center, 1986) shows the degradation of the 
lift/drag ratio above about 80 km (the perigee determined by materials that can withstand the heating 
rates) for the Space Shuttle. The lifting ASTV shows some promise of improvement over the 



Figure 1 1-7.- Space Shuttle maximum lift/drag ratio variation (approximate). 


performance shown, although not a great improvement. Clearly, there is a tradeoff between tolerable 
heating and acceptable lift. At this writing, that tradeoff appears to allow up to a 40° low Earth orbital 
plane change capability by aeromaneuvering. This is a remarkable capability — a propulsive capability 
of 60° orbital plane change requires about as much propulsion as it takes to launch from Earth to low 
Earth orbit. 


For our purposes, we see that the real gas aerodynamic environment at hypervelocity opens a 
new era of flow-field research. Mach number is relatively unimportant; it deals primarily with com- 
pressibility effects of the gas. Velocity is very important; it deals with the energetics of the gas — 
thermal, chemical, and physical effects associated with real gas behavior. These effects are only 
partially understood, and at very high altitudes these effects become increasingly complicated. 


Thus hypervelocity flight, especially at high altitude, requires both computational and experi- 
mental research. Computation must include more than configuration, momentum, simple energy rela- 
tions, and an equation of state. It must include these and much more. We might envision the com- 
plexity of solving a set of state, momentum, energy, diffusion, electric-charge conservation, finite-rate 
chemistry, thermal relaxation, and electrodynamic equations that are nonlinear and coupled to one 
another. We must acquire the basic physical and chemical inputs to these equations, including transport 
phenomena; radiative properties; and physical, chemical, and thermal rate data. The development of 
computational codes may proceed as shown in figure 11-8. The three columns are intended to indicate a 
melding of computational fluid dynamic techniques, real gas aerothermodynamics, and physical and 
chemical inputs. 


Moreover, hypervelocity simulation experimentally requires high-enthalpy research facilities and 
large power supplies. Heating test air in a reservoir by chemical means is not adequate. Gas flow 
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Figure 1 1-8.- Computational aerothermodynamic code development. 
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through a high enthalpy facility must be characterized and understood in detail. Moreover, computa- 
tional developments must be validated by experiment (as well as flight test — fig. 1 1-8). 

In proceeding into these poorly understood areas, we must be very thoughtful. A continual 
assessment of where we stand both computationally and experimentally is essential. Computationally, 
we must order the terms in the equations of change for configurations and flight conditions that are 
appropriate. We must reexamine transport and transfer phenomena and realize that both are affected by 
thermochemical rate processes. We must assess what problems are crucial to solve (some problems may 
be interesting, but unimportant) and make sure there is a correspondence among computation, ground- 
based tests, and flight tests. There must be a correspondence between codes and experiments — or we 
may be attempting to validate codes that are not valid. Experimentally we must carefully consider the 
limits of ground-based experiments. A test in a wind tunnel in which the test gas is near liquifaction 
may be importantly different from flight at the same Mach number. Problems abound. But we are on 
the threshold. The possibilities of space utilization for science, commerce, and national and interna- 
tional goals are enormous. 
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